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Page  B. 

Here  carried  out  100  years  frca  the  birthday  of  Vladimir  Ilyich 
Lenin  - greatest  revolutionary,  leader  of  tie  Communist  Party, 
founder  of  the  first  in  the  world  socialist  state. 

" Hith  name  and  Lenin's  activity  is  connected  whole  revolutionary 
epoch  in  huaanity's  life.  Lenin  gave  answer/responses  to  the  most 
urgent  guestions,  placed  by  the  course  of  historical  development,  he 
thoroughly  developed  the  theory  of  the  Socialist  Revolution  and 
building  Communist  society,  it  arsed  russiac,  all  the  intarnational 
revolutionary  movement  by  scientifically  substantiated  strategy  and 
tactics,  it  headed  the  fight  of  class  for  the  conversion  of  the 
ideals  of  socialism  in  life.  Socialism,  converted  by  Harx  and  Engels 
fren  utopia  into  science  and  enriched  by  Lenin  by  new 
conclusion/derivations  and  discover y/openinge,  was  personified  into 
the  social  practice  of  world-wide  historical  scales,  it  became  basic 
revolutionary  force  of  our  time".  (Theses  of  the  CC  of  KPSS  [‘ 

CPSU]  to  the  100th  anniversary  from  the  birthday  of  V.  I.  Lenin) . 

Lenin,  the  greatest  scientist  in  revolution  and  revolutionary  in 
science,  gave  enormous  value  to  guestiens  of  the  scientific-technical 
progress  of  our  country.  From  the  first  days  of  the  existence  of 
Soviet  state,  all  possible  and  comprehensive  development  of  science 
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and  technology  became  one  of  the  nest  important  and  systematic  I 

i 

directions  of  the  activity  of  the  Communist  Party  and  Soviet  ] 

governnent. 

In  heavy  1918,  when  young  Soviet  republic  in  lethal  fight 
reflected  the  brightness  onset  of  counter  revolution,  Lenin  writes 
the  "sketch  of  the  plan  of  scientific-technical  work"  for  the  Academy  i 

of  Sciences,  in  which  he  scientifically  assigned  the  mission  of 

developing  the  plan/laycut  cf  the  reorganization  of  industry  and  ' 

eccnoiic  lift  of  Russia.  This  plan/layout  strikes  with  its  depth,  1 

newness  of  posing  of  problems,  with  organic  coimunicat ion/connection 

with  life.  Onder  Vladimir  Ilyich's  manageaent/manual  in  the  same  J 

1 

period,  was  developed  the  plan/laycut  of  the  electrification  of 
Russia  - plan/layout  GCELEC  [ State  Commission  for  the 

Electrification  of  Russia],  to  realization  of  which  Lenin  gave  < 

enormous  value. 

Vladimir  Ilyich  paid  great  attention  to  development  of  Soviet 
aviation  and  technology.  Lenin  supported  great  Russian  scholarly 
professor  N.  E.  Joukowski's  proposition  about  the  organization  of 
central  aerohydrody namic  institute.  TsAGI  ['  Central  Institute 

of  Aerohydrody nanics  in.  H.  Te  Zhukovskiy  is  the  authentic  creation 
of  Great  October.  Because  of  the  daily  concerns  of  the  Communist 
Party  and  Soviet  governnent  of  TsAGI,  it  became  the  world  famous 
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scientific  research  aviation  center,  equipped  with  modern  research 
equipment,  disposing  of  the  highly  skilled  scientific 
personnel/f rames. 

Following  Lenin's  legacy,  Soviet  people  under  the 
aanagement/aanual  of  the  Communist  Party  carried  out  an 
industrialization  of  the  country,  they  converted  our  native  land  into 
mighty  socialist  power,  the  reliable  stronghold  of  peace,  progress 
and  socialism. 

Page  C. 

Science  in  our  country  ever  more  and  mere  is  converted  into  the 
direct  productive  force  cf  society.  The  Communist  Party  takes  all 
measures  for  realizing  the  Leninist  precept  about  that,  "so  that  the 
science  for  us  would  not  remain  a dead  letter  or  fashionable  phrase 
...  so  that  the  science  real/actually  would  enter  in  the  flesh  and 
the  blood,  it  was  converted  into  the  component  element  of  mode  of 
life  completely  and  by  present  focmN.  (Coll,  works,  Vol.  45,  page 
39 1)  . 


Published  during  October  1968.  The  resclution  by  the  CC  of  the 
CPSO  and  Council  of  Ministers  of  the  OSSR  "about  measures  for  the 
increase  of  the  effectiveness  of  the  work  of  scientific  organizations 
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and  the  acceleration  of  use  in  the  national  economy  of  the 
achievements  of  science  and  technology"  is  directed  toward  further 
increase  in  the  effectiveness  of  scientific  investigations. 

Soviet  scientists,  accurate  to  Lenin's  legacy,  direct  their 
efforts  for  the  solution  of  stated  before  then  by  party/batch  and 
government  most  important  problems  in  further  increase  in  the 
effectiveness  of  scientific  investigations  for  purpose  of  the 
prevision  for  technical  progress. 


DOC  = 78068001 


PAGE 

(r 

Page  1. 

METHOD  OF  CALCULATION  OF  FLOW  ABOUND  A BODY  OF  BEVOLUTION  OF  ANY  FORM 
DOBING  ARBITRARY  MOTION  IN  IDEAL  FLUID. 

L.  A.  Maslov. 

Is  proposed  the  method  of  calculation  c£  distribution  of  the 
speed,  pressure  and  potential  on  surface,  and  also  in  any  point  of 
space  around  the  body  of  resolution,  which  accomplishes  arbitrary 
motion  in  ideal  fluid.  In  comparison  with  known  methods  in  this  case 
tc  the  form  of  body  of  revolution,  are  superimposed  no  limitations 
and  sufficient  accuracy/precision  of  calculations  is  reached  at  the 
considerably  smaller  expenditures  cf  tine  EVM  [‘  computer]. 

Examples  of  calculations  are  compared  with  known  exact  solutions  and 
with  the  experimental  values  of  pressures  on  the  surface  of  different 
bodies. 

At  present  the  potential  flow  around  the  bodies  of  revolution  of 
any  fora  during  arbitrary  notion  can  be  designed  only  with  the  aid  of 
method  [1].  However,  method  [ 1 ] it  is  very  laborious,  and  more 
effective  proves  to  be  the  method  cf  calculation  of  flow  about  body 
of  revolution  during  arbitrary  motion,  presented  in  works  [2]  and 
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[3].  Method  [2],  furthermore,  easily  it  is  spread  to  the  case  of  the 
calculation  of  three-dimensional  body  during  arbitrary  motion  [4]. 
However,  during  its  use  on  surface  are  superimposed  some  limitations. 

In  this  article  is  proposed  the  calculation  method,  which 
represents  by  itself  the  generalization  of  *ethcd  [2].  In  this  case, 
the  method  of  the  assignment  to  body  surface  makes  it  possible  to 
examine  the  bodies  of  revolution  of  any  for*  during  arbitrary  motion. 

1.  Fundamental  principles.  If  vt , v2 , v3  - projection  of  vector 
v of  the  forward  velocity  of  pole  A of  solid  body,  and  Qlf  Q2,  03  - 
to  the  projection  of  angular  velocity  vector  fl’,  on  the  connected  with 
the  body  axes  of  coordinates  xyz,  then  for  the  potential  of  the 
disturbed  velocities  of  liquid  it  is  possible  tc  write 

6 

<I»  (x,  y,  z,  t)  ~ V v,  (o  (X,  y,  z), 

i _1 

where  are  introduced  designations  v4  * lut  ; v,  = i02 ; v6  = lo3  (1  - 
length  of  body) . 


Page  2. 


For  determining  each  of  six  single  potentials  <M*.  -V.  2 ) it  is 
required  to  solve  the  exterior  problem  of  Neumann 


d<t>L  v,n 

On  . v, 


= 0; 


(<  = 1.2 6) 


(1.1) 
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with  zero  conditions  at  infinity.  Here  v , - the  velocity  of  the 
points  of  body  surface  in  the  appropriate  siaple  notion;  j - surface, 
which  liaits  solid  body;  n - normal  to  the  surface  w,  directed  inside 
liguid.  If  we  the  solution  search  for  in  the  fora  of  the  potential  of 
the  siaple  layer 

f C dot 

!\(Q) 7?  ■ (1.2) 

the  boundary  condition  (1. 1)  is  reduced  to  the  integral  second-order 
of  Fredhola  equation  relative  to  the  intensity  of  layer  p 

^ (P)  + j | (*,  ( Q)  d«=  vt  (P) •n(P),  (|  .:j) 

where  P - an  arbitrary  calculation  point;  Q - current  point  of 

7 7 

surface  w;  H = QP. 

Equation  (1.3)  has  unique  solution,  if  surface  w belongs  to  the 
class  of  Lyapunov's  surfaces  [5]. 

The  longitudinal  X-axis  of  the  Cartesian  system  of  coordinates 

— ’>”7 

xyz  with  unit  vectors  ijk  coincides  with  the  axis  of  the  svaaetry  of 
body  of  revolution.  The  origin  of  coordinates  is  placed  in  the 
leading  edge/nose  of  body  (Fig.  1).  Together  with  Cartesian  system  is 
examined  the  systea  of  cylindrical  coordinates  xr0  (u  = r cos  Q,  z = 
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r sin  0)  . To  the  calculation  point  P appropriate  thenselves 
coordinates  xr0,  to  the  current  point  Q - coordinate  €pd;  R = (x-C) 

♦ (r  cos  0 - p sin  •)  j ♦ (r  sin  9 - p sin  d)  k. 

If  pole  A is  selected  on  X-axis  at  a distance  xA  froa  carrying 

body,  then  for  the  velocity  of  the  points  of  body  surface  during  its 

arbitrary  aotion  it  is  possible  to  srite 

— ¥ -*  — * 

v r sin  0 - vnr  cos  0 vt)  i -}-  \vt(x  — x.\)  — v4  r sin  6 — v,\J  -f 

+ \v4  r cos  0 — v.,  (x  xA)-v3\  k,  (1.4) 

where  the  linear  dimensions  are  referred  to  the  length  of  body 
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Equation  of  generatrix  of  body  of  revolution 


(1.5) 

is  represented  in  the  paraaetric  fora: 


r — r(s) ; x — x(s),  (1.6) 

where  s - an  arc  length  of  generatrix,  calculated  off  the  origin  of 
coordinates.  The  value  of  parameter  s,  characterizing  the  current 
point  Q of  body  surface,  is  designated  by  letter 


For  the  representation  of  the  differential  cell/eleae nts  of 
surface,  are  introduced  the  designations 


r' 


dr  , dx 
ds  ' * dT  ‘ 


(1.7) 
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The  element  of  area  of  surface  is  equal  tc 


d(o  — rdsdd . 

At  each  point  of  body  surface*  is  introduced  the  connected  with 
this  pcint  rectangular  coordinate  systen  whose  unit  vectors  are  equal 

to  rt=  — r' i -\-x'  cos  8/  f x'  sin  0£; 

v — x'  l -f-  r'  cos  0/  — r'  sin  0 k\  (1.8) 

b — — sin  Oy -f-  cos  6 k : 

here  n7  - standard;  vector  it  is  directed  tangentially  toward 

—7 

generatrix  to  the  side  of  an  increase  in  the  arc  length  s;  vector  b 
lie/rests  at  transverse  plane  and  foras  witfc  n and  rright-hand^d 
coordinate  system. 


Unlike  method  [2]*  where  the  surface  is  assigned  in  the  form  of 
{l.5)  and  is  utilized  tie  derivative  dr/dx,  in  this  case  form  ^1.6) 
makes  it  possible  to  present  any  the  curve,  which  limits  simply 
connected  region. 


As  a result  of  substitutions,  fundamental  integral  equation 
(1.3)  takes  the  form 


0)  = vtn(s,  0) 


f f u (,  ~ x'  P cos  (0  — 0)]  p</j  (/H 

JJM’  ’ R*  t-p<  2rpcos(0  »>|*  ’ 
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where  L - complete  arc  length  of  generatrix  from  forepart/nose  point 
x = 0,  r = 0)  to  the  tailed  (x  = j(,  r = 0)  . 

Analogous  expressions  are  obtained  for  velocities  and 
potentials. 

2.  Calculation  foraulas.  as  it  follows  froa  the  considerations 
cf  symmetry,  for  the  characteristic  cf  the  arbitrary  notion  of  body 
of  revolution  in  ideal  fluid,  it  suffices  tc  know  the  parameters  of 
flew  only  for  three  simple  motions: 

- f or  ward/progressive  along  x-axis  with  a velocity  of  v,  * 1 (i 

* 1)  . 

- for ward/progressive  along  y axis  with  a velocity  of  v,  * 1 (i 

* 2). 


- rotary  around  transverse  axis,  for  example  the  parallel  z axis 
and  passing  through  pole  A,  with  angular  is  velocity  v*  * 1 (i  * 6) . 


Page  A. 

The  total  values  of  relative  velocity  and  (Otential  on  the  body 

me r face  of  rotation  it  is  possible  to  present  in  the  form 

U 1 1 u , , v,  -f  Uj,(Vj  co*  h f v,slnt)  w»  ( i>  sin  »||  * 

+ h !«,*{»,  sin  A — v,cost)  f w ,(e,  tint  + rt',|.  (2.1) 

♦ — /|f,v,  + f,  (e,cosl  -i  v,  sin  6)  t f,(t»,cos*  v,  kin  t)|,  (2  2) 
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where  »i'  and  «(»  - dimensionless  components  of  relative  velocities 

cn  exes  ▼ end  b.  - the  dimensionless  potentials  in  the 

appropriate  staple  notion,  which  are  subject  tc  further 
deter ainat ion. 

Solution  is  of  integral  equation  (1.9)  for  the  siaple  notions 

indicated  should  be  searched  for  in  the  fora 

!*i  (s,  «)  — l*i  (*);  (s,  6)  = p*  (s)  cos 6;  (i, (s,  0)  = (i, (s)  cos  0.  (2.3) 

Instead  of  source  strength  ^ (S)  it  is  convenient  to  examine 
other  unknown  functions  gt(s),  connected  with  n,  the 
relationship/ratios 

t,  = 2*rv,'  (24> 

and  all  the  calculations  of  news  in  a dimensionless  fora,  accepting 
as  characteristic  linear  dimension  the  length  of  body  l. 

After  the  substitution  of  values  (2.3)  and  (2.4)  into  equation 
f(l. 9)  and  the  prolonged  conversions,  analogccs  given  in  work  [2], 
integral  equations  for  each  siaple  notion  (i  * 1,  2,  6)  they  take  the 
fora 

L 

Zi  (s)  - fi  0 (s)  - J gt  («)  K,  0 (S,  0)  do. 


(2.5) 
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Here  for  known  functions  fl0(s),  equal  to  the  product  of  the 
nornal  coaponent  of  the  velocity  of  follovicg  to  a radius  of  body  at 
the  particular  point,  are  utilized  the  relationship/ratios 

/l0  = rr'\  f„=-rx'-  /«.  = rjrr' +*'(* -*„)].  (2.6) 


The  nuclei  of  integrals  are  designated: 


+ K,0{s,  °)  = K"(s,o)  = (BHl  + x’Hi)A-, 


L . u _2r(x,(r-p)-r'(jc-  01 

« V(x  - W + (r  + fP  ’ (x  - IjF  + <r  - p)®  • 

G,  = E(0);  C,  = K(k?)—E(k*); 

™ Rl  + ft'1)  t’(ft’)  - 2ft'«  K (fe-)l; 

kL 


(2.7) 


1 [( ) + 3k'*)  a:  (k*)  - (3  + k'»)  /f  (**)|, 


where  K(k*)  and  E(k*)  - conplete  elliptical  integrals  of  the  first 
and  second  kind  with  the  aodule/nodulus 


Page  5. 


£<  _ 4rp 

<*  — 0*  + (/■- f-p)’ 


and 


* 


*'•=!-  k*. 


Analogously  are  obtained  calculated  re laticnshi p/rati  os  for  the 
conponents  of  the  dinensicnless  velocities  end  potentials,  indicated 
in  formulas  (2.1)  and  (2.2): 
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Mi,(s)=/n(s)  I J gi(i)Kn(s,  i)ih\ 

0 1 ^ i w 

«/»(s)=^/l2(s)  I -7-  I gi(°)Kl2(s,  3)rfa; 

' J 
0 

1. 

?1  (s)  •■=  — J g,  (1)  K , a (s,  5)  do;  (2.9) 

0 

fit  — r'\  fc,\—rx' — r'  (a; — .x,,i); 

~ 1 ; /6j  = jc  ; 

, (s,  0)  - fcG,  + G.)  4 ; (s,  ~ G, 

/ , \ (21°) 

(»,  J)  - (s.  0)  - ( Ctf , + - H,  j .4 ; Kn  (s,  0)  = AC62  (s,  o )=  Ht  A. 

where  besides  the  values,  determined  to  foraulas  (2.7) , are 
introduced  the  designations 

r — (x~  ^ 1 r If  ~JL  • q „ 2K  ( k *)• 

(2.11) 

//,=  £-|(l  + k'')K(V)-2E(k% 

It  the  end  points  where  x = 0,  r = 0,  or  x = 1,  r = 0,  are  not 
difficult  to  show  that  integral  of  equation  (2.5)  always  have  zero 
solution  and  do  not  require  special  examination,  as  is  done  in  method 
[2].  Durinq  the  deterainaticn  of  the  values  of  the  velocities  and 
potentials  in  end  points,  it  suffices  to  calculate  the  coaponents  of 
diaensionless  velocity  duricg  transverse  and  rotary  notions  (i  = 2.6) 
along  y axis 


1. 

U -C  f fiW** 

>y  1 J 2 |(x  — E)1  f p»p 


(2.12) 
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and  diaensionless  potential  during  the  longitudinal  flow 


J (Ur-  V+rt1*  ’ 


(2.1.1) 


■her*  C,  • 1,  for  forepart/rose  turn/sharpen  one  should  place  x = 0 
a»drfi  — xa.  but  for  the  tailed  point  x * 1 ;\nd  c6  — x,\ — 1.  Tha  remaining 
components  of  the  velocities  and  potentials  in  these  points  are  equal 


to  zero. 
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The  velocities  in  the  points  of  space,  which  do  not  belong  to 
body  surface,  are  conveniently  calculated  in  cylindrical  coordinate 
systea  *•  For  the  components  of  the  diaensionless  velocity  in  each 

simple  motion,  it  is  not  difficult  to  obtain 


«/,(*.  0 = //: .W+  f r,  a)  do; 

0 

L 

«»,(■*,  r ) -/mW  + f 4Tj(j)K<«(*.  r.  3)^3; 

0 

l 

u, . ( x , r)  = /(,  (x)  -f-  — f & , (a)  \\  (x,  r,  o)  do. 


(2.H) 
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In  this  csss,  to  forwula  for  uiH  it  coincides  with  formula 
(2.8)  for  ulb,  and  for  th«  cue  lei  of  the  integrals  of  first  two 
foraalas  (2.1%)  are  introduced  the  designations 

-^7— Ar14  = (2^ci+i.c,j/i; 

Ku  = = ^r7-°  .4 ; - (2  I ■ /7,^ . 

Here  r*i  * (x-t)*  ♦ (r-p)*,  and  for  renaining  values  are  used 
designations  (2.7)  and  (2.10).  The  dinensiocless  components  of 
velocity  of  following  //»(-*)  and  /,«(•*)  are  egual  to  fl8  = f24  = 1; 

fl4  = ^23  ~ * 0l  fm^XA—X. 


Por  the  calculaticns  of  pressure  coefficients  in  the  case  of 
arbitrary  notion,  one  should  use  Lagrange's  integral.  During  the 
calculation  of  apparent  additional  masses,  it  is  necessary  to  bear  in 
aind,  that  for  a body  of  revolution  independent  variables  and  not 
egual  to  zero  are  identically  four  apparent  additional  masses:  Xlt, 

X 2 2*  X 2 6 and  X*4.  By  analogy  with  work  [3]  it  is  not  difficult  to 
obtain 


L 

- 2icp„  ('  rr' p,  ds\ 

II 

L 

L,2  icp0  f rx'  f3  ds; 

0 

I. 

— "f-0  Jr  \rr'  t x'  (x  xA)\f,ds; 

0 

L 

*««  = - *P«  j r | rr'  + X'  (X  — A"^)]  f,  ds, 

0 I 


(2.15;, 


■ -7;  •»  .- 


where  pQ  - mass  density  of  liquid 


1 
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3.  Methods  of  calculations.  The  scluticn  of  the  fundamental 
integral  equation  of  problem  (2.5)  is  determined  in  the  finite  number 
of  discrete  calculation  points  of  tody  surface.  The  the  intricate 
shape  of  body,  the  greater  the  calculation  points  it  is  necessary  to 
select  for  the  sufficiently  precise  description  of  generatrix.  So,  in 
that  comprised  for  computers  "Hinsk-2"  to  tie  program,  according  to 
which  were  fulfilled  examples  of  calculations,  can  be  utilized  to  160 
calculation  points,  arbitrarily  arrange/located  on  by  generatrix 
bodies  of  revolution.  For  example,  for  smooth  bodies  sufficient 
accuracy/precision  is  reached  at  50  calculation  points. 

Page  7. 

The  proposed  method  allow/assumes  the  presence  of  the  finite 
number  of  salient  points  of  the  enclosures  cf  generatrix.  At  very 
salient  point  for  formal  satisfaction  of  Lyapunov's  conditions,  one 
should  assume  a small  bending  radius.  In  the  process  of  calculations, 
the  rounding  is  re alize/accomplishe d automatically  applying  quadratic 
interpolation  between  calculation  points.  In  the  places  of  an  abrupt 
change  in  the  enclosures,  is  necessary  the  packing/seal  of 
calculation  points. 
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Basic  difficulty  is  the  calculation  of  improper  integrals  in 
equations  and  foraulas  for  velocities  and  potentials.  By  analogy  with 
work  [2]  improper  integrals  are  calculated  vith  the  aid  of  the 
replacement  of  variables: 


3 — s = sign  ( h ) 


0,5  Lh> 


Inverse  dependence  will  be  determined  ty  relationship/ratio 


dh  1 


’ 1 3 ~ * I (1 0 — * I + 0,5  £)*■ 

n L 


Reference  point  alternating/ var iafcle  h coincides  with  the 
special  feature/peculiarity  of  integral,  which  is  located  in  the 
calculation  point  whose  coordinate  s projects  in  this  formula  as 
parameter.  Functions  (3.2)  and  (3.3)  are  continuous  and  different 
frca  zero  everywhere,  with  the  exception/e  1 imination  of  most  singular 
pcint  where  (3.3)  it  vanishes  as  |s-«|V2. 

By  analogy  with  work  [4]  replacement  is  utilized  only  on  the 
section 


| a — a |< 0,054,  (3.4) 

i.e.  near  singular  point.  Integrals  are  represented  in  the  form 


L s—0  05/.  L \)V  \\ 

| z(i)K(s,o)di=  j g(a)Kdi  -r  J g(i)Kth  | j (3.5) 

n o t+n.o!w.  i / V i i 

In  formula  (3.5)  with  s ^ 0.05L  first  term,  but  with  s >,  0.95L 
icond  term  they  are  not  considered,  since  in  these  cases  they  are 
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included  in  the  third.  In  this  case,  tie  liiits  of  the  third  integral 
respectively  change  and  are  calculated  directly  from  formula  (3.2), 
where  cne  should  place  • = C in  the  first  case  and  a = L in  the 
second . 

The  first  two  integrals  of  formula  (3.5)  do  not  have  a special 
feature/peculiarity  and  are  calculated  with  the  aid  of  trapezoidal 
rule  on  the  node/units  of  integration,  arrange/located  in  calculation 
points.  Last/latter  integral  is  calculated  from  trapezoidal  rule  with 
the  constant  space  in  new  a Iter nat i ng/variable  h.  In  this  case 
replacement  (3.2)  provides  the  symmetrical  location  of  the  node/units 
of  integration  relative  to  special  feature/peculiarity  with  an 
increase  of  the  density  of  the  location  of  node/units  in  its 
vicinity,  which  corresponds  to  the  conditions  for  existence  of 
principal  value  for  Cauchy  and  it  makes  it  possible  to  calculate 
improper  integrals  according  to  trapezoidal  rule  with  the  constant 
space  in  new  variables,  necessary  in  this  case  values  g ( ®)  , r<o)  and 
p(«)  in  the  node/units  of  integration,  which  prove  to  be  between 
calculation  points,  are  determined  by  quadratic  interpolation 
according  to  Newton  [6]. 

Integral  equations  are  solved  by  the  method  of  successive 


approximations  according  to  the  following  diagrams: 
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for  i * 1 

gv>  = ~)2f\*\  K'lk  *“/io  - fffi.*  \K„,dr,  frlk  ^ (/?**  + *_i); 

for  i = 2.6 


ft  i o - f i o-  K ik  — ft  o \ gt.  *-i  Ki  o d'i, 

where  k - a number  of  aFproach/appicximaticn. 

Page  8. 


Solution  is  considered  found,  if 

I (Zik  («)  - gt.  * I (s)  n.»  < 0.M5  }/,  0 

i. e.  boundary  conditions  (1.3)  are  fulfilled  with  accuracy/precision 
by  0.5o/o.  The  number  of  approach/approxiaaticns  oscillates  from  4 
for  smooth  bodies  to  10-15  for  the  bodies  of  intricate  shape. 

Information  about  body  is  assigned  by  the  tables  of  values  x and 
r calculation  points.  Derivatives  (1.7)  are  calculated  with  the  aid 
of  Newton's  formulas  from  values  of  x and  r at  three  calculation 
points  [6].  Integrals  in  the  formulas  of  apparent  additional  masses 
(2.15)  are  calculated  from  trapezoidal  rule.  Space 
alter nating/var '.able  h was  selected  as  being  egual  to  Ah  = 

1//TFT5. 


EOC  = 78068001 


PAGE  23 


K.  Examples  of  calculations.  For  purpose  of  checking  the  method 
presented  on  computers  "Hinsk-2”  for  ellipscids  of  revolution  are 
carried  out  the  calculations  of  apparent  additional  masses  and 
relative  velocities  for  which  are  known  precise  analytical 
expressions  [7].  The  applicability  of  method  to  the  determination  of 
pressures  in  real  liquid  is  shown  tascd  on  the  example  of  the 
calculation  of  the  body,  which  has  the  vertical  section  of 
generatrix,  for  which  could  not  be  ottained  the  solution  by  method 
[2],  or  bodies  in  the  form  of  the  combination  cf  the  cone  with 
cylinder,  having  a local  abrupt  change  in  the  enclosures. 

In  Fig.  2 and  3 points  plotted/applied  the  results  of 
calculation  by  the  proposed  method  of  dimensionless  velocities  on  the 
surface  of  the  elliptical  disk,  which  has  the  relationship/ratio  of 
semi-axes  a/b  * 4E  0.1,  and  the  ellipsoid  which  has  a/b  = 9.  By  solid 
lines  are  constructed  precise  values  of  velocities.  Calculation  of 
velocity  field  is  checked  on  the  example  of  the  flow  around  the 
sphere  of  single  diameter  of  f crward/progressi ve  flow  along  X-axis. 
Table  1 gives  corrected  values  of  velocity  u, , the  points  of 
vertical  diameter  (x  = 0.5),  also,  in  the  points  of  horizontal 
diameter  (y  = 0),  arrange/located  on  different  distances  from  the 
surface  of  sphere. 
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Pig.  2. 

Fig.  2. 

Key:  (1).  Exact  solution. 

Fig.  3. 

Key:  (1).  Exact  solution. 
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Fig.  3. 


(2).  proposed  method. 


(2).  proposed  method. 


Results  of  the  calculation  of  the  coefficients  of  apparent 
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additional  masses  for  three  ellipsoids  by  tte  proposed  method  are 
given  in  comparison  vith  precise  values  in  table  2. 

The  comparison  of  precise  velocities  and  apparent  additional 
masses  vith  the  results  of  calculations  according  to  the  proposed 
method  they  testify  to  sufficient  accuracy/precision  of  the  latter. 
The  disagreement  of  precise  and  calculated  values  comprises  less  than 
Ic/o. 


Pig.  4 and  5,  show  the  comparison  of  tte  calculated  and 
experimental  values  of  the  coefficient  of  pressure  p in  the  nose 
section  of  the  body  surface,  which  contain  the  vertical  section  of 
duct/contour.  Fig.  4,  shows  the  enclosures  cf  bodies  of  revolution 
1-3  and  is  given  distribution  p at  zero  angle  cf  attack  ». 

FOOTNOTE  *.  Bodies  of  revolution  1-3  correspond  to  engine  nacelles 
Ho.  25,  85  and  87,  investigated  experimentally  in  work  [8]. 

EHDFOOTHOT  E. 

In  Fig.  5,  is  constructed  distribution  p according  to  one  meridian  of 
body  of  revolution  3 at  the  angle  cf  attack  a = -10°,  when  this 

meridian  is  windward,  and  the  angle  of  attack  a = i 10°,  when  this 
meridian  becomes  leeward.  The  agreement  of  the  results  of  calculation 
with  experiment  proves  to  be  sufficiently  good. 
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table  1. 


(') 

OTCTomme 

OT  flOBepXIIOCTH 

( 3iia«ienHH 

C KOpOCTH  W,r 

e 

OiCTOHime 
“f  IIOHCpJUIOCTH 

A’  in  nut- huh 

CKOpOCTH  U 1 1 

r>) 

no  irpe/iJia- 
raeMOMy 
Me to ^y 

WJ 

TO  MIII4e 

3Hfl4C- 

II  If  H 

no  np<vi ,i a 
racMoMy 
MOTOAy 

& 

TO'IHhll' 

| MiaMfllMH 

(?>) 

1 1 0 BOpTHKa;ib- 
HOMy  /IHflMOTpy 

0,01 

1,470 

1 ,471 

((fl) 

llo  ro|)ii.ioiira.ii>- 
uoMy  anaMPTpy 

0,01 

0,0587 

0,577 

0,1 

1 .288 

1,289 

0.1 

0,4213 

0,4213 

2,0 

1 

1,004 

1 ,004 

2,0 

‘ 1 

0,9020 

0,9920 

Key:  (1).  Distance  froa  surface.  (2).  Values  of  velocity.  (3). 
according  to  proposed  aethod.  (4)  . precise  values.  (5)  . According  to 
vertical  diaaeter.  (6)  . According  to  horizontal  diameter. 


Table  2. 


a 

*.i  | 

Kn 

Km 

b 

| (I)  no  npe. 

.iJiaraeMOMy 

MPTOfly  | 

( Xj  TOUIUC  3 it  a 'ip  h m a 

0.1 

6,130 

0,0751 

i _ 

6,184 

0,0748 

1,0 

0,499 

0,500 

— 

0,500 

0,500 

— 

9,0 

0,0244  i 

0,950 

0,860 

0,0244 

0,954 

0,864 

Key:  (1).  according  to  the  proposed  aethod.  (2).  precise  values. 
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Pig.  4. 

Key:  (1).  the  proposed  method.  (2).  body. 

Page  10. 

The  saae  good  convergence  of  the  calculated  and  experiment.il  values 
of  pressure  is  obtained  for  the  nose  section  of  the  body,  which 
represents  by  itself  to  the  combination  of  cone  with  cylinder  (Fig. 
6). 

It  east  be  noted  that  the  results,  given  to  Fig.  4,  for  a body 
with  the  vertical  sections  cf  generatrix  at  zero  angle  of  attack  can 
be  obtained  with  the  aid  of  aethods  [1]  and  [9],  while  the  case  of 
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flow  at  angle  of  attack  it  can  be  designed  only  according  to  nethod 
[ 1 ].  In  nethod  [9]  by  virtue  of  the  use  of  function  of  current,  can 
be  examined  only  axial  xoticn  of  body  of  revolution.  Examples  of  the 
calculations  of  the  bodies  cf  revolution,  which  have  the  local  abrupt 
changes  in  the  fora  of  enclosures,  given  to  Fig.  2 and  6,  they  are 
encountered  only  of  the  authors  of  work  [1].  However,  in  nethod  [1] 
the  volume  of  calculations  proves  to  be  considerably  greater  than  in 
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Pig.  S (o. 

Fig.  5. 

Key:  (1).  the  proposed  aethed.  (2).  experiment. 

Fig.  6. 
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Page  11. 

HY  DBODYNAM  ICS  OF  THIN  FLEXIELE  BODY.  (Estimation  of  hydrodynamics  of 
rippled  surfaces). 

G.  V.  Logvinovich. 

For  purpose  of  the  explanation  of  the  techanism  of  the  floating 
of  fishes  by  the  aethod  of  flat/plane  sections  is  studied  the 
hydrodynamics  of  the  fine/thin  body  being  deforced.  Are  obtained  the 
siiple  formulas,  which  make  it  possible  to  evaluate  thrust/rod  and 
the  spent  power  during  the  sinusoidal  wave  strains  of  the  axle/axis 
cf  body. 

Are  estiaated  the  hydrodynaaic  characteristics  of  fishes.  Is 
given  the  comparison  of  the  results  cf  theory  with  experimental 
materials. 

Hydrodynamics  of  the  fine/thin  tody  being  deformed,  which 
accomplishes  during  forward/progressive  uniform  notion  small 
undulations,  can  be  sufficiently  simply  studied  via  the  application 
cf  the  aethod  of  flat/plane  sections  with  tbe  use  of  a concept  of  the 
"pierced  layer"  [1].  Twc-di aensiona 1 picbleis  cf  such  kind  were 
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theoretically  examined  by  Sikaan  [2],  that  also  placed  some 
experiments,  and  by  V.  A.  Eroshin  [3],  that  developed  the  theory  of 
the  airfoil/profile  being  deforaed  in  L.  I.  Sedov's  setting  [4], 
However,  the  study  of  spatial  problem,  apparently,  in  larger  measure 
approaches  us  understanding  of  the  aechanisx  of  floating  of  fishes 
and  narine  aniaals,  than  the  solution  of  t wc-di aensional  problems. 

The  use  in  this  case  of  a method  of  flat/plane  sections  is  justified 
by  the  fact  that  the  bodies  of  many  marine  animals  are  very  elongated 
lengt  hwise. 

1.  Let  us  examine  motion  of  slender  body  in  inertial  system  of 
coordinates  xyz,  which  moves  in  unlimited  volume  of  liquid  along 
axle/axis  Ox  with  a constant  velocity  of  v (Fig.  1).  The  longitudinal 
curvilinear  axle/axis  of  the  body  s being  deformed  weakly  differs 
from  axle/axis  Ox,  and  its  strain  in  plane  xOy  let  us  designate 
n(x,  0.  the  abscissas  of  the  ends  of  the  body  let  us  designate  xt  and 
x2;  the  length  of  body  x2  <v  I-n  let  us  assume  that  the  cross 
sections  of  body  are  forsed  by  ellipses  with  the  semimajor  axis  R = 

R (x)  , parallel  axis  Oz. 
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Fig.  1. 

Page  12. 

In  connection  with  assumption  about  the  fineness  of  body,  we  assume 
that  at  the  entice  length  dR/dx  thece  is  a value  small.  Let  us  assume 
that,  passing  through  certain  "penetrable  layer"  motionless  relative 
to  the  quiescent  liquid,  body  in  this  layer  gives  rise  to  the 
transverse  almost  plane  flow,  close  to  the  flow  of  ideal  fluid.  On 
trailing  edge  (on  the  tail  of  fish)  with  x = x,  is  fulfilled 
Joukowski's  condition  about  the  finiteness  of  velocity,  and  after 
body  remains  the  trace,  equivalent  to  the  film  cf  eddy/vortices  with 
elliptical  circulation  distribution  according  tc  spread/scope  z = 

♦ RiUi).  This  model  of  flow  in  essence  represents  the  development  of 
Jones*  known  diagram  in  connection  with  the  low-aspect-r at io  wing 
being  deformed. 
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2.  On  element  of  length  of  body  ds,  acts  ncraal  force  dFn  and 
suction  force  d P,  which  is  caused  ty  so-called  circular  pressure  dQ. 
As  is  Known,  the  apparent  additional  mass,  which  is  necessary  for 
unit  of  the  length  of  ellipse,  is  «*(*)  = p»F2(x),  normal  to 

axle/axis  s the  velocity  of  layer  v ch‘ 

(it  v Ox' 

dFn=  ^ (m*vn)ds.  (2.1) 

Circular  pressure  we  find  as  a result  cf  contour  integration  of 
the  cross  section  of  body  £♦  cf  overpressure  p — p„  */’  /(s*j, 
determined  only  by  velocity  head.  Specific  suction  force  is  equal  to 


OP 

ds 


— —{’(/>  — />„)cos(/i,  x)ds*. 


(2.2) 


Thus,  for  instance,  for  a cylinder  with  circular  cross  section 

overpressure  on  its  surface  at  points  B,  9 is  determined  by 

expression  p — p^—[‘V"  (i  . 4stn-0)  ) P Cos '>  ji_  ^ ^ Normal  to  longitudinal 

» /t!  dt 

axis  force  (2.1)  give?  the  integral 


2* 

— | (p  — Po)cosl)d(M)—  - ((<*  R*  vn), 

0 


circular  pressure  - 


dQ 

ds 


u 

- ({P-P*)d  (/?«)-- 

J 2 

o 


This  circular  pressure  is  negative,  i.e.. 


it  attempts  to  expand 


DOC 


7806800  1 


PAGE  36 


cross  section.  If  body  is  pointed  from  the  front,  then 
slope/ incl ination  toward  axle/axis  Ox  of  generatrix  dR/ds  - dR/dx  is 
alsc  negative  and  on  each  urit  of  length  acts  the  directed  forward 
(along  axle/axis  Os)  force 


Page  13. 


dP  Wl  dR  v-n  dm*  (s) 

ds  2 ds  2 ST 


(2.3) 


It  is  usually  considered  that  suction  force  appears  as  a result 
of  the  action  of  infinite  negative  pressures  on  infinitesimal  leading 
wing  edge;  here  it  is  formed  on  entire  length  of  body  as  a result  of 
acting  the  stagnation  pressures.  Let  us  note  that  in  service  record 
for  overpressure  p-p0,  used  to  the  expanded  opening/aperture  in 
layer,  is  reject/thrown  the  term  (formally  infinite),  caused  by  the 
symmetrical  expansion  of  cylinder.  It  is  possible  to  show  that  this 
is  admissible  during  the  use  of  a hypothesis  of  flat/piane  sections 
for  calculating  the  forces,  which  act  on  slender  bodies  [ 1 ]. 

Por  determining  the  force  in  the  case  of  elliptical  cross 
section,  we  will  use  seccnd-order  of  Lagrange  equation  in  connection 
with  kinetic  energy  in  layer  r m*  (R)  V"  atd  as  generalized 
coordinates  and  velocities,  let  us  accept  for  determination  of 
suction  force  the  semimajor  axis  of  cross  section  H,  while  for  the 
determination  of  lateral  force,  - velocity  vH.  As  a result  for 
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elliptical  cross  section,  will  be  obtained  the  same  expressions  of 

force  gradients,  as  for  circular:  __  ()T  _ _dm*  vi  _ t'vi,  _ dQ 

OR  OR  '1  ~ K 2 ds  ■ 


d OT  d . „ • d 

dt  dv^  dt  im  V"]  dt 


drm 

ds 


Apparent  additional  lass  ■*  in  cur  case  is  a function  only  x (or 
s)  . Poraula  (2.3)  is  valid  for  any  elliptical  cross  section,  and  in 
particular  for  the  case  when  minor  axis  vanishes  and  ellipse 
degenerates  in  the  segment  cf  line  4B  (cn  tail)  . 


Projections  on  the  axis  of  the  coordinates  of  the  elementary 


forces,  applied  to  body,  mill  be 


d/\ 

dl 


- dt, 
4FA 


<hi 

Ox 
dP 

Ox 


dP- 


(24) 


$'\a 

Integrals/  , f d/'x f\  — fdFfi*  undertaken  alcng  the  length  bodies, 

will  give  resulting  forces.  During  integration  is  to  consider 
conditions  R(x2)  = 0 and  E»(xt)  = H,.  Condition  P(xt)  = Rt  realizes 
the  disruption/separation  of  the  vortex  sheet  from  the  tail  of  body 
(Zhukovskiy  condition).  Since  the  forces  are  determined  for  the 
elementary  pierced  layer  which  rests  relative  to  motionless  liquid, 

let  us  note  that  0 i • (i 

dt  dt  Ox 


3.  During  periodic  motion  average  during  period  resulting  pull 
force,  caused  by  "flowing"  frca  tail  acoentum/iapulse/pulses,  it  is 
possible  to  calculate  without  integration  (2.4).  Actually,  liguid  on 
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unit  of  the  length  of  trace  contains  acaentua/iapulse/pulse  (jr, ) c>„ , . 
directed  along  the  noraal  velocity  of  tail.  Projection  on  axle/axis 
Ox  of  the  noaentua/iapulse/pulse,  "which  flews"  from  tail  the  second, 
undertaken  with  opposite  sign,  gives  the  instantaneous  force,  which 
acts  cn  body.  Therefore  pull  force,  obliged  ty  its  origin  to  the 
mo aentua/i apulse/p ulses,  left  in  trace,  will  be 

/=  + «*( Xl)v(dl  — V^\dlt  (3.1) 

\ (it  Ox  I Ox 

Key:  (1).  with. 

Page  14. 


For  the  explanation  of  the  aforesaid,  let  us  give  following 
reasonings.  Let  the  body  be  is  enveloped  by  sufficiently  distant 
ccitiol  surface  Ii,  which  moves  together  with  fcedy.  Another  control 
surface  I*  of  the  saee  configuraticn,  as  vti  is  connected  with 
stationary  liquid.  At  certain  aoment  of  tiae  t,  both  control  surfaces 
coincide.  During  period  r,  control  surface  will  aove  along 

axle/axis  Ox  to  cut  Vr. 

Since  the  flow  of  liquid  inside  at  tor que/moaents  t.  and  t + 
v is  identical,  all  the  increase  in  the  kinetic  energy  and 
aoaentua/inpulse/pulse  cf  liquid  will  be  caused  only  by  track 
segaent,  which  reaained  between  the  rear  walls  of  control  surfaces 
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ill  and  i'2  According  to  the  prerequisite/ptenises  accepted,  after 
"runoff"  fron  the  tail  of  fish  each  eleient  cf  length  of  trace,  being 
deforned,  retains  nonentun/inpulse/pulse  arc  energy;  therefore  for 
calculating  the  composite  force,  it  is  possible  for  these  values  to 
take  the  values  of  specific  inpulse  tbich  relate  to  the 

torgue/aonent  of  descent  fron  tail. 


Thus,  the  average  value  of  pull  force,  caused  by  the  left  in 

trace  nonentun/inpulse/pulses,  it  nill  be  {/I  — ^ | Pit.  Value  I under 

0 

integral  is  calculated  for  point  xt,  i.e. , icr  a tail. 


Another  portion  of  pull  force  is  realized  in  the  forn  of  suction 

force,  which  appears  as  a result  of  the  lateral  flow  around  body. 


P - 


r dm*  (x)  y 

J dx  " 2 


(dr, 
\ dt 


(3.2) 


Let  us  note  that  the  appearance  of  suction  force  P leads  to  the 
fact  that  aonentua  vector  in  the  trace  is  turned  on  certain  angle  so 

that  average  total  pull  force  is  equal  to  { I } ♦ { P },  where 

f P)  - -[-J Pdt. 

Average  active  power  of  pull  forces  is  egual  to  ( { I } ♦ { P }) 

* - ( * )• 


t ron  tail  into  trace  on  unit  of  path  "flows"  the  kinetic  energy 


\_ 

9 


I (h, 

V* 


v 


dr,  y 
,ix  )m.m. 


~ h. 


(3.3) 
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The  power,  spent  on  the  excitation  of  pull  force,  will  be  { K j 
= ({1}  ♦ CP}  ♦ C E } ) V,  and  efficiency  |ri/()  = J^p 

The  given  above  fomulas,  used,  for  exanple,  to  delta 
low-aspect-ratio  wing  with  spread/scope  2R,,  give  the  Known  results 
cf  linear  theory.  It  is  real/actual,  at  constant  angle  of  attack  a, 
velocity  vn  = — i/a  and  fron  (2.1)  is  obtained  ~ m‘V'2a  = pw  R21Vs?a, 

but  frcn  (2.3)  after  integration  fron  xt  to  x2  we  have  P = 1/2 
n*tV2«*. 

Page  15. 


The  integrals  of  expressions  (2.4)  frcn  x»  to  x2  give  induced  drag 


and  lift 


Fx  — m]  V' **  + -J  m]  V*  r*  } p rtf  V 5 a* 

Fy  = F,  + P« = p«Ri  + -!«•+  . . 


Analogously  under  the  sane  assumptions  can  be  examined  the 
general  case  of  the  unsteady  fluctuation  of  low-aspect-ratio  wing. 


4.  Undulations  of  body  will  be  obtained,  if  law  of  strain  is 
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_ / Ct  X-.  — X \ 

assigned  in  the  fora  Ti  — Tl#slnl  ^ f !•  assuming  that  length  of 

body  x,  — x,  =-  /-„  - ‘2ntn  (n  nunber  of  waves,  which  are  placed  at  length 

of  body)  ; C - phase  wave  velocity,  progressive  back/ago. 

The  anplitude  of  the  strain  of  the  axle/axis  of  body  in  the 
general  case  can  be  represented  together 

^Jo  = «0  t-  (x,  — x)  + a2  (x,  - x)3  + ...  . 

Let  us  examine  below  the  simplest  case,  after  accepting 
r(u  = const,  lorsel  velocity  will  be 

vn  = V(C  - l/)  cos  I C/  - 

the  period  of  oscillations  r = 2a^/C^ 

Generally  speaking,  the  law  of  the  distribution  of  apparent 
additional  nasses  a*  (x)  is  essential  during  the  calculation  of  force 
of  P.  The  average  value  of  suction  force  is  calculated  fron  the 

fornnla  r,  A 

(P!  =-  - j-J"  dt  j t]  ~&rdx  = ~ i ( Wdm*- 

0 xt  • 

m\ 

In  oar  case  average  value  K>  - ~J  )* X (C-T)  * is 

0 

constant  in  range  froa  x,  to  xc;  therefore  ( P } * { E } * 

— lV2\m'  the  law  of  the  distribution  of  apparent  additional  nasses 

2 ' 1 

aloag  the  length  l.r  proves  to  be  unessential. 
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Page  16. 

Substituting  these  expressions  in  the  preceding/previous 
formulas,  let  us  find  the  average  values 

|pi-T”';(f)v(-F-'y-<£]’ 

~ 4 ("f ) ''  [ ( T')  “ 1 J ■ ! (4.1) 

Fig.  2,  gives  relative  average  values  ( I } and  { p } of  those 
comprise  of  impulsive  and  suction  force,  obtained  from  (4.1)  by 
division  on  ~2~W*(l#)  In  the  conditio ns/acde  of  floating  with 

high  efficiency,  when  >0,y,  the  portion  of  suction  force  does 

not  exceed  20o/o  of  inpulsive  force.  Fig.  3f  gives  the  kinogram  of 
the  floating  of  nackerel,  borrowed  frca  S.  G.  Aleyeva's  work  [5].  On 
this,  photographs  it  is  possible  to  conclude  that  C/V  ft  2 and  the 
efficiency  is  close  to  0.75;  wavelength  soaewhat  lesser  than  the 
length  of  body  and  *0,4. 


DOC  * 78068001 


PAGE  43 


Page  17. 


If  we  consider  that  the  active  thrust  1/V  { A } overcoats  resistance 

r>V? 

of  friction  u /=C/S then  the  coefficient  cf  friction  drag  will 
be  located  fros  equality  H * 1/V  ( A } . Por  C/V  = 2 and  y —0,4  is 
obtained  the  following  estiaation  of  drag  coefficient:  c,  ^ 0,24 
where  S - the  aoistened  body  surface. 


According  to  experiaents  of  V.  Ye.  Pyatetskiy  [6]  for  bluefish 
42  ca  long  with  the  spread/scope  of  tail  of  approximately  8 ca  at 
velocity  V = 0.55  a/s  were  obtained  values  C/V  = 1.47  and  ^ = 


DOC  = 78068001 


PAGE  44 


0.19-0.25-  Average  efficiency  0,84  In  tfcis  case  »R2,  = 50  cm2  and 
S ~ 640  cn2;  therefore  the  coefficient  of  friction  drag  was  obtained 
by  close  to  1.7-3.  0*  10~3. 

Calculations  according  to  formulas  (4.1)  show  that  at  values  C/V 

* 1.5,  - — 0,2  the  "wave  iotor"  with  spreao/sccpe  1 n at  velocity  v 

* of  10  a/s  develops  thrust/rod  of  apprcxinatel  y 90  leg,  spending 
power  of  a pproxina tel y 11.8  kW. 

The  aethods  of  the  evaluations  of  the  effectiveness  of  "wave 
actors  presented"  can  be  developed  for  the  configurations  of  body  and 
laws  of  strains,  closer  to  those  observed  at  high-speed  fishes.  It  is 
possible  by  the  sane  way  to  consider  tie  effect  of  back  and  ventral 
fins,  variable  along  the  length  of  aaplitude  and  series  of  other 
observed  in  nature  factors.  Integral  estinatiens  according  to  the 
trace,  left  by  body,  are  interesting  by  their  generality.  It  is 
inportant  to  note  that  the  used  nethod  was  checked  in  the  cases  of 
gliding  and  notion  within  the  liquid  cf  solid  bodies  where  the 
results  of  theory  and  experinent  proved  to  fe  very  close.  It  is 
possible  to  expect  that  and  this  estination  cf  the  propulsive 
properties  of  fishes  is  close  to  reality.  This  to  a certain  extent  is 
ccnfirned  by  the  satisfactory  convergence  of  the  calculated  values  of 
poll  force  and  frictional  resistance  for  the  fishes,  inspected  by  V. 
Te.  Pyatetsky. 
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THEORY  OF  UNSTEADY  CURVILINEAR  HOTICN  OF  LIFTING  SURFACE  IN  GAS. 

V.  E.  Baskin. 

Is  examined  the  general  case  of  the  unsteady  curvilinear  motion 
of  lifting  surface  in  gas  (cn  the  kasis  of  linear  theory). 

For  this  surface  are  given  the  formulas,  which  express  gas 
velocity  through  the  dersity  of  distributicr  of  eddy/vortices.  These 
formulas  generalize  Biot-Sa vart *s  usual  law  in  such  a way  that  it 
becomes  suitable  for  arbitrary  transient  vertices  in  gas  (in  linear 
approach/approximation).  Generalization  lies  in  the  fact  thae 
Eiot-Savar t's  usual  formula  is  applied  to  vertex  elements  takinq  into 
account  delay  in  the  fornation/educaticn  of  velocity  to  the  transit 
time  of  sound  signal,  and  appear  seme  suppleaectary  "wave"  component 
of  induced  velocities. 


The  numerous  methods  of  soluticn  of  direct  and  reverse/inverse 
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problems  of  the  flow  around  lifting  surface  of  gas  [1]  - [3]  they  are 
related  to  rectilinear  motion.  The  stationary  helical  motion  of 
carrying  filament  in  gas  is  investigated  by  Frarkl  [4],  and  lifting 
surface  - by  Naykapar. 


Gas  velocities  during  the  curvilirear  notion  of  carrying 
filament  are  examined  in  [5]. 

1.  Let  infinite  gas  to  the  torque/moment  of  time  t = 0 be 
rested,  and  then  it  was  agitated  by  the  motion  in  it  of  lifting 
surface.  Of  perturbation  rates  we  set/assume  much  the  lower  speed  of 
the  motion  of  the  points  of  lifting  surface.  During  the  calculation 
of  velocity  fields,  we  consider  permissible  the  transfer  of  the 
points  of  application  of  force  to  gas  from  lifting  surface  to  certain 
the  closely  spaced  to  it,  permeable  for  gas  surface.  The  vector  of 
the  surface  density  of  the  applied  forces  we  set/assume  normal  to  the 
permeable  surface  indicated,  continuous  at  each  moment  of  tine  within 
certain  chat  driving/mcving  along  the  surface  of  region  and  the  equal 
tc  zero  out  of  this  region.  Under  these  conditions  is  placed  the 
problem  - to  determine  the  rates  of  flew  of  gas,  if  the  motion  of 
lifting  surface  and  force  on  it  are  known. 


Page  19. 
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2.  Disturbed  velocities  v and  pressure  p cf  infinite  gas,  set  to 
motion  by  arbitrary  field  of  external  volume  forces  in  linear 

approach/approximation  it  is  deteriined  by  system  of  equations 


grad  p yj:\  <f-\c'in<livv  0 


and  by  initial  conditions 


v , 0-O,  p , o 


(p0  - density  of  undisturbed  air,  c - speed  of  sound). 


If  we  present  the  velocity  and  pressure  in  the  form 


v =/3  + grad  <p,  P =-  — Po  ()f 


the  systea  (2.1)  it  is  reduced  to  ncnhc  irogeneous  wave  equation  for 


the  potential  ?: 


1 f-l  div  D, 


where 


U « j Fdt 


the  referred  to  p0  vector  cf  the  momentum  density  of 


external  forces. 


Potential  * is  conveniently  searched  for  in  the  fora  tp  = div£, 

— 

which  reduces  to  following  equation  for  vector  E: 


l (P  £ A 

«•*-  ? ....  " 
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Solution  of  aquation  (2.3)  under  the  zero  initial  conditions, 
which  correspond  (2.2),  is  given  by  Kirchhcff's  formula 

E(ro<  tn)  — jjj  -—■b(r,  t„  — 'c  jdfe,  /=r;rg_r  (2.4) 

(integral  it  is  common  on  the  part  of  the  space  where  integrand  is 
not  equal  to  zero)  . 

Vector  E,  accordingly  (2.4),  is  constructed  as  delaying 
potential  of  flows  the  vector  equal  in  momentum  density  to.  Therefore 
let  us  name  this  vector  momentum/impulse/pulse-potential . Passage  to 
the  limit  to  the  surface  field  of  forces  is  converted  (2.4)  to  the 
fora 

*•>  JJ-^r (25) 

w 

where  l(r,  t)  - the  referred  to  p0  vector  of  the  surface  momentum 
density  of  the  external  forces,  which  affected  up  to  tor que/moment  t 
the  gas  at  the  points  of  surface  W. 

3.  Let  us  examine  infinite  quiescent  gas  in  which  with  certain 
moment  of  time  begins  to  move  lifting  surface.  As  has  already  been 
spoken,  let  us  consider  that  the  pcints  of  lifting  surface  little  are 
distant  from  certain  stationary  permeable  sirface  of  W.  Designating 
through  H0(t)  region  on  surface  of  W for  which  is  design/projected 
(along  standards)  at  torque/moment  t the  lifting  surface,  let  us 
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identify  this  region  with  quite  lifting  suriace,  transferring  to  it 
the  points  of  appl icat icn/appendix  to  gas  of  external  forces. 

Page  2 0. 

The  part  of  the  boundary  of  the  lifting  surface  whose  points 
have  different  from  zerc  normal  toward  boundary  conposing  the  rates 
cf  motion,  directed  from  region  W0,  let  us  designate  through  L,(t), 
and  inside  region  - L2  (t)  let  us  respectively  call  these  curves 
leading  and  trailing  edges.  Region  W0  can  te  limited  to  the  pieces  of 
curves,  not  having  the  normal  to  them  comprising  rate  of  motion  (by 
flank  edges) . 

Let  us  describe  the  motion  of  lifting  surface  by  the  parametric 
equation  r = r(u,  v)  of  surface  U by  the  position  of  the  curves  L, 
and  L2  at  different  moment  of  time.  For  convenience  we  consider  L, 
and  L2  at  the  initial  moment  coinciding,  but  in  order  not  to  exclude 
the  case  of  the  instantaneous  emergence  (or  disappearance)  of  the 
section  of  lifting  surface,  let  us  allow/assume  the  motion  of  these 
lines  with  infinite  velocities.  H,  (t)  and  fc2  (t)  - region  on  surface 
of  «,  described  by  the  curves  L»  and  L2  up  to  the  torque/moment  of 
time  t,  but  rt(r),  r2(r)  - the  tor que/moments  of  the  time  when  the 
curves  L*  and  L2  pass  above  the  point  (r)  on  surface  of  V. 
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Let  us  realize  that  the  points  of  regions  W,  and  W2  at  the 
■oient  of  the  passage  above  them  of  lines  L!  and  L2  continuously  emit 
the  acoustic  waves,  which  spherically  diverge  at  a rate  of  c.  Let  us 
fix  certain  point  of  space  (r0)  , not  locating  cn  surface  of  W,  and 
the  torque/moment  of  ti«e  t0.  The  parts  of  regions  Wt  of  andes  W2, 
frcB  which  had  tine  to  reach  point  (r0)  up  to  torque/moment  t0 
acoustic  waves,  let  us  call  the  audible  forms  cf  these  regions  and 

designate  respectively  W\  and  U/j 

The  condition  of  the  determination  of  point  (r)  in  region  W\ 

(cr  w’2)  will  be  »»', (r...  r,  (or  •l\(rUI  r,  f„)  <0),  w here 

T(r0,  r,  tu)  — I rv  — r — r |/„  — t (/>/  and  T »lr,  when  t = t,  (/  = 1,  2).  The  region, 

which  supplements  W]  to  let  us  designate  Wo  Region  W]  can  be 

named  in  an  audible  manner  of  the  film  cf  eddy/vortices,  while  region 
U/*  - audibly  of  lifting  surface.  We  will  be  restricted  to  the 
regular  case  when  in  the  composition  of  boundaries  of  the  region  W\ 
and  W ’•  enter  the  lines,  determined  by  equalities  * 

'«■, (r„.  r,  tu)  0,  V,(r0.  r,  /.)  0.  ( < U 

FOOTNOTE  * . In  the  special  cases,  further  net  examine/ cons idered, 
equalities  (3.  1)  can  occur  for  the  totality  of  the  points,  which  form 
region.  ENDFOOTNOTE. 
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Let  us  designate  the  lines  indicated  respectively  L\,  l‘t  and 
will  call  then  the  audible  forms  of  leading  and  trailing  edges. 
Region  W\  can  be  also  ]i»ited  by  the  secticn  cf  boundary  of  the 
region  H which  let  us  designate  through  /.j*  (Fig.  1).  Since  the  line 
£■!  has  different  from  zero  normal  rates  of  motion  its  intersection 
with  curved  /.,  for  points  (r0)  , which  do  net  lie  on  w,  it  is 
impossible.  Region  W',  besides  lines  /.’  and  Z.J  can  be  limited  to 
the  sections  of  boundary  of  the  region  H whese  totality  let  us 
designate  /»'  (see  Fig.  1). 

Page  21. 


0 


, 

, 


4.  *(r,  t)  - referred  to  p0  vector  of  aensity  of  surface  forces, 

which  act  on  gas  from  the  side  of  points  of  region  w„,  in  direction 

HP 

of  standard  N to  surface  of  W.  Regarding  the  vector  of  the  surface 
momentum  density  of  the  forces 

t 

r ( r , t)  j(r,  t)dt, 

0 

and  the  vector  of  momentum/impulse/pulse-pctent ial  of  the  flow, 
caused  by  these  forces,  accordingly  (2.5)  will  he 

-»  - .•  ,/s  -*  / -*  / \ 

t)- 

Per  the  determination  cf  velocity  potential  *,  it  is  necessary 

•? 

tc  know  derivatives  of  vector  E on  the  coordinates  of  point  (r0) . 


L 
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Let  us  exaaine  the  differential  6E  of  vector  E,  calculated  when 
joint  (r„)  is  displaced  in  the  direction  of  artitrary  unit  vector 
v’  in  distance  6h,  and  time  t0  grow/rises  ty  6t.  Since  region  W' 
will  obtain  during  this  variation  certain  increase  f>W\,  then 
introducing  derivatives  under  integral  sign,  let  us  find 


.(f&'p*- ',)(»-? jiW'. 


. (4. 


Through  V,(r,  t)  is  designated  the  tine  derivative,  a 


indicates  differentiation  ir,  the  sense  of  the  vector 


vith  respect 


r 
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/.*’  displaced  during  variation  position  cf  line  /*.  a r j-5r 
and  r radius- vectors  of  the  arbitrary  infinitely  close  points  of 
these  lines  (Fig.  2).  Regarding  lines  7_*  and  /.^  we  have  the 
egualities 

‘l',  (r0,  r,  t,t)  0,  y,(rn  t vt h,  r + ir,  f0-)-'7)=0, 


difference  in  which  within  Unit  when  Sr  — o,  It-*  0 gives  the  following 
condition,  superimposed  on  l’r: 


T^vn'l'i  t lr  VM‘i 


<Jtti 


where  Vo  and  v - Haailtor’s  operators  on  alternating/variable  r0 


and  r 
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Since 


Vo  V,  (rB  r)  1 r*  — r|,  ?V,=(70  - r) ' 


c Grad 


this  condition  can  be  written  in  the  foci 

(, e + (i)8r  — (t>-v)8/t  — cS/  =0,  - h 


suitable  both  for  a curve  /.*,  and  for  /-?,  moreover 

. _ . . 

c — (r  — rj  r — ra',  |x^=rGradt.  yjthe  particular  case  Zh=lt^-0  equality  (4.2) 
gives  the  following  condition,  superimposed  on  vector  dL  of  the 
cell/eleaent  of  curves  l\  or  Li- 

dL  (t"+  i‘)  =0.  (4-3^ 

Let  us  introduce  in  curves  L\  and  /.*  the  families  of  the 
alter  rat ing/variable  vectors  a,  which  lie  at  tangent  to  surface  of  w 
cf  plane,  but  clear  these  curves.  Assuaing  that  the  tangent  toward 
surface  of  H component  of  vector  Sr  is  directed  along  the  appropriate 
vector  a,  we  will  obtain  frca  (4.2)  with  an  accuracy  to  infinitesimal 
first  crder 

. • (e-v)8/i  4 cU  • (4.4) 

<>r  -=  . * * 1,1 

m ( e + *0 

The  eleaent  of  area  dS',  constructed  cn  vectors  Sr  and  dL,  we 
consider  positive,  if  Sr  it  is  directed  f rca  region  W,  (i  — 1, 2),  and  dL 
with  view  along  standard  H has  region  W‘  tc  the  left  of  its 
direction.  Thun  dS~  — \Slrdi\  and  on  (4.4)  we  find 
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|;V  null.  I 
m (e  i - !*•) 


[(e-  >)lh  | 


clt  |. 


During  replacement  curved  /.’  by  broken  line  from  cel  1/el e men ts 
d?  the  corresponding  areas  ds  • fori  the  region  i'  W\,  which  differs  from 
„\i\  to  the  small  second-order  quantity  relative  to  6h  and  6t. 
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Producing  in  (4.1)  the  replacement  of  region  on  «'  W*  (that  it  will 

net  influence  the  differential  IE),  me  mill  cfctain 


l \ 


. r r as  d / 1 \ r,  / - . I 

«=»JJ  *;(t Tp 

r C (IS  -/*  . / \/.,  Ih  dl 

■'  JJ  -*a f)(“-  , 


[Nmdj\ 

4tt Im  (e  (- 


c ) 


|(£-  > )lh 


cU\. 


The 

■* 

vector  E 


coefficient  of  6h  in  this  expression  gives 
on  the  position  of  point  UD)  in  direction 


derivative  of 
v,  i.  e . r 
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dj: 

d> 


r | NmdL  ](<’•  ') 
*1  r.lm  (e  -|  ji) 


r, 


i \ i m 

e ) c <)<„ 


(4.5) 


Since  * (?,  t)  = 0 with  t < rt  <r),  value  l’/r, t0  — to  curve  /.”  is 
equal  to  zero  and  the  third  integral  in  (4.5)  is  absent.  Calculating 
with  the  aid  of  (4.5)  derivatives  alcng  coordinate  axes  and 
determining  div  E,  we  mill  cbtain  velocity  potential 


f (r O'  (") 


rr  dS 

d ( 1 \ 

V r 1 

1 \ 

JJ  4* 

dNn  { l ) 

1 1 r>  '« 

‘J 

r r <*s 

0 t 

! J 4 rx 

\ 

r, 


c / 


(4.0) 


here  I’  and  # - projection  of  vectors  ■ ’ and  a on  the  going  along 

f) 

them  standard  N , a indicates  differentiation  with  respect  to  the 

position  of  point  (r0)  in  direction  N. 


According  to  (4.6)  in  all  points  of  space,  beyond 
exception/eliainat ion  which  are  located  on  surface  u/,\  potential  0 
exists  and  it  is  continuous,  and  gradient  f determines  the  rates  of 
flow  of  gas. 
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Potential  juap  * with  approach  frox  the  different  sides  of  surface  W 
to  the  located  on  it  point  (r)  on  the  basis  of  the  maximum  properties 
of  the  potential  of  dual  layer  is  equal  to  I (r,  t)\  that  means 
r ( r,  <>- circulation  on  duct/ccntour  G,  which  pierces  surface  of  W at 
this  point  (positive  1 correspond  to  the  intersection  of  surface  of 
H in  the  circuit/bypass  of  duct/contour  G in  the  direction  of 
standard  N) . In  the  case  of  c * « equality  (4.6)  transfer/converts 
into  known  expression  for  the  velocity  potential  of  lifting  surface 
in  the  incoa pressi ble  fluid. 

5.  Assuming  that  function  *(r,  t)  is  differentiated  in  all 
points  of  lifting  surface  and  it  is  final  cn  its  boundaries,  let  us 
determine  derivatives  of  potential  * with  respect  to  coordinates  and 
time.  For  this,  preliminarily  let  us  find  differential  8<p>  when  point 
(r0)  is  displaced  by  vector  and  time  increases  on  6t.  Designating 

through  >>  W*  and  '-W'*  the  increases  of  regions  W\  and  luring  this 
variation,  we  can  write 
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The  entering  here  integrals  on  regions  and  are  converted 

into  integrals  on  lines  /.*  and  /.'  analogous  with  that,  as  this  was 
nade  during  calculation  6E.  Fulfilling  this  conversion  and  taking 
into  account  that  to  curve  /.*  value  1' ( r,  tt—  turns  into  zero,  we 
obtain 


if  & h 


• r ds  d l I 

4*  dNa 


'•  ' )*'  |-  Jw 


r r ds  .a, 

4re  dN„ 


fJi  -| 


where  In  - the  part  of  the  boundary  of  the  region  UPn,  which  goes 
along  lines  l’  and  /.j,  ecreover  positive  circuit/bypass  i'0  with 
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view  along  N aust  leave  region  to  the  left. 
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o? 


Coefficients  of  6h  and  fit  in  (5.1)  essence  derivatives  . and 

r)?  , - 

jf  respectively.  Calculating  with  the  aid  cf  (5.1)  gas  velocity  v 

pressure  p,  we  will  obtain 

Kr^ii';  -r£;(  j-)|  ‘ '.) 


IT  J = COflSt 

e(e,\ ) , 


i | e(el\)  , /- 

+ JJ  4wr»/  *''•  /o 


/ 


/ W/n(e  + ,»)  ^ ’ C )' 


/,(r#-  'o)~  P°JJ  4V^vr(  / )N^  '•  r)+  r 3'(r;  c 

IT* 

0 

+ f | r.-  v.i  s..;.//:.  / 

•'  4«/m  (<?  i |i)  \ c ) 


+ 


(5.3) 


The  first  integral  in  formula  (5.2)  can  be  interpreted  what, 
velocity  in  point  (r0)  of  the  flow  of  the  incompressible  fluid, 
caused  by  the  dipoles,  distributed  over  surface  IT*  with  density 
r*(r)«=  l’(r,  t„  — * ) -\-  ^ i(r,  c -)  and  oriented  alcng  the  normal.  Such 

dipoles  produce  the  same  velocities,  as  covering  this  surface 
eddy/vortices  with  circulation  r*(V)  [function  r*(r)  depends  besides 
r cn  the  tcrgue/moment  of  time  t0  and  of  the  position  of  point  (r„)|. 
Let  us  further  call  i*  audible  circulation.  The  corresponding  to 
audible  circulation  vortex  systea  consists  cf  layer  on  surface  UP'* 
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with  the  vector  of  surface  intensity  ■[*  -/VyOradl'*  and  t he  going 
along  lines  l.\’  and  l.o  discrete  eddy/ vertices,  intensity  AG ♦ of 
which  is  equal  tc  I * on  and  t0  ^ jr  1 cn  Ln.  Utilizing  for 

calculating  first  tern  (5.2)  Biot-Savart* s forwula,  we  will  obtain 
the  following  resultant  expression  for  a gas  velocity  through  the 
intensities  of  the  eddy/vcrtices: 


(r,  t„  1 \dS  f HNi)\NmdL\ 


A!' 


’k.  f» 


Page  26. 

Expressing  eleaent  of  area  dS,  standard  N and  surface  gradient 
froa  l*  by  scans  of  curvilineac  coordinates  u,  v 

dS  = xdudv,  .V«|r'xrj*  1 (*  = l r‘Xf*t!), 

Grad  ! • - * 1 <K,  X /V|  l'l  - | rv  X \\  I',*). 

we  can  present  the  obtained  expression  for  velocity  also  in  the 
following  fora: 


\r.l-  a(u,  V) 


e (Ne) 
Ulc 3 


L O+'l 


r 
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here  A ^ rue ct  u,  H = r„e  + civ,  tut  value  a and  t essence  the  coefficients 
cf  vector  m on  basis  ru,  ru.  As  a result  of  (4.3)  to  curve  occurs 

equality  Adu  ♦ Bdv  = 0,  that  ensures  the  independence  of  the 
determined  by  formula  (5.5)  velocities  from  the  selection  of 
coefficients  a and  b.  That  entering  in  (5.5)  jacobian  from  vector 
r(u,  v)  and  scalar  v)  is  easily  determined  by  direct 

differentiation  of  the  function  !’*(«,  v): 


•Hr,  I’*) 
0 (u,  v) 


d(r,  V (r,  t)i  , I d(r,  a(r,  ft) 

0(u,  v)  , r c d(u,  v) 

t-t»  - - 


In  expressions  (5.4)  and  (5.5)  it  is  possible  to 
transfer/convert  to  limit,  fixing  point  (?0)  tc  surface  of  W.  The 
"direct/straight"  value  of  velocity,  equal  to  the  half-sum  of  the 
limiting  values  with  approach  to  point  (r„)  from  the  different  sides 
H,  will  be  determined  by  these  expressions  (improper  integrals  are 
taken  in  the  sense  of  principal  values) . In  the  most  important  for 
practice  case  when  W there  is  a plane  and  point  (r0)  it  lie/rests  on 
it,  the  third  and  fourth  integrals  in  (5.4)  disappear. 


Formula  (5.4)  solves  stated  problem  of  the  rates  of  flow  of  gas, 
caused  by  the  arbitrarily  driving/moving  in  it  lifting  surface  with 
the  assigned/prescribed  final  surface  load.  It  expresses 
Biot-Savart’s  law,  generali2ed  in  the  case  of  the  compressed  medium. 
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According  to  this  law  durinc  the  unsteady  icticn  of  lifting  surface 
in  gas  Biot-Savart 's  formula  it  is  required  to  apply  only  to  that 
vertex  elements  from  which  had  time  to  reach  the  point  of 
application/appendix  scund  signals.  Besides  in  addition  to  this 
appear  the  supplementary  "wave"  which  comprise  of  velocities,  which 
decrease  as  first  degree  of  distance  of  ed dy/vcrtices. 

6.  Example  1.  To  the  quiescent  gas  frex  tee  torque/moment  of 
time  r = 0 at  the  points  of  plane  = 0,  cctes  into  action  the  field 
of  the  directed  along  axle/axis  Oz  forces  with  constant  surface 
density  a = p.  Let  us  find  the  gas  velocity  in  torque/moment  t0  in 
the  point  of  axle/axis  Oz  with  coordinate  z0. 

Page  27. 

Auditory  sensation  area  r*  will  be  determined  b inequality 
'i  t-ct,,*' o and  will  represent  by  itself  circle  witn  a radius  of 
u Vc'tl  zl  with  center  in  the  beginning  of  coordinates.  The  vector 
of  momentum  density  is  directed  alcng  axle/axis  Oz  and  along  value  is 
equal  to  i (i)^pt.  Hence 


To  this  circulation  corresponds  circular  eddy/vortex  with  a 
radius  of  H with  intensity  - pt0.  It  will  excite  at  point  (z0)  the 
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directed  along  axis  Oz  velocity  r't  Furthermore,  in  the 

expression  of  gas  velocity  accordingly  (5.4)  will  enter  component 


where  integration  it  is  conducted 
B,  which  limits  auditory  sensation 
fl  between  the  radi us- vector  of  the 
the  axle/axis  Ox  and  assuming  that 
we  will  obtain 


e(Ni')  \Nmtll.\ 
4n/rni  (e  \ p ) 


(6. 1 ) 


in  circle  with  radius  of  with  a of 
area.  Introducing  vectorial  angle 
points  of  this  circumference  and 
vectors  t are  directed  radially. 


\Nm<tL  | — AV/0,  Nr  — — cos  a,  me  sun, 


(6.2) 


where  a - a half-apex  angle  of 
and  apex/vertex  in  point  (z0). 
(6.1)  elementary  is  integrated 


•i  i'. 


the  cone  with  circle  V* 
Taking  into  account  (6.2) 
and  gives  v"  = , whence 

z r2 
1 *0 

r Zu 

■Jc  'ibi1  <°>  c ■ 

<« 

0 n pii  /, 


as  basis/base 
integral 


Key:  (1)  . with. 


This  result  can  be  obtained  by  another  way,  for  example  by  the 
solution  of  one-di mensicna  1 problem. 


Example  of  2.  On  the  plane  z = 0 Cartesian  coordinate  system, 
evenly  moves  from  infinity  the  carrying  band  of  final  width.  From  the 
side  of  band  to  gas,  acts  the  field  of  those  directed  conversely  of 
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axle/axis  Oz  of  the  forces,  which  have  constant  surface  density  (*■•/» 

It  is  required  to  determine  the  gas  velocity  at  certain  moment  of 
time  t0  in  point  Q,  which  is  located  on  band.  It  is  arranged 
axle/axis  in  such  a way  that  up  to  the  tor que/mcment  t0  the  forward 
edge  of  band  would  coincide  with  axle/axis  Cy,  and  the  positive  part 
of  axle/axis  Ox  traversed  point  C-  The  functions  rt  and  r2 , which 
show  the  transit  time  of  the  leading  and  trailing  edges  of  the  band 
above  the  point  of  plane  z - C with  coordinate  x > 0,  will  be 


i. 


. (* 
v and  - /0  - y 


where  b - width  of  band,  v - a rate  of  its  motion. 

The  vector  of  momentum  density  for  this  flow  will  be 
r *•*</  -*.)/»  a i with  r,  < t < r2 , t*  = r 2 with  t > r2)  . Areas  Wj  and  w2 
are  half-plane  z = 0,  x > 0 and  z = C,  x > t.  Lines  /.*  and  /. j are 
determined  by  equalities  >r,  / c.*  o and  1%^/  ' (l  “ h)  o;i.e.  represent 

by  itself  conic  sections  with  focus  at  point  Q,  eccentricity  p = c/V 
and  directrices  x = 0 and  x = b. 


Page  28. 

Let  us  examine  the  case  of  supersonic  speed  of  motion,  when  p > 
1,  line  *1  - ellipse,  and  line  /. j it  is  atsent.  Areas  and  #"* 
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ccincide  and  represent  ty  themselves  the  interior  of  ellipse  (Fig. 

3)  . Audible  within  ellipse  circulation  i*  . The  corresponding  to 
this  circulation  vortex  system  consists  of  the  layer  of  those  cover 
the  interior  of  ellipse  of  parallel  to  axle/axis  Oy  eddy/vortices 
with  a density  of  p/V  and  the  concentrated  «ddy/vortex  with  an 
intensity  of  px/V,  that  goes  over  ellipse.  Since  (e  N)  =0  for  point 
Q,  the  rate  of  flow  of  gas  accordingly  (5.4)  are  wholly  determined  by 
Biot-Savart's  formula,  used  to  the  vortex  system  pointed  out  above. 
Introducing  the  vectorial  angle  between  a fecal  radius  and  the 
axle/axis  Ox  for  the  directed  along  axle/axis  Oz  velocities  from  the 
rectilinear  and  elliptical  eddy/vortex,  we  will  obtain 

'I  u 

i I p sin  (I  dx 

1 * J 2n  V (x  — xn)  ' 

0—i i 

0-0 

• f Px 

vp  (co*  0 dy  sin  0 dx),  (6.3) 

0 Si  ft 


where  x0  - coordinate  of  point  Q. 


The  first  of  the  integrals  is  undertaken  in  the  sense  of 
principal  value.  Coordinates  x and  y of  the  points  of  ellipse  and  the 
lengths  of  a focal  radius  z are  connected  with  0 relationship/ratios 
x = Hi,  y — I sin  9.  l = 

w hence  Vl  = — --  f ~ - . This  result 


. - /'( M V'm*-!) 

(6.3)  gives  , v 


(M  ^cosd)  • H * T/c.  Calculation  of  integrals 

/»  \ M>  I 


is  well  known  from  the  theory  of  fine/thin  airf cil/prof ile  in  the 
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supersonic  flow. 


Pig.  3. 
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FLCH  ABOUND  DELTA  WING  BJ  HKPERSCNIC  FLCH. 

V.  P.  Kolgan 

Hock  examines  the  problem  of  the  flow  aiound  delta  ucnj  of 
hypersonic  flow  of  gas  at  lew  angle  of  attack.  Is  used  the  method  of 
the  "sources"  of  the  pressure,  with  the  aid  of  which  th->  problem  came 
to  singular  equation.  There  is  carried  cut  tegu lar izat  ion  of  ♦•his 
equation,  because  of  which  obtained  integral  equation  with  continuous 
nucleus.  The  results  of  the  work  are  illustrated  by  examples  of  th** 
calculations  for  pressures  for  the  disturbed  zone  of  flow. 

This  problem  has  already  been  examined  by  a number  of  the 
authors  [1],  [2],  who  applied  for  the  determination  of  solution  the 
method  of  expansion/decc mpesition  in  series,  re f lect/r epresert ing 
physical  flow  plane  to  certain  fictitious  plane.  In  this  article  is 
proposed  another  approach  tc  this  problem,  which  ensures  obtaining 
solution  immediately  in  physical  plane.  Method  can  render/show  useful 
for  the  solution  of  a series  of  other  problems. 


1.  Pormulation  of  problem.  Let  us  examine  the  lower  surface  of 
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the  delta  flat/plane  wing,  arrange/located  at  low  angle  of  attack  a 
<<  1 ia  the  supersonic  flow  of  perfect  gas  with  large  mach  number 

li  2 7 | - angle,  formed  by  leading  wing  edges.  Let  us  assume  also 
that  the  parameters  M*.  '<  and  yi  satisfy  the  following 
relat ionsh ip/ratios: 

i—  I;  M x-j,»  I.  (Ill 

let  us  search  for  asymptotic  solution  for  a flow  around  of  the 
wing  under  the  made  assumptions.  It  began  the  systems  of  coordinates 
iyz  was  arranged  in  the  spout  of  wing  sc  that  the  plane  xz  would 
coincide  with  wing  plane,  and  X-axis  coincided  with  the  axis  of  * he 
symmetry  of  wing.  Let  us  introduce  the  dimensionless  unknown  velocity 
functions,  pressure  and  density 

V (/„<«/».  v\  W);  />  • iv  :■  (1 -> 

where  i and  f.  velocity  and  the  density  of  flow  in  the 
undisturbed  flow.  By  virtue  of  assumption  (1.1)  leading  winy  edge 
will  be  supersonic,  shock  wave  - plane  and  the  parameters  of  flow 
after  it  by  constants  up  to  the  disturbance  cone,  proceeding  from 
spout  wing.  This  zone  of  flew  with  the  constant  parameters  is 
designated  by  index  1 (Fig.  1)  and  it  will  te  further  called 
exterior.  The  zone  of  flow,  which  lies  within  Bach  cone,  which 
emerges  from  the  spout  of  wing,  let  us  call/naie  internal  area  (area 
2,  Fig.  1)  . 
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Let  us  write 

o ut 

the  solution 

for  an 

exterior.  By  virtue 

of  the 

made  assumptions 

solution  can  be  presented 

in  the  form 

a * 

P* 

a5/’,,  0(»‘);  j 

V* 

o (*>)\ 

h (1 

■'*) 

w* 

<x-  w„  -j-  O (a1); 

Y* (x,  z ) 

a>u(x,  z)  ■ 0(a>), 

where  y = T*  (x. 

z)  - 

a surface  of 

shock  wave,  and  values  with 

zero 

indices  remain  the  order  of  one  with  tendency  a toward  zero. 


Utilizing  relationship/ratios  on  obligue  shock  wave,  it  is 
possible  to  obtain  the  following  solution  ter  the  exterior: 


>'„  A0x  i B0z\  p9  A„  i 1 + — ; - 


Uu  — "n"  A 0 1 lJn 


1 4 A„ 

A, 


; w0 


A, 


(1.4) 


~ 1 1 

where  Av=-  , 1 


IW 


'iMV,  1 
mL«! 


B — - - " 

" , 


> = ‘k  t.; 

-f— adabatic  index 
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Let  us  examine  interior.  Let  us  introduce  the  new  unknown 
functions  and  new  variables  with  the  aid  of  the  relationship/ratios 


u*  1 | o2w„  I 0*U.-\  0(aM;  p*  p„  4-  <xR,  + O (a2); 


v"  a*  V,  + 0(as); 

/it),  | 0(o:l); 

/>*  **p»  I 4- G(«4); 


Y*  aA„ X 4-  a2  Y,(x,  z*) 
z * z/a; 

V*  V a 


I 
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Fig-  1. 


Pig-  2. 
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Then  for  interior  cf  the  equation  of  ges  dynamics  in  this 
approach/approximation  they  seem  in  the  f on 


dU , , r) 


Ox 


Ox 


dv,  d r pa 

Ox  1 dy*  V,  pi, ) 


V P» ) 

d 

Ox 

1 

V 

w 

1 

,,  o , 

< Pi\ 

a/'o'V'i  , <*«’ 

-05 

V Po  ) 

1 r dz:| 

Ov\  0 / p{ 
Ox  Oz*  { p„ 


0;  a\  - t 

Pi, 


(1.6) 


Boundary  conditions  on  shock  nave  for  this  problem  when  y*  -=  Anx 


take  the  form 
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U i I A„  ^ i x*  P\  — ■ | j (1  H A,)  A «'< 


(TH-  -M<>)3 


1 » 


•!  A,,)2 


If  we  introduce  conical  variables  Y—y  and  ? - and  to 

<I„  X II ii  x 

exclude  froa  equations  (1.6)  U,,  V,,  ut,  then  problea  it  is  reduced 


tc  cne  equation  relative  to  pt 


0 ,l"’ 


under  following  boundary  conditions  (Fig.  2): 


on  the  shock  wave  AE 


y = M A"  ; 0 Xt  k (k  171  M'); 

«0 

P\Y  1 j/^O  , 1 A>  '}  72  7 1 1 

/>i/  /<'-'«„/  ° V «?,  lA.(A,  ! 1)  7 1 


A,Z  H;  ; 


on  the  arc  of  unit  circle  BC 


Y |1  Z\  k^Z<A\  /?,  0; 


on  the  surface  of  wing  CO 


>’  0;  0 F t 1;  P\  i — 


(111) 


DOC  = 78068002  PAGE 

75 

on  the  plane  of  symmetry  AO 

Z - 0;  O • Y < fA\  pt  / - 0.  H.J2) 

Function  p,  can  be  analytically  continued  to  entire  polygon  BDEP 
symmetrically,  and,  if  we  starch  for  the  analytical  solution  p,  (Y, 

2)  in  the  fora  of  the  function,  even  relative  to  its  arguments,  then 
conditions  (1.11)  and  (1.12)  will  be  satisfied  automatically. 
Condition  (1.9)  is  spread  to  AF  and  ED  ty  symmetrical  fora,  and  on 
the  arcs  of  the  unit  circle,  is  retained  as  before  condition  p,  = 0. 

Page  32. 


2.  Information  of  boundary-value  problem  to  integral  equation, 
let  us  search  for  pi  in  the  form  of  the  contour  integral 

/>,  - U('t,  Dfi  y.  Z,  ’.)  ds,  (2  I) 


where 


f(y.z,  % g 


i |(i  y ,,  zg»  (i  > ■ -zj)(i-  v ;*n  • 

l/i  1 1 Kt,  z;  id  y1  z- mi  V 


C)  - the  intensity  of  the  "sources"  cl  pressure, 
arrange/located  in  a symmetrical  manner  on  IF  and  ED; 


s - arc  length. 
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The  solution  of  form  (2.1)  satisfies  equation  (1.8),  elliptical 
within  unit  circle,  and  bcurdary  tc  ccnditicns  ( 1.  10) - ( 1 . 1 2)  . For 
this,  in  order  to  satisfy  condition  (1.9)  , let  us  find  the  limiting 
values  of  derivatives  P\*  and  P'  * with  approach  to  shock  wave  of  from 
within  area  BDEF: 


Pt  H M - 0,  j nil,,  Z,  T(,  \)da 


*n  (M, 
k a 


Z) 


p,  /,  (M  0,  Z)  J g(n,  \)y/  (M,  Z,  Tlt  ' ) </ s 

Substituting  expression  (2.2)  under  condition  on  shock  wave 
(1.9),  we  will  obtain  following  singular  function  g (Z)  : 


* k 

. ■/ 1/  tp-  zz  . , r i7  i r 

' A'(Z)  H#i  ) J r / "M 

-A  -A 


«1»  (Z,  '.) 


A.,  A | MA'*, 

MZ|Z(Z-fC)-  2 1 |Z(1  M5)  CA-'|(«,  -4,Z») 

(Z2  j-7J)/e-  I 1 M-  2ZC(I  I M2) 


(2.3) 


Equation  (2.3)  is  related  tc  the  class  of  homogeneous  complete 
special  integral  equations  for  open  circuits. 


Let  us  search  for  function  g(Z)  in  the  class  of  bounded 
functions.  Then  pressure  pI#  calculated  according  to  formula  (2.1), 
will  be  bounded  function. 


Let  us  write  characteristic  equation  fer  equation  (2.3): 


*Z  | k‘  Z* 
A* 


f(Z)  i (/!,  A^J 


(-’.4) 


DOC  = 78068002 


PAGE 

77 

The  index  of  equation  <2.  4)  in  the  class  cf  bounded  functions  x 
= 1,  and  the  limited  solution  of  equation  (2.4)  is  record/written 
with  an  accuracy  to  constant  factor  in  the  form 


l(Z) 


Ht) 


f{Z)  \' k2  — Z.1  b (Z)  exp  I’  (Z)\ 
k-(B,  — A . Za) 

I ZHV-ZZ)  — A , Z3)-|'  ■ 


**  f In  I g*  A'-  '2A.,x‘-[k-  xi)+(Bl-Atxi)(kt  2jc*). 

* J h A:-  x*  -va  (A:1  - jc»j  + A;1  (#,  — A.t  x-f  dx 


rJ) 
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Let  us  conduct  now  for  equation  (2.3)  the  usual  procedure  of 
regularization  by  the  solution  of  characteristic  equation.  As  a 
result  we  will  obtain  the  regularized  equation 

* 

ff(ZH-/A'(Z.  :)*(',)  d:  f(Z),  (2.0) 

0 

where 

K (Z,  ;i  a (Z)c(Z)  '!>,  (Z,  1)  - 

k 

2h(Z)V  k'  Z-  exp  I"  (Z)  /'  re  (T)  ‘h,  (t, 

~ J I ' fe5  - t- Z*)exp  I'(t)  ’ 

u 

Z I k-  Z- 

’ I^M*2  Z ) | *'(»,  -d-.  Z^)-| 

A;2 

CiX>  r.\ZA(k'-  7S)  | At'(/y,  /\,Z')"|'»  ; 

<i>,(Z.  c)  = 'i«(z,  r,)-(-<i>(z,  :). 

When  deriving  the  equation  (2.6)  there  was  used  the  paiity  of 
function  g(Z).  From  the  theory  of  special  integral  equations  [31  it 
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follows  that  equation  (2.6)  represents  fcy  itself  the  integral 
second-order  of  Fredholm  equation  with  the  ccntinuous  nucleus  which 
already  can  be  solved  by  the  approximate  numerical  methods. 

Let  us  note  that  the  procedure  of  regularization  becomes 
inaccurate  when  Mvu  .00,  because  is  changed  the  index  of  equation 
(2.4),  and  coefficient  c(Z)and,  consequently,  a Iso  nucleus  K(Z,£) 
become  disruptive. 

3.  Obtaining  solution  and  its  standardization.  The  problem  of 
finding  the  distribution  cf  source  strength  g(Z)  came  to  second-order 
of  the  integral  equation  of  Fredholm  scluticn  with  continuous  nucleus 
(2-6).  For  finding  the  solution,  let  us  use  the  following  approximate 
diagram.  The  cut  of  the  integration  [0,  k]  let  us  divide  into  n equal 
parts  As  = k/n  and  on  each  section  As,  approximately  let  us  count  the 
function  z,  of  constant.  Then,  replacing  in  (2.6)  integral  by  the 
sum  and  varying  n times  variable  z,  we  obtain  for  the  system  of 

the  linear  algebraic  equations 

Ag=--f,  (3.1) 

where  the  matrix  element  A takes  the  fcrm 

(Iu~l<  (Z,,  Cy)  As  i 5(y; 

column  element  /,  /(Z,);  Z,  ~ the  coordinate  cf  the  middle  of  the  i 
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For  providing  the  necessary  accuracy /precision  during 
calculation  K[7i.'.j)  the  integral,  entering  the  expression  for  K(Z, 

S) , was  calculated  from  the  formula  of  rectangles  with  the  number  of 
separations  n*. 

Solution  g (Z)  , obtained  from  equation  (2.6),  is  determined  with 
an  accuracy  to  constant  factor  due  to  the  arbitrary  selection  of  the 
solution  of  the  characteristic  eguation  f(Z).  Let  us  note  that  from 
conditions  on  shock  wave  (1.7)  it  follows  that  along  shock  wave  is 
correct  the  relationship/ratio  between  w4  and  pt; 


<>P\ 

07. 


A 


, ..  0w{ 

' o Po  (,o  7 ,)y 


(3.2) 


Fage  3h. 


Since  during  motion  along  the  shock  wave  AE  function  w,  with 
continuous  form  changes  from  zero  to  , the  with  the  aid  of 

relationship/ratio  (3.2)  condition  of  the  standardization  of  function 
p,,  can  be  presented  in  the  form 


h 


II 


0 

07 


PiW-  , 


(3.3) 


Substituting  in  (3.3)  expression  (2.1),  we  will  obtain 
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standardization  condition  fcr  function  g(Z): 


(•  dz  r |*2(?  Z*)-AW\g<t)dZ 
Jl  k-  Z-  J l*5(Z*  + Cs)  4M*|2  • 4Z-'C*(1  M2)-  ' 


Ar 

1 f 


dZ 


' g(QdC  4 Ai',f;ot/0 


Z/T*  2*  J C-Z  i-t-'l 


(3.4) 


The  second  integral  in  equation  (3.4)  is  located  with  the  aid  of 
the  known  foraula  of  the  exchange  of  tie  order  of  integration  in  the 
dual  integrals: 


1 


* ft 

I "If 

-k  —k 
k 


nod: 

J?  (0) 

-Z*(C 

Z)  ' 2 ft 

k 

(' 

<1Z 

' J z , 

P Z) 

(3.5) 
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Fig.  3. 


Page  35. 


Last/latter  integral  in  equation  (3.5)  is  equal  to  zero,  since 
internal  integral  on  Z turns  into  zero.  Finally  the  condition  of  th» 
standardization  of  function  g(Z)  takes  the  fori 


k k 


0 


J. 


\kH7  - /:■)  - 4 M"!  d/. 

k1  /.  - \\k-(Zr  -j  - '-’)  * i M**M 

"lr(0)  4 A^,  <i„  ... 

'2  k I ).  1 


4.  Examples  of  calculations.  Employing  the 
was  comprised  the  program  of  calculation  for  cc 
calculations  was  set/assused  y - /.H  and  X = 
with  an  increase  in  the  nuxter  of  separations  ( 


given  above  procedure 
mputers.  During 
1.  It  turned  out  that 
n = 5;  10;  20;  40)  the 
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difference  between  the  consecutively  obtained  solutions  is  rapidly 
reduced;  for  example,  at  the  value  of  the  parameter  M,  a 1 the 
difference  in  the  values  of  pressure  pt  during  calculations  for  n = 

20  and  40  is  exhibited  only  in  the  fifth  sign.  By  this  is  confirmed 
the  effectiveness  of  the  selection  cf  the  calculation  method.  In  Fig. 
3-5,  are  constructed  the  isobars  P!  in  the  disturbed  area 
respectively  for  the  values  of  the  parameter  MOC«  = 0,1;  1;  5.  in  these 
cases  the  calculations  were  performed  with  t = 20.  Given  data  of 
calculations  are  confirmed  ty  the  results  cf  [2],  in  which  the 
solution  is  found  in  the  fora  cf  a series. 


Fig.  4.  Fig.  5. 
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AERODYNAMIC  INVESTIGATION  Of  ELEVONS  ON  LON-ASE ECT-RATIO  WINGS. 

V.  G.  Mikeladze. 

Are  examined  the  aerodynamic  characteristics  of  elevons  as 
crgan/controls  of  longitudinal  and  lateral  control  on 
lcw-aspect-ratio  wings.  Is  presented  the  method  of  calculation  of  the 
aercdynaaic  characteristics  of  elevons  at  sutscnic  and  supersonic 
speeds.  Are  given  the  results  of  systematic  studies  in  the  effect  of 
the  separate  parameters  of  elevons  on  their  aerodynamic 
characteristics.  Is  given  the  flow  pattern  cf  ] ow-aspect-r atio  wing 
with  the  deflected  to  large  angles  elevons  at  subsonic,  transonic  and 
supersonic  speeds.  ,• 

METHOD  OF  CALCULATION  OF  1HE  EFFECTIVENESS  CF  ELEVONS  AT  SUBSONIC 
SPEEDS. 

The  method  of  calculation  of  the  effectiveness  of  elevons  on  the 
wings  of  arbitrary  planfcrt  at  subscnic  speeds  is  instituted  on  the 
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use  of  a reciprocity  tbecrei1.  which  establishes 

ccBBunication/connection  between  aerodynamic  wing  characteristics  in 
direct/straight  and  return  flow  when  the  direction  of  velocity  of 
incident  flow  V0  is  replaced  by  reverse/inverse. 


FOOTNOTE  * . The  use  of  a reciprocity  theorei  for  the  evaluation  of 
the  control  effectiveness  was  suggested  by  A.  I.  Golubinskiy. 
ENDFOOTNOTE. 


The  derivatives  of  the  lift  coefficient  c"iB,  of  the  coefficient 
of  pitching  moment  m,/a  and  cf  the  i ol  li  rg- me  men  t coefficient  m ts»  in 
the  angle  of  deflection  of  elevon  8,„  can  be  presented  as  follows: 


/// . 


js ; 

(1) 

Vll 

(2) 

^ in 

s.  Jfp^.  dS-, 

H 

(3) 

here  S - an  area  of  wing; 
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p pressure-drop  coefficient 
of  wing,  driving/moving  at  angle  of 
turned  flow; 


between  lower  and  suction  sides 
attack  without  rotation,  in  the 


■> 


L 


j 


* 
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Carrying  vortex/eddy  surface  is  simulated  ty  connected  Vikhrev's 
series  cords.  Each  cord  is  replaced  ty  several  cblique  horseshoe 
vortices,  which  consist  of  the  bound  vortex  with  constant 
intensity/strength  along  spread/sccpe  !',  ard  free  vortices. 
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Page  38. 

Design  diagram  is  represented  in  Fig.  1.  Wing  v.ith  the  fracture  of 
leading  edge  in  the  examined  case  of  determining  the  effectiveness  of 
elevons  was  replaced  by  four  Vikhrev  by  the  cords  from  chord,  each  of 
which  Has,  in  turn,  was  replaced  ty  twelve  ctligue  horseshoe  vortices 
along  the  semirange  of  wing. 


Thus,  on  each  half  of  the  wing  it  is  arrange/located  of  48  bound 
vortices.  In  each  cell  of  the  formed  grid,  the  bound  vortex  coincides 
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with  the  line  of  1/4  chords  of  cell,  and  the  distance  between  free 
Vikhrev  by  cords  is  equal  to  the  spread/sccpe  of  this  cell.  The 
boundary  conditions  of  nonpassage  are  satisfied  for  each  cell  at  the 
point,  arrange/located  cn  the  Biddle  of  the  line  of  3/4  chords. 

During  the  calculation  of  the  effectiveness  of  elevons,  it  is 
necessary  to  have  values  of  the  dimensionless  circulation  of 
eddy/vortex  during  the  action  of  wing  at  angle  of  attack  without 
rotation  ?«, , dinensionless  circulation  during  the  rotation  of 

wing  of  relatively  axle/axis  Oz  and  di lensicnless  circulation  for 

the  wing,  which  rotates  around  axle/axis  Ox. 


Satisfaction  to  boundary  conditions  gives  three  independent 
systens  of  equations  which  make  it  possible  to  determine  three  the 
dimensionless  circulations  indicated  above 


V 

i i 


V 

i 1 


<U/V  t 

' AUV  -n 

<«%  i 

*“/ 

— 2ir  ; 

K 

2 it  ' 

(J  1 

> 2,  ■ • /«; 

m = Nn), 

(4> 


where  U/v  - diaensicnless  velocities,  caused  by  the  oblique  horseshoe 


vortex  i at  point  j. 


V 

V'n  I. 


x0/,  z0  - coordinate  of  the  junction/unit  calculation  point  j 
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at 


- diaensionless  circulation  of  eddy/vor tex;  /„  - distance 


between  free  vortices  (see  Fig 


S#  = — 

. I);1'  l,2\  z - dimensionless 


I '» 

I M - 


coordinates  of  point  H ; / - sweep  angle  of  the  bound  vortex; 

AU/>jy  - additional  diaensionless  velocities  which  appear  at  point  j 
frca  the  eddy/vortex,  which  is  located  on  left  half  wing  and 
syaaetrical  to  eddy/vortex  i, 

/>• 

Page  39. 

The  calculation  of  all  diaensionless  velocities  was  conducted  on 
the  formulas 


V -r 

U/  _ 


where 


, - (i#,  /), 


W^v(?n.  /)=!/,({,,  /)  | V/v(50.  C,,.  /), 


(5) 


(fi> 


C„,  /)-  the  velocity,  caused  by  the  bound  vortex  at  the  arbitrary 


point: 


/]  cos  / i Co  sin/ 


sin  / C-,  cos  /-[ 


1 

COS  / 


cos  / 


i tK'/.i'  t (I  ■ Cy)- 

C„sin/  ! C0  cos  / 


(7) 


I (5«  — in  /)“  1 

l/, (c0»  *•*  /)-  the  velocity,  caused  by  free  vortices  at  the  arbitrary 

pcint: 


( '0  * %l>>  / ) 


j r 


>n/ 


1 Co 


I (So  1 tR/)*  | (I  c,,)-  J 
- <n  / 

i (C„  ikz)2  id + w I ' 


(«) 


L 
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After  solving  system  of  equations  (4)  and  after  calculating 
values  p * , P-,  and  p«,x_.  we  determine  total  wing  characteristics  with 

the  deflected  elevens  on  the  formulas 

*1 

cj— */ 

Z, 

m*’"-  / | f/2; 

= >■  J •h„,, 


(12) 

(13) 

04) 

here  x* - coordinate  of  the  leading  edge  of  wing  section  in  the  turned 
flew;  b„  - relative  chord  of  eleven. 

Page  40. 


n 

where  \ = ~ - wing  aspect  ratio; 


y„  2<57.|  J ^ 


h»f  2.57  3 -1  t/x; 

' * ^ JH 

^ *• 

p,„x_dx, 


Vi, 


. '2-57.3 


(!») 

(10) 

(11) 


For  an  example  Fig.  2,  gives  the  results  cf  calculations  with 

( 

the  aid  of  the  computers  of  the  effectiveness  of  elevons  on  wing  with 
the  fracture  of  leading  edge.  The  ccupaiiscn  of  calculated  and 
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experimental  data  shows  that  computed  values  of  coefficients  c'y*“,  m1** 
and  are  higher  experi uerta 1.  Ore  cf  the  reasons  for  this 
disagreement  is  the  not  considered  fcy  theory  presence  of  the  slots 
which  are  formed  between  the  stationary  part  of  the  wing  and  the 
deflected  elevon.  The  effect  of  slot  on  the  effectiveness  of  elevons 
is  shewn  on  Fig.  3.  The  introducticn  of  the  empirical  coefficients  of 
k - 0.85  into  computed  values  of  derivatives  c'v<" , rn  •“  and  m\w  makes  it 
possible  to  obtain  the  values  of  these  derivatives  with  an  accuracy 
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Fig.  3. 


Key:  (1).  Ill  Version  of  elevon;  4 sections.  (2)  with 
partition/baffles.  (3)  without  partiticn/baffles. 
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METHOD  OF  CALCULATION  OF  EFFECTIVENESS  AND  HINGE  HOHENTS  OF  ELEVONS 
AT  SUPEHSONIC  SPEEDS. 

Effectiveness  of  elevons.  The  nethod  oi  calculation  of  the 
effectiveness  of  elevons  cn  low-aspect-rat ic  Mings  at  supersonic 
speeds  is  instituted  on  the  linear  theory  of  supersonic  flows.  During 
the  calculation  of  the  effectiveness  of  eleicns,  it  is  assuned  that 
the  rotational  axis  of  elevon  is  supersonic,  Mach  lines  do  not 
intersect  the  wing  chord  and  elevon. 

Elevon  with  the  adjacent  sections  of  wing  is  divide/nar ked  off 
into  zones  (Fig.  4) : zone  I is  linited  by  the  Each  lines,  which 
proceed  fron  the  point  of  intersection  of  rotational  axis  with  the 
root  chord  of  elevon,  and  the  trailing  edge  of  elevon;  zone  II  is 
linited  by  the  Hach  lines,  which  proceed  frcn  the  point  of 
intersection  of  rotational  axis  with  the  rcct  ard  end  chords  of 
elevon,  rotational  axis  and  the  trailing  edge  of  elevons;  zone  IU  is 
United  by  the  Hach  lines,  which  proceed  frcn  the  point  of 
intersection  of  rotational  axis  with  the  end  chord  of  elevon,  and  the 
trailing  edge  of  elevon. 


The  pressure-drop  coefficients  in  these  zones  are  egual  to 
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Pi' 


p\sr 


4 

«-V 

‘ JC 

(15) 

« l P* 

-t  K*y". 

1 

p\r 

4 

(l«) 

l I5’  — /,» 

4 

arccos 

-«* 

‘ X 

(17) 

: l fT 

i ' 

her*  Pi  - juap/drop  in  the  pressure  coefficient  between  the  lower 
and  upper  surfaces  of  elevon  in  the  i zone  during  the  deviation  of 


tK/. 

P 


elevon  of  1 rad;  a — S'"  (where  tg  /»•  - sweep  angle  along  the 


azis  of  the  rotatlos  of  elevon);  * VM- — I;  o„  - angle  of 
deflection  of  elevon,  determined  in  section  throughout  flow. 
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Fig.  4. 


Page  42. 


Foe  a low- aspect-ratio  wing,  including  for  wings  with 
alteraating/variable  sweepback  on  leading  ccge  and  the  trailing  edge, 
perpendicular  to  the  axle/axis  of  syaaetry,  expression  for 
derivatives  e*»\  and  m they  take  the  fera 


c !'.<»- 


4S,1ge 


57,3  S YfP  — tg2  / , 


m ••  — • 


0o\»  — !>«  ,»)  (V  1— a* — 1)  , 

tg  X,. 

, - I 

+ (*»  b*.  JM  l-mu) 

25,  /,  tg  £ 

2X..  > ' 

h 3 \ I , 

t ( hu  i n |~  h*.  jh)  /.h  Zo  „ ■ 


(18) 


f | 1 tl-  1 | .7  ' i .'l  s 

■ , - I - — (Oil  D K ,H  j iO  »» 

57,35/ igJX»  1 


-.'i  7 thu  ,h -,») 


1 

3 


m. 


57,3 


m. 


.)  4j ; 

(19) 

'/III  (^0  JB  ^ 

;■  i /*/2  1 

bA  \ 

, (20) 
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where 


* 2 b0i  Siig2*  (a  1 +]/l  a’Jfto®, 

m A"  ^ 


.S'  /JB  I p-  - tg*z. 


(21) 


m ’ft 


fto  A -Sa  tg2  e 


1 3ft.  b'Q\*-{b'n„-bi  p/;.  — 


,n  3 bA  Sp  Kp2-tgsx.„ 

‘ (*0  3B  Sll)2  ' ^0  3H  (^A)  >H  “I'  3/>..  3B)|; 


/« 


8 feto*$*tg*e  (fl  -1  |-  /1-a2  Ak. 


Alll 


bA  s tgz,.l'p*--tg*z 


(22) 

(23> 


/■  I H - 

S in  — 


4.S'a  /»kJ,b  t?_L (1  .{  rt—  /I  - <i2). 


Slg/,,  V'P5  — tg*Z. 


(24) 


In  these  foraulas:  c***  - lift  coefficient  of  wing  froa  the 
deviation  of  elevons  on  both  halves  of  wing  to  one  side  on  1°;  m\-«  - 

rolling-aoaent  coefficient  during  the  deviation  of  elevons  on  both 
halves  of  wing  to  opposite  sides  on  O®;  »»**•-  coefficient  of 

longitudinal  noaent  with  respect  to  z axis,  passing  through  the  point 
of  intersection  of  the  loot  chord  cf  elevon  with  its  rotational  axis 
(Pig.  4)  during  the  deviation  of  elevons  on  both  halves  of  wing  to 
one  side  on  1°;. 


Page  43. 


•Sa,  « - area,  root  chcrd,  span  and  seaiapex  angle  of  base  delta 


alag  (Fig.  4) ; 

bn  tit 

7,  a2  ’ 

to*,  tt  i 


bo .» — - 

b' 


7a  2 ’ 


: 

*0-— I7j2 


— relative  chord,  the  span  and  the  position  of  the 
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beginning  of  elevon,  in  reference  to  the  seairange  of  base  delta 
Ming. 


In  the  case  when  elevon  has  the  constant  absolute  chord  (z», =>0) 
of  expression  for  derivatives  4.  m\  and  m)  substantially  they 


are  sinplified: 


a,_  H VNA« e 

>•  ""  57,3  Sp 


i 2.S ± /j^3r  lan  /')  i * , 

= — g7  ;1  v/p  ■ — (2z0sb  I 

4.S'i  bn  4 tgJ  «/ 

w*  ~ 57,3  Sb a P 


Those  obtained  by  the  calculation  of  the  characteristic  of  the 
effectiveness  of  elevons  are  sere  experimental  values.  Processing  the 
results  of  experinental  data  and  their  comparison  with  the  results  of 
calculations  shows  that  for  the  evaluation  cf  the  effectiveness  of 
elevons  at  the  noderate  supersonic  velocities  cf  the  value  of 
derivatives  c*p,  m\'*,  /»*«,  determined  according  tc  formulas  (18) -(20)  it 
is  necessary  to  multiply  by  the  empirical  coefficient  of  K = 

0.  85-0.9. 


Hinge  moments  of  elevons.  The  calculation  of  the  hinge  moments 
of  elevons  at  supersonic  speeds  is  done  with  the  same  limitations, 
which  were  accepted  for  the  calculation  of  the  effectiveness  of 
elevons.  The  hinge-aoaent  coefficient  relative  to  the  rotational  axis 
of  elevon  during  its  deviation  on  1°  is  equal  tc 
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here  -S, 


»yt 


m 


in  i 


m)r  - '57^r  ■*  + mmni)’  (28> 

...  >*--•*'  _ I +v^7'-  (4 

3S„  bA  „ tg  /,» V p*  - tgv.M  l 2 v - * y 


arccos<i  _j_  1 — «‘ 


(2!») 


*7^“  Jfco  +h  | $ — 4&o  *h  1 ( ^0  iP  " bK  m) 


4 tg  * 

S.B  fr/t  .H  1 tg‘  /9 

(fr0  911  fr  k 9»)*  |/'9B  4-  (fr»  9B  fr«.  »»)1)  Ijs  /ill  ^.B  ^fl  »»  I^B 

fro  BB  ~{~  2 ( ^0  9B  frtt  9b)1  ^9B  (fro  9B  fr*  9B ) (9B  ~f“  (fro  9B  ~ fr*  9ll) 
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III  III 


bo  b k . vh  (b[)m  b ^ •■♦p ) 

(^n  9p 

fr«  9»  H l 

(30) 
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Hfrh;'.B  tg  s 
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W J J * 

1 </-'  ( 1 

30  >„  Ija  iU  tg  / ,„  | ;i*-  tg’  /,„ 
arecos  a 

! —ir 

2 

2 \ 2 

r-  J 

1 

(31) 

•S\.  . fr*,. 

Vi  , fr^  ..  ! J 


- relative  area  and  the  lean  aerodynamic 


chord  of  elevon  to  rotational  axis. 


Page  44. 


In  the  case  of  the  constant  atsolute  chord  of  elevon,  the 


hinge-aoaent  coefficieat  is  e^aal  to 


m 


57 


4tg>fr»  /j  4fr„  \ ..... 

3 sL*a~>  xm  Vj‘ 


As  shoved  the  cospariscn  of  the  calculated  and  experiaental 
values  of  hinge-aoaent  characteristics,  calculated  hinge-aoaent 
coefficients  prove  to  be  soaevhat  overstated.  Therefore  during  the 
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PIGE 

estimation  of  the  hinge  moments  of  elevens  at  the  aoderate  supersonic 
velocities  coaputed  value  cf  coefficient  .determined  in  formula 
(28) , nust  be  multiplied  by  the  empirical  coefficient  of  k 0.85. 

Physical  flow  pattern  of  wirg  with  the  deflected  elevons,  obtained  by 
the  method  of  pressure  distribution. 

The  studies  of  the  physical  flow  pattern  of  wing  with  the  ' 

deflected  elevons  were  conducted  on  the  model  of  low-aspect-ratio 
wing  with  alternating/variafcle  sweepback  on  leading  edge  (Pig.  5). 

The  analysis  of  the  effect  cf  each  number  of  the  incident  flow  on 

flow  around  of  the  wing  with  the  deflected  eleven  is  carried  out  ^ 

based  on  the  example  of  the  examination  of  the  air-load  distribution 

in  wing  section,  arrange/located  approximately  in  the  middle  the 

spread/scope  of  elevon. 

Fig.  6,  depicts  the  diagran/curve  of  pressure  distribution  along 
wing  chord  in  the  presence  of  the  deflected  elevon  with  Nach  number  = 

0.6. 


XT' 

If 


T 
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Key:  (1).  Section.  . 

Page  45. 

It  is  evident  that  the  deviation  of  elevon  causes  the  redistribution 
of  pressure  along  an  entice  ving  chord  both  on  upper  and  on  lover 
surface.  Hith  Bach  nunber  = 0.83,  deviation  of  elevon  of  negative 
angles  is  no  longer  caused  the  redistr ifcut icn  of  pressure  along  an 
entire  ving  chord  cn  lover  surface1. 

PCOTWOTE  1 . Analogous  results  on  a irfcil/prcf ile  vith  control  at 
transonic  speeds  were  obtained  by  G.  P.  Svishchev  into  1948. 
EHDPOOTHOTE. 


During  transition  to  snail  supersonic  velocities  (R  * 1.05,  Fig 
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7)  the  deviation  of  elevon  of  negative  angles  causes  the 
redistribution  of  pressure  on  lower  surface  only  along  the  chord  of 
elevon.  On  suction  side  of  wing,  the  zcne  of  the  effect  of  the 
deviation  of  elevon  of  pressure  distributicn  is  spread  forward,  the 
further,  the  greater  the  angle  of  deflecticn  of  elevon. 

The  analogous  phenomenon  is  observed  with  large  mach  numbers  of 
the  incident  flow.  The  greater  the  mach  number,  up  to  smaller 
distance  elevon  is  forward  from  spread  its  effect  on  the  side,  turned 
tc  flow  (Fig.  8).  During  transition  frcm  the  subsonic  to  supersonic 
speeds,  changes  the  form  cf  the  diagraa/cur ve  cf  pressure  on  elevon 
itself.  If  at  subsonic  speeds  the  form  of  d iagr am/curve  is  close  to 
triangular,  then  at  supersonic  speeds  it  is  close  to  rectangular. 

The  enumerated  above  special  feature/peculiarities  of  a change 
in  the  character  of  flow  arcund  of  the  wing  with  the  deflected 
elevons  explain  the  reasons  for  an  incidence/drop  in  the 
effectiveness  of  elevons  and  an  increase  in  their  hinge  moments 
during  transition  from  the  subsonic  to  supersonic  speeds. 
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A.  P.  Bazzhin,  0.  N.  Trusova,  and 
I.  P.  Chelysheva 


The  calculation  results  of  a flow  around 
a family  of  elliptic  cones  by  a flow  of  Ideal 
gas  at  large  angles  of  attack  were  presented 
in  works  [1,  2].  Subsequently,  several  variants 
of  flow  were  calculated  taking  Into  account  the 
real  properties  of  air,  which  are  In  a state  of 
thermodynamic  equilibrium.  These  calculation 
results  permit  one  to  evaluate  the  effect  of 
real  gas  properties,  which  proves  to  be  In- 
significant for  the  variants  of  flow,  examined 
in  works  [1]  and  [2j.  The  first  part  of  this 
paper  Is  devoted  to  this  problem. 

Calculated  aerodynamic  characteristics  of 
elliptic  cones  over  the  angle  of  attack  range 
from  30°  to  50°  in  the  case  of  an  Ideal  gas  are 
presented  In  the  second  part  of  this  paper.  The 
comparison  of  these  results  with  the  calculated 
and  experimental  results  of  other  authors  [3,  *0 
has  confirmed  the  validity  of  the  results  ob- 
tained ty  the  calculation  method  with  large 
angles  of  attack. 
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the  effect  of  real  properties  of  aih 
on  parameters  of  flow  near  an  elliptic 
cone 


A calculation  was  made  of  the  flow  around  an  elliptic  cone 
having  a cross-section  axes  ratio  of  6 « 2,  half-angle  aperture 
of  the  cone  in  a horizontal  plane  6h  = 15?,  with  the  angle  of 
attack  a » 30°.  The  air  was  examined  as  a three-component  gas 
consisting  of  78.08?  nitrogen,  20.95?  oxygen,  and  0.9755  argon, 
and  its  thermodynamic  functions  were  calculated  according  to 
the  standardized  program. 

The  incident  flew  velocities  were  equal  to  2350,  3356, 
and  6713  m/s,  which  corresponded  to  the  numbers  equaling  7,  10, 
and  20  (with  the  speed  of  sound  aa  = 335.6  m/s).  The  main  bulk 
of  calculation  data  in  works  [1]  and  [2]  was  obtained  at  M^  * 7. 

Figure  1 shows  the  position  of  shock  waves  near  the  cone  in 
the  perfect  and  imperfect  gases  at  velocities  ■ 2350  and 
6713  m/s.  The  difference  in  distance  from  the  body  to  the  shock 
wave  in  the  symmetry  plane  of  the  flow,  in  the  case  where  « 7 
(V„  * 2350  m/s),  comprises  about  10?.  The  absolute  shock  wave 
displacement  arising  when  considering  the  real  properties  of  air 
has  a negligible  change  through  out  their  duration.  The  same 
thing  applies  also  to  the  case  of  the  flow  with  * 6713  m/s. 
Transition  lines  II  near  the  lower  surface  change  together  with 
the  change  in  the  position  of  shock  waves;  however,  the  points 
of  transition  on  the  body  surface  are  displaced  very  little. 

Change  in  the  relative  distance  from  the  body  to  the  shock  wave 
e/t(  (c.j  - distance  from  the  body  to  the  wave  in  the  symmetry 
plane)  are  plotted  in  Fig.  2 as  a function  of  central  angle  u 
(see  Pig.  1).  The  range  uT  angles  w corresponds  to  the 

lower  surface  of  the  cone.  In  the  range  w < 80°  all  values  of 
quantity  e/Cg  fall  on  the  line  having  a width  of  not  more  than 
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0.0**.  In  other  words,  in  the  examined  area,  dependence  c/Sq  * 
can  be  represented  by  a curve  pertaining  to  the  perfect  gas,  with 
an  accuracy  to  within  '■*% . 


Functions  —=/,(“*)  and  —=/>(“)  have  a similar  nature 

Po 

(Fig.  3).  Vaiues  pQ  and  pQ  (pressure  and  density)  in  the  symmetry 
plane  of  the  flow  (on  the  wave  and  body)  are  referred  to  P„vmax 
and  pa,  respectively  (see  Table  1).  First  of  all  we  should  note 
the  extremely  slight  effect  of  the  real  properties  of  air  on  the 
magnitude  of  relative  pressure  when  - 2350  m/s  (M^  « 7).  The 
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Fig.  2. 


relative  pressure  on  the  body  surface  and  shock  wave,  when  0 < 

< oj  < 100°,  Is  virtually  Independent  of  the  real  properties  of 
air.  Changes  In  the  p^  value  of  the  body,  to  which  the  pressure 
on  the  body  Is  referred,  are  also  negligibly  small.  This  means  that 
the  values  of  aerodynamic  coefficients,  calculated  at  ■ 7 In 
the  case  of  a perfect  gas,  will  also  be  valid  with  high  accuracy 
In  the  case  of  an  Imperfect  gas  as  well. 

The  effect  of  the  real  properties  of  air  on  the  relative 
density  when  2350  m/s  (Mb  « 7)  Is  also  slight;  however,  the 
change  In  values  of  these  parameters  In  the  symmetry  plane,  with 
the  consideration  of  the  real  properties  of  air,  comprises  about 
8*  (see  Table  1 ) . 

Ch  .nge  In  the  relative  pressure  and  relative  density  on  the 
body  surface  remains  slight  when  considering  the  real  properties 
of  air,  even  at  velocity  * 6/13  m/s  (Mb  “ 20).  This  variation 
does  not  exceed  several  percent.  On  the  shock  wave,  especially 
on  the  upper  section,  the  change  In  the  relative  values  Is  more 
noticeable . 
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6713 

*0  pi'l  ■.  X 

h» 

5 0.918 

0,6*0 

0.630 

• 

*o  pert . 

g*s 

0,915 

0,835 

0.623 

If  we  assume  that  the  Indicated  nature  of  change  in  relative 
values  is  valid  not  only  for  the  examined  variants  of  flow,  but 
also  in  the  case  of  other  variants  close  to  those  examined,  then 
it  is  possible  to  propose  the  following  method  of  approximate 
calculation  of  the  effect  of  properties  of  an  imperfect  gas. 

Parameters  on  the  shock  wave  and  body  surface  in  the  symmetry 
plane,  when  calculating  the  real  properties  of  air,  vary  with  an 
accuracy  to  within  several  percent  (see  Table  1).  This  variation 
is  easily  obtained  when  the  slope  of  the  shock  wave  is  known. 

Then,  using  the  distribution  of  relative  parameters  obtained  for 
perfect  gas,  it  is  possible  to  obtain  the  real  distribution  of 
gas-dynamic  parameters  along  the  surface. 

With  regard  to  the  determination  of  the  shock  wave  inclination 
or  the  distance  from  the  shock  wave  to  the  body  in  the  flow  of 
an  Imperfect  gas,  as  a result  of  the  calculations  it  was  revealed 
that  the  ratio  of  distances  to  the  shock  wave  in  the  plane  symmetry 
in  the  case  of  an  Imperfect  gas,  to  the  corresponding  distance 
in  the  case  of  a perfect  gas,  is  equal,  with  high  accuracy,  to  the 
inverse  ratio  of  densities  on  the  shock  wave,  as  this  can  be  seen 
from  Table  2.  Consequently,  if  the  calculation  data  are  available 
for  a perfect  gas,  then  it  is  possible  to  approximately  determine 
value  eQ  « eQ  per.f  gas  on  the  shock  wave  in  the  symmetry  plane 
in  an  imperfect  gas;  then  to  find  value  Q0  perf.  gas,  refine  the 

PQ 

inclination  of  the  shock  wave  in  an  imperfect  gas,  and  find  a 
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■op*  precise  value  for  p0  on  the  shock  wavs  in  an  Imperfect  gas 
and  tha  new  value  Tor  cQ . Then,  using  the  available  dependence 

for  perfect  gas.  It  la  possible  to  determine  the  location 

of  the  shock  wave  near  the  lower  surface  of  the  cone  In  an  Imperfect 

gas. 

AERODYNAMIC  CHARACTERISTICS  OK  ELLIPTIC 
CONES  AT  LARUE  ANOLES  OK  ATTACK 


Aerodynamic  characteristics  of  elliptic  cones  were  calculated 
for  a perfect  gas  with  • 7 over  the  angle  or  attack  range  of 
30  to  50°.  The  error  In  determining  the  flow  parameters  on  the 
cone  surface.  In  particular  the  distribution  of  pressure  according 
to  the  carried  out  estimations,  comprises  the  value  on  the  orler 
of  IS.  The  error  In  calculating  the  aerodynamic  coefficients 
should  be  on  the  magnitude  of  the  same  order.  The  forces  acting 
on  that  part  of  the  upper  body  surface  where  the  flow  was  not 
calculated  were  not  considered,  when  calculating  the  forces  and 
moment.  However,  It  Is  entirely  obvious  that  If  the  streamlining 
occurs  without  a break  in  the  flow  over  the  upper  surface  of  the 
body  the  forces  ate  very  slight. 

Aerodynamic  coefficients  were  calculated  using  the  formulas:" 

normal  force  coefficient 


Cm 


* v ■»«  r . 


axial  force  coefficient 


(1) 


(2) 
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(t! 

coefficient  of  longitudinal  momentum 

e-“  + at>e»’  (3) 


where  ; S - cone  area  In  the  plan. 

Limits  of  Integration  correspond  to  the  bypass  of  the  cross- 
section  contour  from  the  symmetry  plane  to  the  last  point  at 
which  the  solution  Is  known,  i.e.,  the  change  in  the  variable  Is 
first  from  zero  to  b and  then  from  b to  n#.  In  range  (2)  t(b)  ■ 

■ 0.  In  the  vicinity  of  this  point  the  Integration  was  carried 
out  by  means  of  variable  for  which  the  following  substitution  was 

made  «*  ~ **  l*“ l(‘>  or  E — 1(0—  equations  of  transverse  elliptic 

section;  t and  n conical  variables]. 

The  coefficients  of  the  lift  and  resistance  forces  were 
determined  using  the  c^  and  cT>- 

c,  — Cn  cos  » — Cj  si. i a.  c , - c.v  sin  » -f  cj  cos  ». 

Figure  U shows  the  aerodynamic  coefficients  cx>  cy , cm  and 
aerodynamic  quality  K as  a function  of  the  angle  of  attack  of 
cones  with  opening  semiangles  6^  » 10°,  15°  and  20°  and  the  ratio 
of  axes  of  the  cross  section  6 » 1;  2 and  3-  Values  cx,  c^ , cm 
and  K ( 6 » 1)  are  given  over  the  entire  range  of  the  angles  of 
attack  from  zero  to  50°.  The  solid  lines  indicate  the  calculation 
data  from  work  [3]  at  small  angles  of  attack.  Solid  lines  In  the 
angle  of  attack  range  from  30°  to  S0°  indicate  the  results  obtained 
In  this  work.  The  experimental  data  from  work  [*0  are  plotted 
by  different  points  pertaining  to  air.  The  axes  indicate  the 
experimental  data  obtained  by  authors  earlier  at  ■ 6.  Such  a 


8 


comparison  of  the  various  data  had  one  purpose  In  mind  - to  arrive 
at  a concept  concerning  the  nature  of  change  of  the  aerodynamic 
characteristics  of  elliptic  cones  over  the  entire  range  of  the 
angles  of  attack  and  to  determine  the  validity  of  the  calculation 
data  obtained  by  us  over  the  angle  of  attack  of  attack  range  from 
30°  to  50°.  As  can  be  seen,  as  a whole,  there  Is  good  qualitative 
and  quantitative  agreement  between  all  the  results  presented.  The 
dashed  lines  In  the  Intermediate  angle  of  attack  range  can  be 
considered  as  a possible  Interpolation  of  the  aerodynamic  coefficient 
values  In  this  area.  The  remaining  curves  In  Fig.  ^ represent  the 
aerodynamic  characteristics  of  elliptic  cones  In  the  range  of 
large  angles  of  attack  at  different  values  of  6 . 


/W_-  7 


The  effect  of  6 with  dirrerent  constant  parameters  Is  shown 
In  Fig.  5.  We  will  note  that  the  aerodynamic  quality  of  cones 
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always  Increases  when  passing  to  greater  elliptlcity  of  the  cone 
cross  section.  The  coefficient  of  longitudinal  momentum  over  the 
range  of  large  angles  of  attack  changes  almost  linearly  with  the 
angle  of  attack.  According  to  formula  (3),  the  position  of  the 

pressure  center  of  an  elliptic  cone  Is  determined  by  value  (1  + 

2 

+ a ),  i.e.,  only  by  the  opening  semiangle  of  the  cone  In  the 
symmetry  plane  of  the  flow. 
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ST  DDT  CP  THE  PLOW  OF  GAS  IN  A CYLINDRICAL  CHANNEL  DORING  THE  SUDDEN 
EXPANSION  OP  SONIC  FLO*. 

G.  P.  Glotov,  E.  K.  Moroz. 

Is  carried  out  the  study  of  flow  during  the  sudden  expansion  of  \ 

scnic  airflow  in  axisyaaetric  cylindrical  channel,  changed  relative 
length  and  the  area  of  chanrel  in  ranges  / * 1. 5-4.5;  P = 1. 5-3.0. 

Are  investigated  the  special  feature/peculiarities  of  flow  in  the 
area  of  the  connection  cf  flow  to  the  wall  cf  channel  and  is 
establish/installed  the  existence  of  the  single  condition  of 
connection . 

The  problea  of  the  connection  of  turbulent  supersonic  flow  - one 
of  the  basic  with  solution  of  which  we  encounter  in  a series  of  the 
gas-dynaaic  eguipaent/devices:  air  intakes,  ejector  nozzles,  the 
caaera/chaaber  of  Eiffel,  etc.  One  Of  the  piobleas  in  this  case 
consists  of  the  deterai nation  of  pressure  cf  stagnation  zone.  The 
experiaental  investigation  cf  pressure  in  stagnation  zone  at  the 
large  lengths  of  axisyaaetric  cylindrical  channel  (in  connection  with 
ejectcrs  with  the  zero  coefficient  of  ejection)  was  for  the  first 
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t time  carried  out  by  G.  L.  Grodzovskiy  et  al.  into  1953  [1]. 

Subsequently  similar  data  were  obtained  in  a series  of  the  works  (for 
example,  see  [2]  and  [3]).  The  effect  of  the  length  of  channel  on 
pressure  in  stagnation  zone  was  for  the  first  time 

establish/installed  in  the  experiments  cf  G.  L.  Grodzovskiy  and  V.  T. 

Zhdanov  whose  results  were  presented  in  work  [3]. 

The  beginning  of  theoretical  studies  cf  the  problem  of  the  i 

connection  of  the  turbulent  flow  was  placed  in  work  [4]  and  it.  is 
continued  in  [5],  [6]. 

In  this  article  are  investigated  the  tasic  physical  phenomena, 
which  appear  during  the  connection  cf  turbulent  supersonic  flow  to 
wall,  and  the  condition  of  the  connecticn  ci  separating  flow  line. 

For  this  purpose,  was  carried  out  the  experimental  study  of  the 
flow  of  turbulent  supersonic  flow  in  cylindrical  channel  with  sudden 
expansion.  The  schematic  of  the  model  of  channel  with  designations 
and  the  geometric  parameters  of  the  investigated  versions  are  given 
to  Fig.  1. 

In  experiments  discretely  changed  the  relative  length  of 


channel.  The  range  of  the  lengths  of  channel,  in  reference  to  the 
height/altitude  of  step,  the  equal  to  the  half-difference  of 
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w diameters  earner a/chaibers  and  nozzle  throats 

i 1.5-4. 5,  but  the  area  of  channel,  in  re 
critical  nozzle,  changed  in  the  range  F = 1.5 

Page  54. 

Boundary  layer  thickness  in  nozzle  throat,  re 
height/altitude  of  step,  it  was  equal  to  6 = 
number,  calculated  according  to  critical  thro 
Be  = (2-7.2)  x 10*.  Testings  were  conducted  a 
in  precoab ustion  chamber,  equal  to  3-8  atm  (at 
1„  = 290°K. 

During  testings,  besides  the  total  pressure  in  precombustion 
chamber,  were  measured  the  pressure  in  stagnaticn  zone  p,  and  the 
distribution  of  pressure  according  to  the  wall  cf  channel  p with  the 
aid  of  static- pressure  probes  0.8  is  in  diameter,  arrange/located 
with  space  1.5  am.  The  accuracy/precisicn  of  the  determination  of 
relative  pressure  was  *1o/o. 

The  picture  of  flow  at  output/yield  frem  channel  was 
photographed  by  Toepler's  instrument.  With  the  aid  of  oil  film 
(aixture  of  oil  and  carton  black)  was  visualized  the  picture  of  flow 
cn  the  wall  of  channel  and  in  the  aeridian  plane  of  stagnation  zone. 


h ~ ~ , comprised 

ference  to  the  area  of 


f erred  to  the 
tc  0.10-0.25.  Reynolds 
at  diameter,  comprised 
t the  pressure  air  flow 
s.),  and  to  temperature 
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Fig.  2,  gives  the  typical  dependence  of  relative  pressure  in 
stagnation  zone  P!  (pt  = Pi/p0»  where  p0  - averaged  according  to 
expenditure/consumption  total  pressure  flow  in  the  section/shear  of 
sonic  nozzle)  on  relative  nozzle  pressure  p0  (pc  = Po/P2»  where  p2  - 
aabient  pressure).  Are  isolated  three  characteristic 

conditions/aodes:  1 - ccndi tions/aode  of  the  connected  flow,  which  is 
characterized  by  constant  quantity  of  relative  pressure  in  stagnation 
zone,  2 - transient  conditicns/mode,  3 - separating  cond it ions/mode. 
Further  analysis  of  the  obtained  experimental  data  is  conducted  for 
conditions/aode  1. 
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Fig.  2. 


0 2 t l 


Pig.  3. 


Key:  (1)  foreward  stroke.  (2)  back  strcke. 


Page  55. 

The  generalized  dependences  /' i /(/)  at  the  different  values  of 
the  relative  areas  of  channels  F are  given  to  Fig.  3.  As  can  be  seen 
frca  this  curve/graph,  the  value  of  relative  pressure  in  stagnation 
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zone  pj  remains  co nstan t/ in var iab le  fee  this  value  of  F during  the 
decrease  of  the  relative  length  of  channel  / to  certain  value  which 
let  us  designate  Me  will  call  the  area  cf  the  lengths  of  channel 

l i the  region  of  self-siiilar  flow  alccg  the  length.  During 
further  decrease  of  the  relative  length  of  channel  to  / /.,  value 

Pi  begins  to  decrease  (region  of  ncn-self-simulating  flow)  . Pig.  3, 
gives  the  boundary  of  the  region  of  self-siiilar  flow. 

Comparison  showed  that  our  data  on  pressure  in  stagnation  zone 
for  the  region  of  self-similar  flow  will  agree  well  with  other 
authors's  data,  obtained  during  the  discharge  of  sonic  flow  into 
cylindrical  channel. 

Mas  also  carried  out  the  estimation  of  the  known  criteria  of  the 
connection  of  turbulent  flow. 

One  of  the  most  successful  criteria  is  the  examined  in  works 
[5],  [6]  condition  for  the  angle  of  the  correction  of  flow  V 
where  H0  - mach  number  cn  the  boundary  of  inviscid  jet.  At  angle 
is  understood  the  angle  of  incidence  with  the  wall  of  the  channel  of 
the  boundary  of  the  inviscid  jet,  constructed  by  method  of 
characteristics  according  to  the  measured  in  experiment  sense  of 
pressures  p0/Pi  (see  schematic  in  Fig.  H)  . As  is  shown  comparison, 
the  values  «r  calculated  according  to  the  results  of  this  work. 
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in  the  region  of  self-siailar  flow  virtually  coincide  with  the  values 
1'.  given  for  the  appropriate  flow  in  work  [5]  (Fig.  4).  However,  in 
the  region  of  non-self-simulating  flew,  angles  <1  grow/rise  and  for 
each  value  of  F is  obtained  its  dependence  t'  >MM„>  Thus,  in  the 
exaained  case  the  criterion  of  the  connection  of  turbulent  flow  in 
the  fora  of  single  dependence  <i  t(M„)  is  valid  only  in  the  range  of 
self-siaulating  flow  and  is  not  spread  to  tfce  non-self -simulating 
region  (this  observation  is  related  also  to  the  correlation  parameter 
x.  introduced  in  work  [5]).  Therefore  is  necessary  the  search  of 
other  aore  coanon/gener al/tctal  criteria  of  connection. 

For  purpose  of  the  explanation  of  the  physical  picture  of  flow 
in  the  region  of  the  connection  of  flew  to  the  wall  of  cylindrical 
channel,  was  carried  out  the  visualization  of  tfce  picture  of  flow  on 
the  wall  of  channel  and  in  the  zone  of  nixing  with  the  simultaneous 
measurement  of  static  pressure  distribution  on  wall  and  the 
photographing  of  flow  at  output/yield  with  the  aid  of  Toepler's 
instrument. 

The  examination  of  the  obtained  photographs  of  oil  film  and 
their  comparison  with  tfce  diagram/curves  of  the  distribution  of 
pressure  on  wall  nade  it  possible  to  present  the  real  picture  of  flow 
in  the  region  of  connection  (Fig.  5a).  In  photographs  are  visible 
three  zones.  In  zone  I (we  examine  fren  nozzle  edge)  oil  remained 
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intact.  The  coaparison  of  this  zone  with  the  air-load  distribution  on 
wall  which  in  this  zone  it  is  constant,  shews  that  the  flow  here 
either  entirely  is  absent  cr  so  weak  that  it  dees  not  act  on  oil 
fill.  In  zone  II,  are  observed  the  longitudinal  overflows  of  oil  and 
an  insignificant  change  in  the  pressure,  which  indicates  the  presence 
cf  weak  current  in  the  liaited  region  towards  nozzle.  In  zone  III  oil 
is  washed  off  coapletely  cn  an  entire  wall,  except  the  narrow 
transverse  band  with  a width  approxinately  1 ci  (line  of  the 
connection  P) . Flow  in  this  zone  is  acccapanied  by  sharp  pressure 
increase  on  wall  tc  certain  aaxiaut  value  ,/p, 
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This  testifies  to  the  presence  of  powerful  flow  with  interface  (line 
II),  to  the  left  of  which  flow  is  directed  toward  nozzle,  and  to  the 
right  * toward  the  section/shear  of  channel.  It  is  logical  to  assume 
that  the  interface  represents  by  itself  the  feint  of  rendezvous  of 
separating  flow  line  (line  I)  with  the  wall  of  channel. 

Pig.  6,  gives  the  photograph  of  oil  fill  to  longitudinal  plate 
and  the  corresponding  schematic  of  flow.  On  figure  are  noted:  1 - 
separating  flow  line,  2 - boundary  of  zero  longitudinal  velocities,  3 
- boundary  of  flow,  4 - boundary  of  inviscid  jet,  5 - point  of  the 
connection,  by  6 - duct/contour  of  the  plate.  In  photograph  clearly 
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is  outlined  the  rotation  of  flow  line  in  viscous  layer, 
arrange/located  lower  than  separating  line,  into  stagnation  zone  and 
the  foraation/education  of  the  reverse/inverse  flow  near  the  wall. 

The  lass  of  the  reflux  gas  returns  to  the  aein  flow,  foraing  the 
local  eddy/vortex  between  wall  and  boundary  of  flow.  The  longitudinal 
size/diaension  of  this  eddy/vortex  virtually  coincides  with  the 
extent  of  zone  II  (Fig.  6). 

The  part  of  the  viscous  layer,  arrange/located  higher  than 
separating  flow  line  after  aeeting  with  wall  turns  to  output/yield 
frca  channel.  During  this  rctation  in  flow,  appears  the  system  of 
characteristics.  Intersecting,  they  create  the  oblique  shock  wave, 
seen  at  output/yield  free  tie  channel  (see  Fig.  5). 

Is  of  interest  the  ccaparison  of  calculated  and  deterained  in 
experiaents  in  the  positions  of  separating  flow  line.  In  work  [5]  as 
separating  line  is  accepted  the  boundary  of  inviscid  jet.  For  the 
region  of  self-siailar  flow,  it  is  possible  to  note  the  satisfactory 
ccnforaity  of  the  calculated  and  experiaental  results  (see  Fig.  6). 


DOC  = 78068003 


PAGE 


Only  of  wall  Itself  begins  their  disagreement,  as  a result  of  which 
the  calculated  boundary  of  inviscid  jet  neets  the  wall  of  channel  at 
the  point,  which  lies  approxinately  tc  7-  lOc/c  further  from  nozzle. 


than  the  real  line  of  connection 


In  the  case  of  snail  length  of  channel  (ncn-self-sinu lating  zone 
of  flow)  the  calculated  boundary  of  the  inviscid  jet  A (Fig.  5b)  can 
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even  exceed  the  Units  of  channel,  although  in  actuality  the  flow  is 
that  connected.  For  such  channels,  as  shows  the  visualization  of  flow 
in  aeridian  plane,  separating  flow  line  B considerably  differs  fron 
line  of  denarcation  of  inviscid  jet  and  aore  steeply  it  turns  to  wall 
(Fig.  5b). 

The  specific  in  experiments  positions  of  the  line  of  connection 
on  the  wall  of  channel  aade  it  passible  by  the  aeasured 
diagraa/curves  of  pressure  to  deteraine  the  pressure  at  attachaent 
point,  equal  to  the  total  pressure  on  separating  flow  line  (Fig.  5). 
The  coaparison  of  this  value  of  pressure  with  pressure  in  stagnation 
zone  shows  that  their  sense  for  all  investigated  values  of  relative 
areas  and  lengths  of  channels  apprexiaately  is  constant  and  is  equal 
tc  Pulp,  s to  1.9  j^O.05  (Fig.  7a). 

The  results  of  processing  given  works  [6}  - [9]  show  (Fig. 

7b) that  during  the  flow  around  flat/plane  step  is  observed  certain 
tendency  toward  an  increase  in  value  P«IP\  Ipa/pi  * 1.7-2  with  d0  = 

2.  1-4. 4)  . In  the  first  apprexiaatien,  this  sense  can  be  accepted 
equal  to  1.9.  k change  of  the  relative  aaxiaua  pressure  on  wall 
depending  on  nuaber  n0  both  for  the  flow  in  channel  and  during  the 
flow  around  flat/plane  step  does  not  in  practice  affect  value  pn/p, 
(fig.  7). 
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Fig.  6. 

Key:  (1).  Region  of  self-siailar  flow.  (2).  Zone. 

Page  58. 

The  relationship/ratio  pu/p,  1.9,  obtained  fee  flow  in  a xif.y  ametric 
cylindrical  channel  and  during  the  flow  around  flat/plane  step,  can 
used  in  the  analysis  of  the  connection  of  supersonic  turbulent 
flew  on  wall  in  the  range  of  nuabers  N0  ;s  2-3.5. 

The  conducted  investigations  aade  it  pcssible  to  also  explain 
the  aechanisa  of  the  effect  of  the  length  of  channel  on  relative 


\ 
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pressure  in  the  stagnation  zone  (see  Fig.  5).  In  accordance  with  a 
change  of  the  basic  flow  parameters  in  the  region  of  connection,  it 

is  possible  to  isolate  three  reference  lengths  cf  the  channel: 

I ^ l < /„ 

^ and  I lK i,  (this  tc  some  degree  of  analogous  the 

introduction  critical  points  for  the  case  of  the  flow  around 
flat/plane  step  [5]).  The  specific  above  relative  length  of  channel 
/.,  (see  Fig.  3)  it  corresponds  to  xaxisum  pressure  increase  on  wall 
in  the  region  of  self-sisilar  flow.  During  the  decrease  of  the  length  \ 

cf  channel  to  /a  (position  of  the  section/shear  of  channel  1 and  2, 

Fig.  5a)  relative  pressure  in  stagnation  zone  remains  constant.  Tn  j 

this  case,  remain  without  change  the  position  cf  the  line  of 

connection,  maximum  pressure  on  wall  and  the  angle  of  the  slope  of 

the  resulting  shock  wave,  observed  in  output/yield  from  the  channel 

(see  Fig.  5a),  it  which  indicates  the  invariability  of  flow 

disturbance  in  local  region  after  the  line  of  the  connection  (the 

value  cf  slope  angle  very  weakly  depending  on  value  of  F) . The 

observing  when  / decreases  of  pressure  on  the  wall  (see  Fig. 

5a)  testify  to  the  presence  on  this  section  of  the  accelerated  j 

supersonic  flow. 

iith  the  decrease  of  the  length  of  channel  into  the  region  of 
non-self-simulating  flew  / /„  tc  certain  value  which  let  us 

call/naae  critical  P,  the  fora  of  the  distribution  curve  of  pressure 
and  the  position  of  the  line  of  connection  remain  without  change.  In 
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this  case,  maximum  pressure  on  the  vail  of  channel  and  the  angle  of 
the  slope  of  the  resulting  jump  at  output/yield  decrease.  The 
decrease  of  flow  disturbance  after  the  line  of  connection  leads  to 
the  decrease  of  relative  pressure  in  stagnation  zone.  In  this  case, 

hiIiui  by  pressure  on  wall  is  greater  than  pressure  environment 

ip I).. 


During  the  decrease  of  the  le 
4,  Pig.  5b)  occurs  the  shift/shear 
line  of  connection  to  nozzle,  as  a 
stagnation  zone  it  decreases. 


ngth  of 
of  the 
result 


c tannel 
air-load 
cf  which 


/ (position 
distribution  and 
the  length  of 
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a)  • ?,0  S*  - <t*  0,055+0,1* 

■»  Z 5 * ft  - t?,S mm  Hf]; 

• 3,0  -O'  0,005\  0,9  ,ft -</S;25  mm 

Fig.  7. 

Key:  (1).  Axisynaetric  flow. 

Page  59. 

On  section  froa  attachaent  point  tc  the  point,  which  corresponds  to 
critical  length  (In  I I, ,.) , flow  in  wall  viscous  layer  subsonic  [9], 
Therefore  when  / /,,  the  picture  of  flow  either  aust  be  broken  as  a 

result  of  the  report/coaaun ication  of  stagnaticn  zone  with 
environaent  or,  at  the  sufficiently  large  piessure  flow,  the  line  of 
connection  aust  aove  froa  the  section/shear  of  channel,  that  also  is 
observed  in  experiaent.  In  this  case,  aaxiaum  pressure  on  wall 
continues  to  decrease,  reaaining  acre  than  pressure  environaent,  and 
with  respiect  it  decreases  pressure  in  stagnaticn  zone. 
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Thus,  for  the  investigated  co edit ions/iodes  of  the  connected 
flew  pressure  environment  does  not  affect  value  pt . Pressure  in 
stagnation  zone  depends  cn  Baximum  pressure  on  the  wall  of  channel 
after  the  line  of  connection,  that  it  is  necessary  to  consider  during 
the  development  of  the  calculation  method. 
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flow  of  gas  in  a flat  duct,  caused  longitudinal  gradient  of  the 
temperature  at  Knudsen's  arbitrary  number. 

H.  N.  Kogan,  N.  K.  flakashev. 

For  the  kinetic  model  equation  of  Boltzmann  in  linear  setting, 
is  solved  the  problem  of  tenperature  creep  in  flat  duct  for  the 
arbitrary  values  of  Knudsen's  number. 

Obtained  approximate  analytical  scluticn.  On  the  basis  of  this 
solution,  are  made  evaluations  of  possible  faults  of  measurement  of 
pressure,  for  example,  the  heated  gas  v it h the  aid  of  "cold" 
instrument . 

As  she  was  noted  already  by  Harwell  [ 1 ],  if  along  vail  is  a 
gradient  of  temperature,  then  the  ccming  into  contact  with  it  gas 
moves  relative  to  mall.  This  motion  calls  thermal  slip  or  creep.  The 
gas  flow  in  this  case  depends  on  tbe  number  of  Knudsen  K n,  equal  to 
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the  ratio  of  lean  free  path  X to  the  width  cf  channel  or  the  diameter 
of  tube  d.  If  by  the  unevenly  heated  tube  are  connected  two 
containers  with  different  temperature,  then  the  equilibrium  (zero 
expenditure/consumption  through  the  tube)  steps  at  certain  the 
pressure  differential  which  also  depends  on  Knudsen's  number. 

These  phenomena  can  exert  the  essential  influence,  for  example, 
during  low-pressure  measurement  the  heated  gas  by  "cold"  instrument; 
in  the  porous  media  they  can  cause  flow  or  the  pressure 
differentials. 

Vsfiou-ii  edit*  cl  the  F/cu.  cl  c?  cj<3c,  e ¥-prcb5ec(  by  ie i»pe rcliarf. 

qi'  Hct  / f i‘~t  cl  till  u ' c5//*  . O " 2.  2 Z&r-  / e gd  i/7  £'2  7-  £ 7 J7  * 

Different  cases  of  the  flow  of  gas,  caused  in  flat/plane  duct  by 

the  gradient  of  the  temperature  at  the  arbitrary  values  of  Knudsen's 
nuxfcer.  Approximate  soluticr  of  the  model  equation  of  Boltzman  is 
obtained  in  analytical  for*. 

1.  Let  us  examine  gas  between  two  infinite  parallel  motionless 

plates.  The  temperature  cf  walls  Tw  is  changed  on  z.  Let  us  consider 

that  this  change  is  small,  so  that 

Tw  7„  1 1 t(z)l;  t (z)'  1;  T„  Tu  :<>),  (11) 


and  problem  is  linearized.  Let  us  assume  also  that  the  walls  reflect 

molecules  according  to  Maxwellian  law;  the  temperature  of  the 

molecules  reflected  is  equal  to  the  temperature  of  wall 1 * * * * *  7 V 
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For  the  solution  of  prcblen,  we  will  use  the  model  equation  of 
Bcitzaann  (for  example,  see  [8]): 

■.<.£-*«</.  /».  /. 

where  f - distribution  function;  n,m,i  (itf  \t)  - respectively 

numerical  density,  mass  and  the  velocity  of  nolecules;  u and  T - 
■acroscopic  velocity  and  temperature;  k - Ecltzmann  constant. 


Condition  on  the  wall: 


/(•  "w  ( '2r.k  Tw  J W'C-m  1-  ^ 

where  mu,  it  is  determined  from  the  condition  cf  nonpassage. 


Solution  let  us  search  for  in  the  form 


/ 1 1 + *(*.  «)l; 

i " Y „n  L 


/ m \‘J  I m t 
/,in  Ink  7,,  J CXP  f IkTa 


where  n0  = n(0,  0)  ; * - small  addition  squares  of  which  we  disregard. 
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Linearizing  equation  (1.2)  and  after  aaking  it  dimensionless,  we 


will  obtain 


v , <>'- 
a dx, 


V . <>? 


a <>Z , 


V I 2v 


I “ Tj  j 'l-1 


-?•  (1-5) 


Lr«  here  introduced  designations: 


j i j ; A«  ~^kT  ’ ? 2 ^ - v7 

n w„(l -|  v);  T — 7',, (I  ! *);  a - A nud  V h„  i Kn 

ii  , u,  l /iu  - ' f 'f  X',  dv\  v - - ft-  -u  ; 

r.  - n 2 

- 1 f 

T ..  — j — ^ e v ■ yd  v — v. 


Eoundary  condition  (1.3)  takes  the  fern: 


K.(  1 t ■ 


2.  Let  temperature  of  wall  change  linearly: 


Page  71. 


Then  the  solution  of  equation  (1.6)  it  is  pcssitle  to  search  for  in 
the  form 


y 1 12  V-  V), 


in  this  case 
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p hnttTJ\  2 azj. 


(2  U 


and  for  •}>  are  obtained  following  equation  and  the  boundary 
ccndi tion: 

v,  , <id  , . .. 

~i  H-l  + a ^ » 1 J"': 

•>(;*  . «,  •>')  <>• 

Set/assuaing  temporarily  u,  by  known  function  and  integrating 
(2.4),  we  obtain  integral  equation  for  V 


(2-4) 


•M*i.  vx  °)  a 


2m,  (s) 


(/</  v‘‘ 


<■,  <Js 

V, 


(2.5) 


Multiplying  this  equation  on  t^exp  ( - V1)  and  integrating  by 

> 

velocities  v,  we  will  obtain  integral  equation  for  ut(x1): 


«i 


V' 


ad  . 

"!<•')  a ^ 


(«  t,  - s ,)  rf*  - 7 - J A (»  v i ' ) 


V 


•M»)  \ • i'!  )•  ,,  1 1 

A,<*>  - ax 


(2.0) 


Function  ./ _ i j jc , s|)  has  logarithaic  special 

feature/peculiarity  with  s * x,/  Therefore  approximately  it  is 
possible  to  assume  u, (s)  * u(  (xt)  and  to  reaove  u,  ( x t ) from  under 
integral.  In  this  approximation  of  solution  of  equation  (2.6)  is 
obtained  in  an  explicit  fen: 
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ml 

a 


(A„  V r.)  f ^ ( A, 

A„ 


1 


x,']  t Jn 


,(  1 


y Jn  ^ *|J 


(2  7) 


Analogous  approach/approximation  was  used  in  [ 8 ] for  the 
examination  of  Poiseuille  flow  in  the  sale  setting.  Problem  in  [8] 
differs  froa  that  exaaine/ccnsidered  ccly  by  the  fact  that  in  it  '«• 
it  is  set/assuaed  by  constant  and  is  assign/p rescribed  pressure 
gradient.  In  this  case 


H 1 p 


b 


bd  A„  | r. 

‘2*  A, 

1 dp 
Pn  dz 


C2.H) 
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The  coaparison  of  the  obtained  solution  with  precise  solution  of 
Cherchin' yani  [9]  showed  its  satisfactory  accur acy /precision  (Fig. 

1) . It  is  possible  to  expect  the  saae  accuracy/precision  of 
approach/approxiaation,  also,  in  the  problem  in  question.  Accordingly 
(2.1)  and  (2.3)  solution  (2.7)  corresponds  to  the  gradient  of 
temperature  a * (l/To^dT/dz)  and  to  pressure  gradient  b = 5/2  a.  Since 
in  linear  setting  it  is  valid  superposition,  let  us  exclude  with  the 
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aid  of  (2.8)  from  solution  (2.7)  the  part,  caused  by  pressure 
gradient.  Then  with  dr/dz  = 0 we  obtain 


a i / 


nil  I Aj 


(2  9) 


3.  Knowing  distribution  of  velocities,  it  is  easy  to  determine 
volumetric  flow  rate: 


I 

1 

Qi  ,,  d \ ul7  p(x,)  dx,\ 
*1 


('ll) 


Results  of  its  numerical  calculation  are  represented  in  Pig.  1. 
Theme  results  were  obtained  with  the  aid  of  the  tables  of  integrals 
J„.  given  in  [10].  At  Pciseuille  flow,  as  is  known,  has  the  minimum 
of  ezpenditure/consumption  (Knudsen's  paradox).  As  can  be  seen  from 
Fig.  1 of  the  flow,  caused  ty  temperature  gradient, 

expenditure/consumption  is  changed  aonotonically,  after  grow/r ising 
with  the  decrease  of  pressure.  Zero  expenditure/consumption  is 
established  with  a specific  ratio  between  the  gradient  of  temperature 
and  pressure  gradient.  Set/assuning  <Ji  , accordingly  (3.1)  we  have 


V 

P 


K («) 


A T 
T 


K(  a) 


i/i  («) 


CU) 


To  dependence  K(a)  it  is  represented  in  Fig.  2.  During  the 
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measurement  of  temperature  on  the  value  of  the  order  of  magnitude  of 
temperature  itself  the  pressure  differential  can  coaprise  to  59o/o  of 
its  average  value.  This  aaxinun  value  K (a)  is  reached  in  free 
aolecular  condi tio ns/aode  when  Kn  - and  it  was  obtained  from  the 
solution  of  the  equation  of  Boltzaann  for  this  case.  Thus,  for 
instance,  if  instruaent  and  the  aeasured  volume  are  connected  by  the 
tube  with  a diaaeter  of  1 at  and  is  aeasured  pressure  order  0.  1 mm  Hg 
with  AT/T  = 0.3,  then  the  error  in  readings  will  comprise 
approximately  8o/o. 


0 t 2 J t*  S 6 u 


Fig.  1. 

Key:  ( 1)  . work. 

Page  73. 

Fee  a comparison  Fig.  1 and  2 by  dotted  line  show  results  for 
Na vier-Stokes  eguation  with  slip  conditions  on  boundary. 

With  large  a the  accuracy /precision  of  obtained  approximate 
solution  falls  and  it  gives  inaccurate  asymptotic  behavior  when  a 9 
i.e.,  when  Kn  9 0.  With  Knudsen's  small  numbers,  it  is  possible  to 
utilize  Sa  vier-Stokes  eguations  with  the  conditions  of  temperature 
slip  on  wall. 
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According  to  the  solution,  obtained  in  [8],  this  slip 

1 is  equal 

ii  . p riti 

since  for  here  the  aodel  equation  in  question  the  coefficient  of 
ductility/toughness/visccsity  p = kT/A. 

PCOTHOTE  * . in  work  [7]  is  cbtained  the  value  cf  coefficient,  equal 
to  0.383.  ENDPOOTNOTE. 


Consequently,  the  rates  cf  flow,  caused  by  creep,  are  of  the  order 
pa,  and  the  inertia  and  viscous  teras  of  Na vier-Stokes  equations  - an 
erder  p2a2,  i.e.,  the  saae  erder  as  scae  of  the  additive  terns, 
entering  the  equations  of  Barnet.  However,  it  is  possible  to  show 
that  for  here  a saall  linear  gradient  of  teiperatures  <t  r and  a << 

1)  in  question  the  solutions  of  the  equations  of  Barnett  and 
Havier-Stokes  coincide  (with  p = ccnst,  i.e.,  with  b = 0) : 

„,(v.  I)  const  O.v"’!  : (,n 

Q,  0.12  ar-  n»  <*r» 

1 ft  »*  < 


Bork  [7]  shows,  that  for  the  axisyanetric  case  solution  in  the 
approach/approx iaations  cf  Ha vier-Stckes  and  Barnett  they  do  not 
coincide  and  in  the  appropriate  expansion/decoaposition  of  type  (3.5) 
eaters  tera  0 ( a'2).  Thus,  in  here  the  flat/plane  case  in  question 
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Ha vier-Stokes* s asyaptotic  behavior  (3.5)  aust  te  satisfactory 
already  with  not  very  large  a (see  Fig.  1)  . 

Let  us  note  that  Navier-Stokes *s  asyaptotic  behavior  for 
Poiseuille  flow  is  such  worse  (see  Fig.  1)  . It  is  real/actual, 
accepting  for  the  rate  cf  slip  (see  [8]) 

I llu, 

“>»  — * ."*2  ---  ,fjt|  u . 

we  have 


< i iii 
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Fig.  2. 

f/rj/Z  Zfj 


Fig.  3. 

Page  74. 

Bith  zero  flow  rate  the  gas  of  wall  flews  to  one  side,  and  in 


center  another  (Fig.  3) 
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It  is  interesting  to  note  that  under  tie  condition  of  the  zero 
flew  rate  through  the  section  and  at  the  sufficiently  high  values  of 
number  Kn  the  velocity  profile  is  such,  that  afccut  wall  gas  flows  in 
the  direction,  opposite  to  the  gradient  of  the  temperature  of  the 
walls  (see  Fig.  3)  . At  the  sane  tine  from  sclution  this  same  of 
problem  with  the  small  numbers  Kn,  obtained  from  Navier-Stokes 
equation,  it  follows  that  the  gas  velocity  cf  wall  has  another 
direction,  i.e.,  during  a change  ir.  the  number  of  Knudsen,  under  the 
condition  of  the  zero  flow  rate  through  the  section,  gas  velocity  of 
wall  reverses  the  sign. 
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OPTIMIZATION  OF  THE  FLYING  EANGE  OF  VEHICLE  IN  THE  ATMOSPHERE  TAKING 
INTO  ACCOUNT  LIMITATION  TO  COMPLETI  CVEFLOAI. 

V.  V.  Dikusar,  A.  A.  Shilov. 

Is  examined  the  problem  of  the  determination  of  maneuverability 
capabilities  of  the  space  vehicle,  which  possesses  lift,  during 
reduction  in  the  atmosphere  talcing  into  account  limitation  to  phase 
coordinates.  Problem  is  solved  with  the  use  of  classical  principle  of 
L.  S.  Pontriagin's  maximum.  Are  given  numerical  examples. 

Great  practical  interest  represents  the  application/use  of 
methods  of  the  optimization  of  trajectories  in  the  presence  of 
limitations  during  the  function  of  phase  coordinates.  To  theoretical 
guestions  of  these  problems  are  dedicated  works  [ 1 ] - [ 5 ].  Tn  works 
[ 1 ],  f2],  [4],  [5]  the  principle  of  maximum  is  demonstrated  for  the 
case  «uen  in  the  optimum  trajectory  in  question  everywhere  is 
retained  the  local  effectiveness  of  control.  This  case  is  called 
regula  r [ 3 ]. 


The  basic  difficulty  of  applying  the  principle  of  maximum  is 
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connected  with  the  need  for  the  solution  of  the  bounda ry -val ue 
problem  which  is  coaplicated  upon  consideration  of  limitations.  In 
the  present  work  are  examined  the  systematic  special 

feature/peculiarities  of  the  solution  of  placed  problem  which  make  it 
possible  to  overcome  the  difficulty  indicated. 


1.  Setting  and  the  analysis  cf  problem. 

Let  us  examine  the  problem  of  the  selecticn  of  the  angle  of 
attack  control  of  the  vehicle,  which  is  braked  in  the  atmosphere,  in 
flight  to  minimum  and  maximum  distance  takirg  into  account  limitation 
to  the  value  of  the  complete  overload  whose  solutions  make  it 
possible  to  determine  tie  maneuverability  capabilities  of  vehicle 
(Fig.  1). 

Expression  for  a complete  overload  takes  the  following  form: 

n V ft  -t-  cl  q fS.  -V  . * 

where  g = pV2/ 2 - velocity  head  [kgf/m2];  p - atmospheric  density 
[kgf«s2/m*];  V - velocity  [m/s];  cx  - drag  coefficient;  c>  - lift 
coefficient  S - characteristic  area  of  vehicle  [m2];  G - weight  of 
vehicle  [ kg  ]. 
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Proa  (1.1)  it  is  evident  that  n cleanly  depends  on  steering 
function  rv,  and  the  limitation  in  question  belongs  to  class 

<t>, (.v,  ti)  <>  (see  [ 4 ])  . 

Let  us  assuae  that  the  aerodynaaic  forces,  which  act  on  vehicle, 
are  characterized  by  the  pciar  of  the  fern 

t s (■  1 1)  "t-C  i' , 

where  r,„  - a drag  coefficient  of  zero  angle  of  attack;  , 
parameter  of  polar. 

The  use  of  the  dependence  indicated  Bakes  it  possible  to 
sufficient  siaply  explain  the  physical  sense  of  optinum  solution. 

To  value  ry  (i.e.  to  the  value  of  angle,  cf  attack)  are 
superimposed  the  limitations: 

mm  l'y  f *'««»• 

For  the  developaent/detection  of  a broader  class  of  solutions 
and  rcle  of  the  superimposed  for  value  liaitations  the  parameters 

of  polar  and  value  cY  mln.  let  us  select  sc  that  t v(Km,x)  1/ 
would  be  inside  of  cut  |c here  K - lift-drag  ratio; 
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A The  equations  of  the  plane  aoticn  of  vehicle  in  the 

ataoaphere  take  the  fora 

1/  </  ‘ K sin  0; 

' ' m 

, q S (-  V W 

v ' hi  V ' H • n V ) | (12) 

II  V sin  'i. 

R l/co.s  0 
1 /<*(/!  ’ 

where  g = g0  R2/(R  ♦ h) 2 - acceleration  of  gravity  [ a/s2)  ; R - radius 
cf  planet  [a];  h - heiqht/altitude  of  vehicle  [a].;  g„  - surface 
gravity  of  placet  [a/s2];  9 - local  flight  path  angle  [rad];  L - 
flying  range  [ka];  t - tiae  [s];  a - mass  of  vehicle  [kg*s2/m]. 
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Fig.  1.  Key:  (1).  Boundary  cf  the  atacsphere. 
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Point  designates  differentiation  with  respect  to  t. 


Let  us  count  the  ataosphere  of  isotherial  f<  ;v  ■’*,  where 
ataospheric  density  on  the  surface  of  planet  [kgf«s*/»4)  ; P - 
cf  the  exponential  in  fcraula  for  density  [!/■]. 


Po  ~ 
index 


Let  the  descent  vehicle  coae  from  initial  state  into  final 
optimally  in  the  sense  cf  aaxiaua  cr  ainiaua  of  distance.  Let  us 
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assuie  that  in  optimum  trajectory  when  n V (,V  - limitation  to  the 

On, 

total  overload)  is  satisfied  the  condition  of  regularity  (/  ' / o 
(see  [3])-  In  this  case  for  the  solution  of  stated  problem,  it  is 
possible  to  use  the  mathematical  vehicle,  developed  in  nor ks  [1]  - 

[5]- 


Designating  those  con jugate/ccabined  tc  0,  h,  V and  L 
alternating/variable  of  variational  problee  for  system  81.2)  through 
Pi#  p 2 , Pa,  p«#  let  us  write  the  expression  of  the  Hamiltonian  of  the 

expanded  system 


H 


rv  pl/.V 


jm 


V K \ 

- l'/"”* 


! p V St  II 0 


(C.t  V-S 
\ " ’ 2 m 


\ 


^slnO  j 


p,  R l/cosft 
R h 


(1 


Si  nee  system  (1.2)  is  autonomous  during  the  period  of  the 
descent  no  limitations  are  imposed,  the  H = 0 in  all  interval  of 
action. 


In  accordance  with  [ 1]  and  [4]  the  system  the  adjoint  from  81.2) 
eguatiens  must  take  the  fera 


P, 


OH 

ox, 


HO 


On 

Ox, 
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Then 


P\  r, 


v 

IT  • h 


0 />„  UcosO  -|- 

/t’l’sin  0 


/>..  K cos  »i /?,  ^ /t 


?■<-,>  VS 

1 ’ cos  fl 

2^f  cos  0 j 

- Ps 

In, 

(ft  : /')■' 

Vs  1 

1 

'V,  u V " s 

2/,- sin')  . 

R V cos  f)  | 

P:,\ 

2 m 

R , h ' 

Pi  (R  ; •/!)-  ! <i 

i 

3-,  IMS’  , , 

' w1 r;  CI’: 

/c,  :-.v  cos  6 tfcosO  , 

* *(ia  k+i  v ! 

f,?r6'  ft  cos  4 ? l'-s  . -. 

f /’•- ^ 7,1  ~A>  /t  +-/(  > 2w 1 c*  <v  ’ 


where  X ( t)  - Lagrange's  factor,  moreover  M0<«  .V)  - 0, 

„ - "■"■')  , , 

V Ml  ♦ •J/fr)  > y 

for  systea  (1.2)  are  assiga/prescr ibmd  the  initial  conditions  v(t0)  => 
VO.  0 (t0)  * 0°,  L(tg)  * 0,  t0  * 0,  h(t0)  = fc°.  It  is  required  at  the 
specific  fixed/recorded  height/altitude  h‘,  is  sufficient  snail,  so 
that  the  distance  flight  it  Mould  te  possitlc  tc  consider  finished, 
the  provision  for  naxiaun  cr  lininun  cf  flying  range.  Inasnuch  as  8 
and  V at  the  end  of  the  flight  are  net  fix/ceccrded,  then  at  the  end 
pcint 

d',"  n\"  «*.  (15) 


_o>  „o> 

P’  P i 


Conditions  (1.5)  are  boundary  for  systea  81.4).  Pron  condition 
p4  * 0 and  p^  * -1,  it  follows  that  p4  * -1  in  an  entire  trajectory. 


Thus,  stated  proxies  is  radsced  to  tmo-poiat  boundary- value  problem 
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fcr  the  systea  of  ordinary  differential  equations. 

If  we  assign  p/0)  and  then  cf  the  condition  H = 0 it  is 

possible  to  deteraine  pa<°),  and  the  nuiber  ot  those  controlled  at  the 
end  of  the  trajectory  of  functions  p and  g3(l)  it  coincides  with  the 
nuaber  of  paraaeters,  assigned  at  initial  point.  With  this  program  of 
control  it  is  deterained  free  the  conditions 

H, , * min  itpn  L( '»  • mo x iijih  H(y  * max  npu  L <'>  — min.  (1.6) 

*»  r* 

The  constancy  of  overload  is  provided  ty  a change  in  absolute 
value  i’y(t)  in  accordance  aith  the  condition  cf 

coaaanicat ion/connection  n-  N,  and  the  sign  of  function  c,u)  is 
deterained  by  sign  pt  accordingly  (1.6).  If  in  the  process  of  notion 
along  liaitation  /i.  N in  soae  point  cy  ()  and  g ft  0,  then  in  this 
trajectory  at  the  subsequent  torque/acaent  is  possible  the 
disturbance  of  the  assigned/prescribed  liaitation.  This  is  connected 
with  the  fact  that  the  local  effect  of  control  cn  the  aaouV.  of 

g-force  is  already  exhausted  <>r'1  o)  During  the  solution  of 

\ or,  I 

boundary- value  problea  by  iterative  aethods  the  fact  indicated  is 
iaportant  first  of  all  because  in  soae  test  trajectories  can  occur 
the  disturbance  of  linitaticn.  Simultaneously  with  this  appear 
coaputational  difficulties  in  the  constructicn  cf  the  iterative 
calculation  aethods,  since  of  fl.4)  it  follows  that  \{t)  g - when 


i 
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So  that  the  iterative  process  would  not  have  the  special 
feature/peculiarity  indicated,  but  the  trajectory  of  the  expanded 
system  of  equations  with  the  disturbance  of  limitation  they 
continuously  transf er/ccnve rted  in  trajectory  without  the  disturbance 
cf  limitation,  let  us  artificially  restrict  tie  decrease  of  value 
<■„  fron  below  by  value  e.  h-e n upon  reacting  of  value  | the 
value  A (t)  is  limitedly  cn  top  .by  value  <>/  1 j Tnasauch  as  by 
hypothesis  on  unknown  optimum  trajectory  / o,  then  with 

sufficiently  small  £.  and  ; it  is  possible  to  satisfy 

boundary  conditions  (1.5).  This  makes  it  possible  to  carry  out  the 
regular  iterative  process  of  the  sclution  of  boundary- value  problem, 
without  exceeding  calculation  grid  ETsVM  ^digital  computer]. 

Page  79. 


Let  us  examine  the  systematic  special  feature/peculiavities  of 
the  solution  of  problem  taking  intc  account  limitation  nt  A During 
motion  in  open  domain  n,-  A accordingly  [1]  - [4]  A(t)  = 0;  in  this 

case  value  n corresponds  to  value  r».  determined  according  to  the 

& ’v 

principle  of  maximum  (mqffimum)  : 

uN  . (A  . 

ll.'IH  fv,  171 1'  Cy 


< 

£ 


* v ^ v min'  * y max 


/'l 

V' 
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Key:  (1).  or.  (2).  where. 

At  certain  torgue/nonent  function  " will  teccne  more  than  N,  and  by 
this  "intersection"  will  te  deternined  output/yield  to  limitation. 
Curing  notion  along  linitation  (">.  A),  in  order  to  determine 

the  torgue/nonent  of  descent  fron  linitation,  simultaneously  with 

we  conpnte  cy(pt)  without  the  account  cf  limitation  in  terms  of 
the  instantaneous  values  of  pulses  r , [see  (1.4)]  from  condition 
(1.6).  Then  the  torgue/nonent  of  descent  is  determined  by  their 
intersection  : Cyipi)^  cy(N). 


let  us  examine  the  trajectory  phase,  which  adjoins  the  end 
point,  determined  by  boundary  conditions  (1.5).  On  this  trajectory 
phase  when  L(d  - max.  carried  out  condition  « A/,  then  it  is  possible 

to  deternine  the  character  of  optimum  control  in  the  vicinity  of  the 
end  of  the  trajectory.  For  this,  let  us  examine  2*p~V'  h * 

functions  pt(»)  and  pjW  decrease  to  zero.  According  to  1*  flopital's 

pin 

rule  liinr*v=  From  condition  H = 0 ard  system  (1.4)  we  obtain 

2x  p.\  V' 11 

the  following  relationship/ratios: 


in 

Pi 


/i<")  sin 


_<n 

Pi 


n, geos  ft*11 
P'  \/in* 


in  W" 
P'  (/in 


(1.7) 


(IS) 


With  sufficiently  snail  h(*'  usually  sin  < 0,  since  the  vehicle 
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loses  altitude;  then  pt(V  > 0.  Hence  it  follows  that  pt  < 0 with  t = t. 
-A  (A  > 0)  . 

Value  p 3 ( 1 1 ) , as  P3(tt),  is  egual  to  2ero.  Then  from  conditions 
p,(t,-A)  < 0 and  0 > 0 - »/ 2 follow  p3  (t,-A)  < 0,  i.e. ■/,)  - v 

and  nrl(t,)  cy,,,  bat  when  -*  < 0 < »/2  fellows  >'VU  I,)  • and 

r ,,  , , in  accordance  with  (1.6).  In  the  case  of  flight  to 

maximum  of  distance  L(t,)  > 0 £id  cos  0 > C;  therefore  during  the 
soluticn  of  problem  on  iJ^-frsax.  the  last/latter  case  drops  off. 

However,  in  problem  on  if1'  min  both  the  case  (-v/2  < 8 < 0 and  -n  < 

6 < - */2)  are  possible. 

In  the  initial  stace  cf  aotioc  in  the  presence  of 
ccimunication/connection  '■  and  rv  it  is  a priori  unclear,  that  it 
is  batter  in  the  sense  of  the  maximization  cf  gliding  distance 
r,.  C1' or  f y rv ( k,„.,x)  with  increase  <v  increases  increases 

and  distance  with  seme  initial  Vi®’,  9°  can  decrease  as  a result  of 
premature  speed  loss.  The  selection  of  contiol  is  clearer  at  the  end 
of  the  trajectory  when  the  effect  of  the  instantaneous  value  of 
velocity  is  small  (momentum/impulse/pulse  p3  - this  the  influence 
coefficient  of  a variation  in  the  velocity  cn  distance)  and  distance 
can  be  increased  because  cf  an  increase  in  the  positive  lift,  i.e., 

<,  For  the  last/latter  phase  of  trajectory  ii*  -t  min  with  sin  0 


< 0 and  cos  « ► 0;  according  to  to  H - max.  fellows 


' V «»|*1 


1 v nil 


and 
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<v  «,  with  sin  0 < 0 and  cos  9 < 0. 

Page  80. 

In  the  case  of  the  trajectory  of  BininuB  range,  control  cy<0-  11 
with  - w/2  < 0 < 0 contributes  to  th€  decrease  cf  flying  range.  But 
if  -w  < 9 < w/2,  then  flight  is  accompl  is  he  d in  the  direction, 
opposite  initial,  and  for  tie  minimization  cf  distance  L(t,)  it  is  \ 

necessary  to  increase  the  duration  of  the  last/latter  trajectory 
phase.  This  occurs  when  'y  'y  . which  corresponds  to  the  fact  that  { 

the  lift  is  directed  against  weight. 

N 

The  made  analysis  Bakes  it  possible  tc  solve  two-point 
boundar y-value  problem  takir.g  into  account  limitation  «>.  A',  if  in 

optimum  trajectory  / o.  Case  0 is  examined  in  separate 

work. 

2.  Procedure  and  the  results  of  the  numerical  determination  of 
optimum  trajectories. 

For  practical  determination  of  optimum  trajectories  the  presence 
of  limitation  to  the  value  cf  complete  overload,  it  is  necessary  to 
rumerically  solve  system  of  equations  (1.2),  (1.4)  under  boundary 

conditions  ^1.5).  The  boundar y-value  problem  of  the  selection  of 
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initial  momen  tum/inpulse/pulses  for  satisfaction  of  conditions 

(1.5)  was  solved  by  Newton's  method.  It  was  reveal/detected  that  the 
sensitivity  of  solution  tc  changes  />"'  is  very  great,  and  surfaces 
//'  (//’  ) have  very  complex  structure.  Por  the  search  of  the  first 
approximat ion,  was  suggested  the  procedure,  based  on  with  the  aid  of 
the  principle  of  aaxiaum  the  functional  ii*1  in  QUESTION  can  be 
expressed  as  function  of  parameters  p\'\ 

» 

In  optimum  trajectory  are  fulfilled  conditions  pt(*)  = p3'»  = by  0 
and  Ii^  = aax.  (ain  1 l*>)  with  appropriate  p\ 1%,  but  at  other  values 

p\\  and  £<»>•  /!'!,  (/-“>  ■ / «,!..)■ 

In  stated  problea  only  d2L/dt2  contains  clearly  control  1 v, 
surface  lin(p i"'>*  aust  be  smooth  or,  at  least,  have  simpler  structure, 
than  p\u3 ( pi'"i)  During  search  by  the  method  cf  the  gradient  of  seguence 
pT1  which  ensures  Ii*)^  aax.  (U>) min)  will  automatically  decrease 
values  and  when  as  a result  of  the  flatness  of  surface  l.0,tp' ,%) 

the  convergence  of  the  method  of  gradiert  in  the  region  of  extremum 
it  will  deteriorate,  it  is  possible  to  refine  values  by  Newton's 

method.  This  approach  can  be  used  also  with  the  larger  number  of 
unknown  parameters. 

During  limitation  n V this  method  effectively  was  utilized  and 
provided  the  rapid  convergence  of  iterative  processes.  Nith  the  use 
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of  a gradient  method  of  the  search  of  initial  conditions  p'!\  and  of 
quadratic  extrapolation  for  Known  solutions  p^'3  (Nt)  with  i = 1,  2,  3 
where  Nt  - assigned  sizes  cf  the  g-force#  was  fulfilled  the  search 
cf  solution  /'*  i ( -V, ).  which  was  being  more  precisely  formulated  then  by 
Newton's  method. 


in  practice  during  the  numerical  realization  of  problem  by  ETsVH 

, . . V h 

instead  of  alternating/vanable  systems  (1.2)  were  utili  zed  v ^ . h j{ 

dimensionless  variables  / ' 

Page  81. 


Per  an  example  let  us  give  the  optimum  lines  of  descent  in  the 
vehicle  in  the  atmosphere  with  the  use  cf  a lift-drag  ratio.  During 


calculations  was  accepted: 


m It) 

. /'I  Mt-CrK * / , 


0,55 


k 111  il  X 

./(Hi 


■ i"v  nun  4,4 

0,302  ii/vff  cv  0,331); 


'KtU 


Mpi’K  , 


V*  0,00311  K joJh  /"’  > max; 


10  i '// 


A, 


Hi 


0,|2!»  Ktc  frK’1  i/*, 
,1  0,000137  \/m\ 

/i1"’*  31,0166  P, 

R 6371  I On  «. 


ku  - o.aotrfif,  mi 


k. 


■/. 

V' 


Key:  (1).  kg»s2/m*.  (2).  with.  (2a).  m/s.  (3).  for 
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For  the  regularization  of  problem  during  the  search  of  the  regular 
optimum  programs,  close  tc  irregular,  it  was  accepted  £ <»,oi 

During  the  calculation  of  zero  furcticns  p,, 
n v h h (".  cv  <v  they  were  determined  with  an  accuracy  to 

10'T-10  * and  values  /'I0,  were  more  precisely  formulated  during  the 
solution  of  boundary- value  problem  tc  10~4-10~9. 

Fig.  2 and  3,  give  the  results  of  the  calculation  of  the  optimum 
trajectories  U*1  max.  without  limitation  for  overload.  Let  us  focus 
attention  on  fracture  curved  n (6),  that  occurs  at.  the  value  of  the 

angle  at  which  the  teight/altitudes  of  first  and  second  maximum 

n'~  coincide.  From  Fig.  3 it  is  evident  that  dependence  ry(/) 
oscillates  about  value  ,,  <>,;{,  which  corresponds  >y  Let  us 

note  that  fluctuations  cAf)  contribute  to  damping  fluctuations  h(t) 
and  "'Jt). 
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The  obtained  solutions  were  used  for  tie  sclution  of  problem 
taking  into  account  limitation  n X Fig.  4,  gives  the  results  of 
calculations.  It  turned  out  that  during  decrease  of  N prograa  <„(/i 
changed  so  that  to  some  degree  ,it  compensated  for  the  losses  of 
distance  from  the  action  of  limitation.  Wher  the  possibilities  of 
compensation  were  exhausted,  beginning  the  noticeable  decrease  of 
distance. 
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Fig.  4 . 


Key:  (1).  a/s.  (2) 
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Key:  (1) . b/s.  (2)  . s. 
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For  the  same  parameters  of  vehicle,  were  determined  solutions  I !*) 
Bin  at  different  rates  of  entry  V°)  and  a.  Findings  were  used  during 
the  determination  of  solutions  taking  into  account  limitation  to 
overload.  Pig.  5,  gives  the  optisum  trajectory  lS^  -f  min,  for  which  N 
* 19.35  with  Vf°l  = 7900  m/s  and  * -3°. 8.  It  is  evident  that  the 
possibilities  of  decreasing  the  overload  by  local  variation  in  the 
ccntrol  are  almost  exhausted;  there  is  a maximum  size  of  the  q-force 
V-  at  which  optimum  trajectory  still  remains  regular  for  a problem 
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cn  ii1)  - Bin.  With  N < »,  optima  trajectory  will  be  irregular  and 
this  fact  must  be  considered  during  the  numerical  solution  of  stated 
problea. 
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THERMOPLASTIC  STRESSES  AND  EEFORMATICNS  OF  FUEL  TANK  IN  THE  PROCESS 
CF  ITS  EMPTYING. 

* 

H 

Y.  fl.  Harchenko. 

Is  examined  the  quasi-stationary  c ne- d imensional  tas*  oi 
thermoplasticity:  the  deternination  of  the  stressed  anu  of  tne  states 
of  strain  of  circular  cylindrical  shell  duralumin  fuel  tank  with 
horizontal  axis,  that  appear  during  its  uneven  Seating  in  the  pioce^s 
of  continuous  emptying.  Is  considered  the  Eauschinger  efrect  and  the 

variability  of  plastic  deformation  (cn  space  cccrdinate)  during 

■* 

discharging.  For  the  linear  law  of  hardening  oltained  exact  solution. 

1.  Formulation  of  the  problem.  Easic  assumptions 

The  position  of  the  current  point  cn  tie  circumference  of  the 
cross  section  of  median  surface  of  shell  is  assigned  by  angle  *,  the 
position  of  the  fuel  level  in  tank  - by  angle  >|  (Fig.  1). 


>> 
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The  determination  cf  deformations  and  stresses  is  conducted 
separately  on  the  reyicts 

0 9 < it  •]>',  (1.1) 

f - *■  (1-2) 

In  each  reyion  to  the  assigne d/pr escr i le d and  unknown  values 
appropriates  itself  index  i (i^.  = 1,  2),  that  coincides  with  the 
number  of  this  region,  to  relative  and  similar1  values  - a line  above 
le  tter ing. 

FCdNOTE  l.  Similar  are,  for  example,  boundary  values  <t  for  elastic 
solution  <h  'U  m 1 ■)  and  boundary  values  <1  for  elasto-plastic 

solution  i«l  i|m,  m i,  2.  . ■>)  BBDFCCTNOTi. 

Temperature  T:  is  asscaed  to  be  constant  in  the  i region, 
moreover  T2  > tt.  Let  us  designate  the  moduli  of  elasticity  through 

H,\  < E,  -eh:.,  pure  thermal  deformation  «<  *i  G|  <>;  “A-  1 

where  a2  - coefficient  cf  linear  expansion. 
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As  a result  of  an  abrupt  change  in  the  teaperature  cn  boundary 
cf  the  region  (1.1)  and  (1.2)  (when  ^ all  defor mations: 

elastic  *,/  f/  . irreversible  plastic  v pure/clean  thermal 
®<V  and  also  thereal  stresses  3/  and  3/  change  abruptly. 

The  complete  relative  deformation  cf  slell  e (■!/,  ■i)i  continuous  at 
all  values  *,  we  set/assume  by  that  obeying  the  law  of  tne  flat/plane 

sections: 

*(•?,?)  F ('l»)-j  (•!.)  cos?  (1.8) 


[values  F(^)  and  <I> (p)  are  subject  to  determination]. 
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/•i® 


it  is  accepted,  that 


■(+.,)  “V  0 -I  V:<». 


/>' 

E/-<  /J 


£,/> 


(14) 


and  that  in  each  region 


?)  *,/ -I  «,,(  + ■*-■  °/- 


(!■">) 


thermal  stresses  o*  are  subordinated  to  the  conditions  of 
self-fcalaqce  to  which  it  is  convenient  to  give  the  form 


2 


~ J 1 Op,  f ) £/'f  (I;  ^ j 1 (•(<,  ft  COS  fil'f 


0. 


(!-•*) 


integration  in  (1.6)  is  realize/accom p J ish ed  separately  on 
regions  (1.1)  and  (1.2). 

Equation  (1.6),  and  also,  therefore,  tfe  subsequent  solution 
they  dc  not  depend  on  radius  and  thickness  cf  shell  (usual  or  given, 
if  shell  it  is  supported).  Ihermal  stresses  taking  into  account  the 
power  stresses,  which  depend  on  thickness,  in  this  work  on,  are 


exaained. 
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The  surfaces  of  the  yield  (or  the  plastic  deformation  [2])  for 
regions  (1-1)  and  (1.2)  they  are  accepted  ic  the  form 

/(»„»„.  7\) -3,  sign  3,  n\*pt  sign?,  n,  ( 1 • ' 1 

where  /,'  the  plastic  aodule/aoduli  £ slope  tangents  or  straigi.t 
liqes,  that  connect  points  (>*.  i,,)  and  (sm>  ’/) !• 

■>„  yield  points  for  elongation  < for  compression). 

He  will  be  restricted  to  the  case  wheq  /•;)  depend  only  on 
temperature,  i.e.,  comm unication/ccqnection  (1.7)  is  linear. 


Page  66. 

The  conditions  under  which  occurs  the  resistrve  load  / 0)  ana 

discharging*  (d ,,  ii),  if  we  ascribe  the  role  cf  time  to  variable  > and 


■ 
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tc  consider  that  7',  const,  they  are  record/written  as  follows: 

l D ( l>  ,) 

J.  . / <*  r>  ill  / 0 M - , i/Stgll  3,  <>; 

y,  ' ( i)  . 0x^«  , <»■«» 

ji  a,  t'r'fH  / <•  ih6o  ecnti  y ft  n ^ sign  0. 

Kef:  41).  if.  (2).  and.  (3).  or  if. 

tGCl  NOTE  *.  Not  oi)ly  discharyiny,  tut  also  passive  loading. 
fcNDFCCTNOTE. 

Biffe  rentials  d r quotients  (ci)  ^). 

Accordingly  (1.9)  during  "discharyiny"  variables  ep/  can  be 
functions  a (but  not  ■>)•  The  reyuiriny  in  ttis  case  supplemental y 
r e lat  icnship/ratio  for  deter  linaticq  in  the  case  when  j,,..  <• 

(which  occurs  in  the  beyinniny  of  unluading  caused  by  emptying), 
us  introduce  as  follows. 

Let  the  region 

it  — i,  f » — C (ll(t) 

is  an;  of  the  regions  cf  »c rotonicit y wholly  belonging  (1.2), 

V;’ 

when  ^ 


let 
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Sot  these  values  of  '?  we  assume  that  in  (1.10) 

e/<-  (?)  r 


V(?> 


COS  /, 

' COS  'f 

( V 

«/>i  <r 

ecjm 

cos  ^ 

COS 

s 

COS* 

j COS'f 

CO 

ZP'.  <* 

cos  r, 

cos  ^ 

ec  Jin 

(l.n; 


Key:  )1).  if. 


For  values  of  •/  that  are  sulcrdinated  to  inequality  C < 
assume  that  (1-11)  it  ie  correct  crly  in  the  region 


}n  both  cases  must  he  v >0. 

Let  us  assume  also  that  dete r s i na t ions  above  the  variable  e 
satisfies  (1.8)  (condition  f = 0)  in  boundary  feints  of  region 
(1.10)  . 

Subsequently  let  us  utilize  the  designations 

f 1 /,  It  1.2);  <*(l  /,)  <•;  ^ '*  V (lid) 

2.  Eqaations  for  determining  of  P and  (j> 


L 


<i,  we 
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is  let  us  see,  in  certain  ranee  cf  charge  0)  value  irl  u. 

a are  defined  for  all  * in  (1.2).  Ecr  these  values  of  .y 

accordingly  (1-3)  - (1.5),  in  ranges  (1.1).  and  (1.2)  respectively, 
snst  lie 

, a.  ( V,  ■*)  ) 

1 e h-  | (‘j.n 

n\r(-'t)  ti>  ( p)  cos  <p!  - -t  ■ vW  | 


Eage  87. 

Jf  v = 0,  then  of  (2.1),  (1.4),  (1.6)-  it  fellows 

F(  0)  d>(0)  = 0;  , 

’>(»-  ?)  = 0’  I ('_ 

:>2  (0,  Tr)  — £2  |a2  o.,  -{-  i («)|;  1 

Cn  the  other  hand,  accordingly  (1.7), 

32  (0,  *)—  — 3,2  4 /:'i  ..  (*). 

Qf  two  expressions  for  s2  (0,  •)  with  the  aid  of  designations 
(1-131  *e  will  obtain 

•p  iW  - >•*  ( V ) ; I 

, ' , , (2  3) 


' v 
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and  therefore 


a2  0, 

~~D~ 


1 - j A j X.,;  5,,=  1 — (1  — >-_■)  X ,. 


(2.4) 


Iaicing  into  account  these  relaticnsh if /rat ios  and  (1-4) 
equalities  (2.1)  will  fce  remitter  as  fcllcts; 


o,  = e (F -\-  d>  cos<p); 

o2  --=  F -I  <P  COS  f — 1 — >"J  x,  lr  2 (-f).  I 

introducing  (2.5)  in  (1.6)  and  designating 

1 * ! 
p~— -J  P-i**  + «,*(?)!<*?: 


(2.5) 


Q=  _ J | >-2 **  f epJ(v)|cos'ft/?, 


(2.6) 


let  us  arrive  at  tlje  equations 

Fl  e ~ (e  - 1) -r  ! d >(r  ])  — -L  £ ) p; 

sin  o | , . . / 0 sin  2'1 


/'  2 ( e 1) 


f <\>  | e - (< 

l 


• i)U  ) 


2^+y. 


(2.7) 


LCC  = 78068005 


E*Gi  IjO-  , 

lO* 


Solving  (2.7)  , let  us  find 


Mt)  ' R(W\e  - 


i ) r ( >) 


+ p 


p - (r  1)^ 


sinJi  \ 

s ) 

, -I 


Q(e  1)  S'nJ 


2P(e  1) 


sill  0 I 


(2.8) 


here 


Af(-i) 


— <?(<* 

rtw 


»)sl'  (•!»); 

y-  Sill  2’1,/  ()  Sill'-' ^ 


(2.9) 


Face  88. 


let  us  designate 

3 (■',,  r >1  P I j)  — ‘l*  ( V)  COS  ■!>  --  £y  ('!•). 


(2.10) 


Variable  ®>C5)  i flays  important  role  during  the  determination 
cf  plastic  deformations. 


Eet  us  note  that  accordingly  (i.  5).,  i.e.#  for  it  is  sufficient 


L 
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si  all  •]>, 

fi C2.ll) 

Tc  us  is  known  the  £ unction  ^(i).  Let  us  designate  through 
the  rgct  of  the  equation 


If  when  < •;» , 

~H<*>  (2.18) 

and  3| (1*.  ~ 'M  in  range  (1.1)  is  the  greatest  stress,  tnan  of  (2.13) 
and  (2.11)  it  follows  that  in  range  (1.1)  e,,,  0.  Under  the 

formulated  above  conditions  of  the  fcmula  cf  this  section,  are  valid 
for  values  ■/.  satisfying  the  inequality,  entering  in  (2.13). 

3.  Purely  elastic  deformations  and  cf  the  stress  in  shell. 

This  case,  otherwise  than  in  us,  it  was  examined  in  work  [1].  We 

will  chtain  it,  after  assuming  /:,  _ oc.  Then,  accordingly  (1.8),  0 

with  any  that  leans  and  in  the  zones  of  discharging.  With 

. d, 

this  = D = JT|  92»  E * C = 0 equation  (2-£)  is  determined  F and 
but  equations  (2.5)  - with  -V  ^ *. 
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Accordingly  (2.5),  (2.6)  and  (2.10) 


AM*) 

(<’  or (■»!  m<*. - 


/ ? , _slin  27 
in'-  ‘2* 

*)  + 1 - 

(3.1) 


3i  (>.  51  *) 

e 


Curve/graph  with  e = 1-113;  -P2  P Q-=<>)  is  shown  to  by 

lower  curve  Fig.  2.  This  curve  intersects  frcm  straight  line  „ 

at  the  points  for  which  j — •?„  (n  = 1,  2,  4),  and  from 
straight  line  - x,  at  pcint  ^ — Function  e.v  (on  lower  curved) 

reaches  the  nininun  with  .j,  = .jT,  (-l2  < '>s  ' 
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c 
2 

Key:  /I).  Taking  into  acccutt  plastic  def ccaaticns.  (2).  without 
acccunt. 

Page  89. 

From  Fig.  2 and  re latio nshi p/r atio  (3.1)  it  fellows 

nith  0^:  ^ and  with  •*, « •'<  ^ 

r ’i>) 1 ^ 0-2) 

•Ji ; '!*  it 


when 
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Consequently,  with  increase  y for  its  values,  which  correspond 
(3;2)  or  which  exceed  j,  (brt  it  is  later  ic  the  process  of  emptying, 
j,)  in  so»e  parts  of  range  (1.2),  and  with  these  correspond  (3.3), 
in  range  (1.1)  or  its  part  they  can  arise  tie  condition  of  forming 
the  irreversible  plastic  deformations. 


f 


Y 


y 


a 
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The  .curve/graphs  of  stresses  fcund  from  formulas  (1.4), 

(2.5),  (2.8)  when  <>,  are  shewn  by  Ircken  lines  uy  curves  on 

Fig.  i.  They  testify  that  with  increase  of  j in  range  (1.2)  occurs 
the  discharging,  and  in  ranee  (1.1),  beginning  from  ‘load 

beyend  yield  point.  From  them  it  is  evident  alsc  which  a t ((■;,,  r — ■ j i 
ir  range  (1.1)  has  great,  ard  o2  <o,  r.—  ',i  in  ranee  (1.2)  - a small 
value.  Then,  as  it  follows  from  Pig.  2 and  relationship/ratio  (J.1), 
with  . -'j,  and  with  .j,,  occurs  the  inequality- 

13:.  (•’.,,  r.  'j)  3*  . O-4) 

a zone  s .•  and  _ ..  _ r either  their  part  when 

these  zones  or  their  parts  belong  (1.2),  are  kncwingly  tree  frcii 
plastic  deformations  = 

The  made  analysis  cf  elastic  solution  substantially  facilitates 
finding  the  zones  of  plastic  detonations. 

4.  Elasto- plastic  deformations  and  of  the  stress  in  shell. 

Firs^case.  Let  us  expand/develop  the  solution,  obtained  earlier 

'r 

for  values  of  J>  satisfying  inequality  (2.12),  at  which  6 First 

let  us  construct  s ..  Accord ing  to  results  Section  1,  it  is  possible 
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Jfi 

tc  expect  the  f oca  ation/e  ducation  c£  three  ranees  of  non  at  onic  ity  . (•* ! 


and 


< ’ * 


two  zones  in  which  w ten 

’h 


» 


■?:  < * < » 


•Vi  and  wnen  n-fs  ^ * 


(4.2) 

Let  us  designate  through  4m  (a  = 1,,  2.  4)  values  *,  that 
satisfy  condition  sfl„  .( ■ ~ f ) = 0,  through  ^ 4 to  the  conditions 

sp  i (••»:..  * 'tO  -=  0 ( " £+  ('Vi)  = xl.  ^ ^ 

Key:  (1).  and. 

through  V-  v.«  — to  conditions 

**  ('>:<)  = min  s+  (4)  (<V,)  = 0.  (4.4) 

«y 

Fcr  determination  '>m  let  us  net*  that  g — 1 when  r ,%  > ,, 

2 S ep  *J  \TC  .m'  u> 

and  then  according  to  tfe  seccpd  ecuaticns  frea  (2.4)  and  (2.5),  we 
will  obtain 

i\ (■'jJ  = %,  (m  — ) , ‘2,  4).  (4.5) 


* 
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finding  from  (4.5)  makes  sense  only  when  >0  | *.  <)..  < 

when  real/actually  there  exists.  When  < 0 everywhere,^ 

As  follows  from  Secticij  1 in  the  first  cf  ranges  (4.1)  must 
appear  dependence  - <'f >,  ccrrespo cdi ng  tc  the  first  case  from 
(1.111,  ■oreover  for  it  accordingly  (4.1),  (2.3)  and  (1.11) 

0;  t v,;  , -E 

Page  9 1. 


Jn  the  second  of  ranges  (4.1)  it  is  necessary  to  expect  the 
dependence  e, ,*(?),  corresponding  tc  seccrid  case  (1.11),  Lor  which 

^ 'i  and  acccrdirgly  (1.5),  (2.4),  (4.4)  and  (1.11) 

E:,  (‘{'l*  ~ ' + *2  *5  — V- 


®n  the  other  hand,  frca  (1.8),  (1.11).  and  (2.4)  it  rollows 


o.  ('(*3>  ~ I**)  1 1 ' 


fjj 

F, 


Thus,  taking  into  accctnt  (1.13) 
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In  the  third  of  ranges  (4.1)# 
(1.11),  for  the  first  case  ; = 


acccrdiqc  tc  tnat  presented  and 

•'jt\  v-  E/»  ? ( Va)  * i ' 1 i : i 


Consequently,  for  that  lie  at  ranges  (4.1)  and  (4.2),  when 
these  ranges  or  at  least  their  parts  fcelonc  (1.2),  function  e„ (f) 
necessary  to  assign  thus: 


cos  'j. 

*r  2 " l?  j 

*;s=o 

Zr  ..  (/,  - s:1) 

e;,  , = ) ,(x2  — s,) 

e,.2'  () 


-(  COS  9 

cost, 

cos  ■}>.,  -f 
cos  y,  — 
COS  ’j,  + 
CMS 


uptl  " 


cos  9 

CMS  \ 

cos  9 
COS 


o) 

nim 

W 


II  pi! 


'?l 

'b  ' 

'•'■i 

If* 

’h 


: (4.c>) 

i 


Key:  41)  . with. 


Then  consecutively  it  is  possible  to  find  everything  'i»„.  and  s,, 
which  will  make  it  possible  to  ccnpleteiy  sclve  our  task  tot  () 


It  is  real/actual,  0 ■?  ' "f\  Then  in  field  (1.2)  we  assign  7n, 

according  to  first  equation  (4.6)  and  cd  (2.6)  we  compute  P = P x (j ) 
and  Q * Q»  ('/)  [P  and  Q it  is  convenient  to  assign  the  inaex,  which 
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corresponds  to  the  nusber  of  cons  iron  (4.6),  used  for  tneir 
detersinat  icn  ].  On  (2.8),  we  find  F acd  <]>,  cn  (2.10)  (y)  and 

free  (4.5)  with  a = 1 we  determine  Further  similarly  are  located 

/’.  <l>,  =■  a Is  c , from  (4.5)  that  satisries  the 

condition 


Vs* 


(1.7) 


With  observance  (4.7)  we  find  <*>,(•».  /'.  <l>.  and  from 

equations  (4.4)  and  £|  If  it  seems  that  E /,,  then  we  compute 
■(  :>).»  Q« ('»  and  so  forth,  and  iron  (4.5)  Finally,  through  p5ty), 

we  find  7.  <l\  e+  also,  frem  equation  (4.3)  If  it  seems 

that  it  region  (4.7)  there  are  no  rccts  of  equation  (4.5),  then  it 
■eans  that  min  it  lie/rests  above  straight  line  E;  *•.  and 

therefore  \p%  0 everywhere,  besides  the  first  of  regions  (4.1). 

If  It  turns  out  that  e - x2j  then  the  mentioned 
straightline  is  tangent  to  ^(yl  in  the  point  of  the  .y  , •,  ~ t 

everywhere  except  the  first  area  (4.1). 

Incidentally  are  calculated  stresses 


Second  case.  If  0 > y,,' 
defoliations  «„,><),  in  the 
region  (1.2).  For  (1.2)  they 


then  iq  region  (1.1)  will  arise 
process  of  emptying  converting  into 
will  become  the  deformations  of 


discharging  >0). 


that  seans  that  in  the  cospletely  esptied  tank  when 


>2 


(4  H) 


COS  j , cos  y 

1 cm  J* 
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According  (4.  8) 

*..)=  •,«<*- 'W  ~0’ 

■,.10)- >-v,(°)=  ):■  l3>’  0)  3"'- 

(4.9) 

Let  us  designate 

H-10> 

and  let  us  note  that  accordingly  (1.5)  and  (4.9), 


5,  (-.  U) 

€ 


(4.11) 


from  (4.11)  and  (4.9)  taking  into  acccmt  (1.13)  we  find  that 
after  which  (4.8)  it  is  copied  as  fellows: 

■ hen  0 f " _ 'V- 

. V COS  'j5  I-  cos  * 

*„■:  '\  , • 


(4.12) 
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*cw  for  the  completely  emptied  tuk  (•>  = ■)  we  rignt  tc  use  the 
second  of  equations  (2.5),  and  means  ty  formulas  (2. 6) -(2.  8).  Taking 
into  account  (4.6)  and  (4.  12)  the  defcrnaticn  ,,,(■+)  is  determined 
everywhere  and  computations  using  formulas  (2.6)  do  not  encounter 
difficulties.  Computing  ty  (2.6)  P(»)  and  C (»)  « and  then  on  (2.8) 

F (it)  mod  «i>  (w)  m .in)  wf  will  obtain  linear  equation  relatively 

(when  vr.  < * must  <be  Let  us  return  tc  case  by  < v.  The 

stresses  #,  ((•{«,»  #)  , examined  in  Section  3,  and  means  in  the  case  iq 
question,  with  decrease  h,  they  decrease.  Therefore  with  V 

generally  speaking,  not  in  entire  region  (1.1) 

lot  •(>..>•>»)  he  that  value  of  parameter  i,,  at  mhich  iq  (1.1) 

*, * <’f ) \ / i>;  moreover  t/M<0)  = °-  Then  for  + with  the  aid  of 

(U5)W  ( 1. 8) , (2.4),  (1.13)  we  obtain 

»i  (■('•  *'(/■'  -*  4*  cot  9)  t-  i ; I (4.13) 

J . ( , t))  - t *1*  (08  I — A>j  *2  — 

moreover  in  region  f-  ; ^ ^ by  beiqg  part  (1-2),  7p,(f)  we 

assign  on  formula  (4.11).  Bith  the  aid  of  (4.13)  and  (1.6),  again  let 
us  arrive  at  equations  (2.7),  if  we  in  thea  replace  e by  e and  to 

asmumd 

I r. 

K 

sin  <1  2 (•  — 

<2  ( ')  - r.  • -K  \ I'***  I 

r* 


(4.M) 
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Consequently,  as  solution  (2.7)  for  the  case  in  question  they 
also  serve  (2.8),  (2.9),  if  we  in  thea  replace  e by  e,  and  P and  Q to 

calculate  according  to  (4.14).  Let  as  acts  that  this  solution  aith  > 
= y leads  to  the  sate  values  of  F (*)  and  *t<(a),  which  enter  in  (4.10). 

the  first  of  equalities  (4.13),  aith  4 ~ C and  «|  i^,|i  0)  3jit 
reduces  to  the  equation 

M!-.) + ♦<;•«>  v <4-l5> 

By  the  aethod  presented  it  is  possible  to  obtain  equations  tor 
F,  o>.  i,  and  in  the  case  k ot  >«  Howevex,  sith  ...  close  to 
(which  usually  occurs),  this  ipterval/gap  change  it  is  possible  not 
to  exasine. 

Page  93. 

the  given  fornulas  wake  it  possible  tc  solve  states  piotlca  in 
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criterial  Cora.  Criterion  *,  &2  takes  it  possible  even 

before  calculations  to  judge  the  possibility  of  the  eaergence  c£ 
plastic  de for aatio ps  in  tag*.  Critericf  £;1  solves  the  problem  about 
the  naaber  of  plastic  flow  areas  during  discharging. 


Calculation  ifas  carried  out  by  the  author  for  stresses  and 


straiss  of  teak  frqn  D 1611  vith  f clicking  k)  initial  data: 


3Q0°K| 


ft,-  250°;  *,0,  8,53;  °3 , ■*)/>  mc/mm2; 

= 6540  kzc  mm2,  e --  1,113;  E J - lN'.K)  kzc  mv: 


o,j—  25  kzc/mm 
ti\{  1615  kzc  ! m m ' . 


T,  = 


AlC  - K 


the  results  of  calculation  are  given  tc  Fig.  2 and  3. 

the  nunerical  values  of  the  stresses  and  strains  in  regxon  (1.2) 
satisfy  the  conditions  of  discharging  (1.9).  Ibe  account  of  plastic 
detonations  leads  to  a noticeable  reduction  in  the  stresses. 
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HEIHOB  OF  SUCCESSIVE  AE FBCX IH ATIC NS  IN  EHOEIEM  CF  TRANSIENT  CREEP  AND 
CF  NONLINEAR  ELASTICI1X. 

I.  I.  Pospelov. 

In  work  is  given  the  method  of  the  solution  of  the  problems  of 
theory  of  creep  that  is  the  development  of  the  aethod  of  elastic 
solutions  [1].  Method  allows  physically  the  nonlinear  task  of  the 
theory  of  creep  to  reduce  to  the  sequence  ci  linear  tasks  and  tc 
describe  the  redistribution  of  the  stresses  in  construction  in  the 
process  of  creep.  Unlike  work  [2]  ccsplebe  strain  is  represented  in 
the  fpra  of  the  sua  of  instantaneous  deformation,  by  nonlinear  form 
voltage-sensitive,  and  creep  strain,  nonlinear  voltage-sensitive  and 
tiae. 


The  behavior  of  aaterial  during  creep  is  described  by  the  theory 


cf  fie*.  1a  given  an  example  of  nuacrical  ccapatation 
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1.  Fundamental  equations  c£  the  theory  of  creep. 

It  is  assumed  that  the  components  o£  the  deviator  of  complete 
strain  etl  are  accumulated  frca  the  components  cf  the  deviator  of 
instantaneous  deformation  e'n>  depending  only  cn  stresses,  and  creep 
strain  Pn<  that  voltage-sensitive  agd  the  tine: 

eu  = e'ij  + Pu  u’ i *•  3)>  0-') 

where  here  -,/  a strain  .tensor;  i(/ — 1 with  i = j and 

*'v“w  with  i ^ j. 

f or  describing  the  process  of  creep,,  are  .utilized  tne  equations 
cf  the  type 

Pt)  - 7f  Stj’  ( 1-2) 

where  sij—  a stress  deviatcr;  plt  - derivative  of  the  strain 
deviator  ef  creep  on  the  nodified  tine  r(t) , which  is  the  function  of 
the  physical  tine  t;  3«.  P*-  the  stress  intensities  to  of  the  rates 

cf  creep  strain,  moreover 

3 
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(h«re  ’</ — stress  tenser;  • - average/iean  stress). 


Eage  95. 


the  relationship/ratio  between  intensities 
cf  creep  strain  aqd  tbe  tiee  is  accepted  in  the 


o£  stresses  and  nates 
following  form: 


p»  -/<5 *)• 


(1.4) 


During  numerical  ccnputations  is  utili2«d  /(3„)  ,l3».  Function  f (■>„) 

and  the  aodified  time  r - r (t)  are  determined  frem  tne  grid  of 
curves  cf  creep,  obtained  experimentally  with  elongation  under 
conditions  of  constant  temperature  [2]. 


The  components  of  the  deviator  of  instantaneous  deformation  are 
determined  by  the  eguations 


T 


(l-r>) 
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Coaaunication/connecticn  between  intensities  of  i nsta ntaaeous 
defcraation  and  stresses 

c'n  ?(«■)  (I‘>) 


let  ua  accept  in  tbe  fora 


(1.7) 


where  \>  - shear  aodolus;  a,  ••  - constants  cf  aaterial  which  can  be 
deterained  froa  diagraa 


Coaaunication/connecticn  between  .the  aterace/aean  stress  3 = -,u 

and  the  awerage/aean  strain  £ }- tn  ia  expressed  by  the  equation 

z-KH,  <!•«) 

/• 

where  A'  ^ ( is  (here  E - Young's  aodtlus,  Poisson 

ratio}. 


Frca  equations  (1.1),  (1.2),  (1.4), -(  lit)  we  will  ootain 


deu  ..  < d fhH) 
d~.  2 d-  a.,  s l( 


< /K) 


( 1 .it) 


Equations  (1.9)  describe  the  tehavior  cf  aaterial  both  during 
instantaneous  nonlinear  defcraaticc  and  transient  creep.  They  not 
linear  and  their  use  during  the  solution  of  prcfelens  is  connected 
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uith  ocnsiderable  mathematical  difficulties. 
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Equation  (1.4)  is  convenient,  isclatiqc  linear  part,  to  present 
in  the  fori 


p„  1 - Ti)- 


(1.10) 


where 


v ' ' - !M. 


(1-11) 


it  is  pcssible  to  subordinate  to  condition  C ^ 'i  < 1 by  selection  ot 
Key's  constant 


[)  (,,'d  (*J-  > 0; 

Oh  ihii  "•'h 


(1).  «ith 


here  o„„l3>  certain  conditional  nuiter- 


I n calculations  D,  it  was  calculated  from  the  formula 


f (°m  max) 


(1.12) 


Analogous  to  aquation  (1.6)  we  represent  in  the  form 


<>H  “0  +*»). 


(1.13) 
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where 


f (3„)  - 


!‘i  — I1  !li  f ?«\m  ' 

11  ' 3°  l O0  J 


(1-14) 


it  is  possible  to  subordinate  to  condition  *1  ' « ^ o by  the 


selection 


I + ? fi 

II  3"  ^ 


’ ( m “ V 

• V ) 


for  pi  it  is  possible  to  accept 


1 ;1  / °M  maV  \*—  1 ’ 

p 1 3°  l J,‘  I 


which  corresponds  to  secant  module/mcdulus  to  the  diagram  £?m~3m  for 


mat*  Or 


ii  + " a"(  ,»  ) 

that  corresponds  to  the  tangent  modulus  on  diagram  *h~8h  for  7Hltla>. 


{n  the  region  where  communicaticn/concecticn  between  the  stress 
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and  the  strain  beats  linear  character, 

'<  /,  \m  -I  ^ 

t j U |i  n,  ,i. 

3"  ^ 9"  y 

Key:  11).  and. 

Taking  into  account  (1.10)  and  (1.13)  equation  (1.9)  let  us 
present  in  the  forn 

S//  + 3;* ,/>*„  — -M-i  <¥>,>-  0-17) 

cr  in  the  integral  fora 

Sir  -K//  * ///+/<;-  ^I  I8^ 

uhece  ;iz  a linear  operator. 


3ji*  w><  ~t„i  j’ 

(1.19) 

P,e  rJ  ■3'  r,,z'- 

(1.20) 

■ 2 Dr‘*‘<  < •2Mt(S"U,,: 

(1.21) 

,(t„)  2,M  f,y(t0)|c 

(1.22) 
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Iguations  (1.18)  describe  the  kehavioi  of  aaterial  uoth  during 
instantaneous  nonlinear  defornation  and  during  creep.  For  the  elastic 
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medium  for  which  is  — ll.  - (l-  -s/,lT»>  '-'M’t/Mtzc*  equations 

( 1-  19J- ( 1. 2 1)  we  have  /,;  f Lj  ",  n j».  § Ecr  the  linear  medium, 

descrited  by  Maxwell's  modal,  *x  - '.  = C,  t = 0 and  from  equations 

(1.181  ue  will  obtain 

s i i “!“'//• 


Srca  comparison  (1.18)  and  (1.13)  it  is  evident  that  the 
functions  t,  ago  w wten  = p characterize  the  deviation  of  the 
(rcpecties  of  body  froa  the  properties  cf  Paxwell's  model  with  the 
■ edified  time  r = r (t)  . 


If  we  place  D = 0,  then  equations  (1.1b)  will  describe  the 
plastic  deforaation  of  material  during  the  active  process  cf  leading 
or  the  nonlinear  elastic  tetavicr  cf  material.  In  this  case  t,  it  can 
serve  as  the  parameter  of  the  loading  of  cccstr uction. 

2.  Eqaations  of  the  theory  cf  cree(  in  displacement/movements. 


Froa  the  equations  cf  static  equilibria,  expressed  into 
stresses, 

^'  + P/V-,1  (/.y-1,2.3).  (2-1) 


(wherm 


F,  vector  of  aass  forces). 


1 (efu, 
'2  \0.x/ 


du , 
d x , 


from  >Cauchy  formula  - 

(2.2) 
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(here  ut  displacement  vector),  and  from  equations  (1.  Id)  we  will 
obtain  the  equations  of  equilibrium,  expressed  in 
displacement/aoveaents,  in  the  form  cf  Lama: 


nv  u i + (:*•  W > 


-(f >/  /-  *,l. 


where  ).z  linear  operator. 


>z  - ( K 


5 


^ ( /<  > ' f > / ' 


The  boundary  conditions,  expressed  in  the  stresses 


Jr-  F < 


(where  lt—  direction  cosires  of  normal  to  the  surface; 
forces,  assiqn/prescrihed  on  boundary  surface  with  the  use  of 
relat ienship/r atios  (1.18)  and  (2.2)  they  are  converted  to  the  fori 


A),,  <>u . \ 


where 


•l*-.  (ft,  r 
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* 


i 


h 


for  the  elastic  medium  for  which  <•,  n,  s(/(t(1)  (t0),  f)  (), 

and  consequently,  <t>,,  <1>,  o,  n r,  / the  system  of  equations 

(2.3)  and  (2.6)  is  reduced  to  the  known  equations  of  Lame  and 
boundary  conditions.  Equations  (2.3)  and  (2.6)  formally  coincide  with 
the  equations  of  the  theory  of  small  e lastc-plastic  deformations,  to 
which  let  us  use  the  method  of  elastic  solutions.  \ 

Ihe  procedure  of  the  calculation  of  the  stressed  ana  state  of  I 

strain  of  body,  which  is  located  under  conditions  of  transient  creep, 
by  the  method  of  successive  approximations,  consists  of  rolloving. 


Bor  determining  the  first  approximates,  me  set/assume  that 

*he( 


«<")  o.  Then 


in) 

"n  / 


/.?  ",  /„,  c pf 


*fi, 

Ox . 


f ill  be  the  known 


functions,  determined  by  initial  conditions  and  which  are  the 
solution  either  of  elastic  or  nonlinear  elastic  problem.  Equations 
(2.3)  and  (2.6)  become  the  equations  of  the  linear  theory  of 


viscoelasticity.  1th*  solution  of  these  equations  with  the  specified 
initial  conditions  «»"  we  take  for  the  first  ap froximatioa . From 
equations  (2.2)  1st  us  find  that  s',/’,  e(,)\  c i equations  (1.18) 

= ' f,r  fros  equations  (1.3),  (1.11),  (1.1k),  (1-21),  (1.20), 

(2m  4)  r (2.7)  - respectively  ’l,'’,  V",  f‘}\  cl);'1,  and  for  the 

determination  of  the  second  approach/approximation  we  again 


A 
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utilile  a system  of  linear  equations  (2.3).  and  (2.6)  witn  the 
converted  eight  sides  which  can  be  interpreted  as  fictitious  unsteady 
external  forces.  Ccntincicq  this  the  process  of  the  solution  of  the 
sequence  of  the  tasks  of  the  linear  theory  ct  viscoelasticity  with 
the  introduction  of  fictiticus  external  forces,  we  will  obtain 
required  the  accur acy/precision  of  results.  The  rate  of  the 
ccpvexgence  of  the  sequence  of  apprcach/apptcxiaations  will  decrease 
in  the  course  of  time-,  therefore  fer  the  purpose  of  the  savings  of 
machine  time  during  the  calculatioE  that  stressed  and  of  states  of 
strain  one  should  utilize  pcint-by-pcint  method,  i.e.,  at  each  space 
cn  time  to  solve  system  of  equations  (2.3)  and  (2.6)  with  initial 
conditions  and  constants  £ and  p, , calculated  cn  the 
precedinq/previous  space.  Hcrecver  fer  conditional  value  31)1Ilax  is  used 
the  maximum  value  of  the  stress  intensity  it  tedy,  calculated  in  the 
precedinq/previous  space:  U,  x,  x,  x'„),  mrltiplied  on  a,  where  a 

1.  On  given  one  a in  the  process  satisfying  the  inequality  jwu-.  a. a-, if' 

we  process  of  the  space  A-.,  for  U ...  o,  o t, <' T, 

In  general  form  a questicn  concerning  the  convergence  ot  the 
approach/approximations,  obtained  ty  the  set-ferth  method,  requires 
supplementary  investigations.  The  examined  ielcw  example  of  numerical 
computation  testifies  tc  the  sufficiently  high  cate  of  the  sequence 
cf  a p | r cach/approximat iens . 
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3.  Elongation  of  rod  with  a given  rate  c£  deformation. 


Let  us  examine  the  red,  which  is  deforied  with  speed 
Then  taking  into  account  the  incompressibility  cf  material  3„  — ’i  si  ^=‘ 
where  and  E|  - stress  and  complete  strain  of  rod,  from 
equations  (1.10),  (1.11),  (1.13)  and  (1-14)  we  cbtain 


3;i,C|  <lli  ^3iTl  — //r  (V">- 


The  solution  of  eguaticn  (3.1)  for  the  k iteration  can  be 
presented  in  the  fora 


where 


o'*'(t)-  C — -,.||(,0)  1-  (o' ' l)-i- 

X 

t :i !»,/->  f (v.)*"1  — (5|«,)*-|e'1:,|W'-  ’ I 

x„ 


a,  — |i  3 i*i  / °i  v 

' H + 0" 


/(®l) 

Dh 
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This  task  allcw/as  suites  the  exact  solution,  expressed  in 
cuadratures, 

t, 

T-T.+  \ ■?(5|)*l.  <:V3> 


where 


•M3i) 


W 


r,  - /(’,) 

/(». 


/Is; . 


(3.4) 


The  results  of  the  calculation  cf  a change  in  the  stress  in 
tixe,  obtained  on  formulas  (3.2)  and  (3.3L  with  A^.o.in  io  '■(  h*‘ , 
A-  9.9  p = of  2117  kgf/aa*,  «o  - 46  kgfyaa*  <and  n = 3.66,  coincide  and 
are  represented  in  Fig.  1.  For  accelerating  the  process  of  the 
convergence  of  approach/approxiaaticns,  the  tiie  interval  in  question 
was  divide/narked  off  into  the  cuts  tilth  space  tr  = 0.02,  in  each  of 
which  the  stress  was  deterained  frca  focaula  (3.2)  with  initial 
condition  and  constants  D and  p|  ty  deteraited  equations  (1.12)  and 
(1.151  or  (1.16),  calculated  on  the  prccedii g/previous  cut. 
Conditional  constant  nunber  sma,  was  determined  froa  formula  = «(t0). 


Has  accepted  a 


1.2 
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Ke > : /I),  kgf/aa2.  (2).  r [sin]. 

Fage  ICQ. 

In  the  process  of  count*  satisfying  inequality  3m„,  we 

deteraine  the  value  of  space  Ar  via  its  f ragaentation  in  such  a May* 
that  aculd  be  observed  the  conditions  by  - 1 < w ^ 0*  0 ^ v < 1-  The 
results  of  calculation  testify  to  the  sufficiectly  high  rate  of  the 
convergence  cf  the  sequence  cf  ap p rcac h/a p p ioxi «ations . For  the 
calculation  of  stresses  uitfc  an  accuracy  tc  fifth  significant  digit* 
is  not  required  aore  than  eight  apprcach/apprcxiaations. 
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If llCt  OF  HEATING  ICB- F f ESS U BE  GAS  IN  A SHOCK  TUBE  ON  INCREASE  IN  THE 
ATTAINABLE  TEHPERATUBE  Cf  STAGNATION. 

G.  1.  Grodzovskiy. 

Ie  investigated  the  problem  of  an  increase  in  the  attainable 
teaperature  of  stagnaticn  cf  gas  flew  in  shcck  tube.  It  is  shown, 
that  fee  this  purpose  is  advisable  heatipg  lew- pressure  gas.  Is 
carried  out  the  analysis  of  the  effect  of  beating  low-pressure  gas  on 
the  increase  of  the  attainable  tenperature  cf  stagnation  in  flew 
behind  shock  wave  for  tbe  case  cf  a cne-diaphraga  cylindrical  shock 
tuke . 

To  gas  dynamics  of  flows  in  shock  tulies  devoted  large  number  cf 
investigations  (for  exaaple,  see  [1]  - [3,]>.  Primary  attention  in 
these  investigations  is  deveted  to  the  problem  cf  the  achievement  of 
maximal  relation  to  the  rate  of  the  shcck  wave  0 to  the  speed  of 
scand  in  the  motioalesa  gas  before  wave  a4 : 

M,  ~ P.  . 

"i 
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where  V*uift,  l static  teaperatrre  cf  low-pressure  yas  in 
front  of  the  nave,  *,  and  P,  - adiatatic  icdex  and  the  ja  s corstant 
cf  this  gas. 


for  the  simplest  cylindrical  c qe-c  iaphi ay ■ shock  tube  (Fiy.  'll 
■axiaua  value  of  n uabec  fll#  as  is  kncwn,  is  reached  at  au  infinite 
pressure  differential  cn  the  diaphragm 


M, 


(v,  t- 1 ) «»;_  _ i 1 i /*.  r< 
(».,  i)  nt  *i  i ' 7 1 K\  7 i 


(i) 


where  by  index  4 are  noted  the  paraaeters  cl  high-pressure  yas. 


In  accordance  with  relationship/ratio  |1)  nuaber  Mi  increases 
with  an  increase  in  the  teaperature  of  the  t igb-pressure  gas  T4  and 
with  an  increase  in  its  gas  constant  E«.  Therefore  considerable 
attention  in  works  [ 1 ],  [3]  it  is  given  to  the  probleas  of  heating 
high-pressure  gas  with  the  cse  cf  high-pressure  gases  with  light 
aclecilar  weight.  Froa  these  positions  low-pressure  gas  (Tt)  was 
exaained  cold. 

By  us  is  investigated  the  protlea  of  ap  increase  in  the 
attainable  teaperature  cf  stagnaticq  T0  of  gas  flow  in  shock  tube.  It 
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is  shgwn,  that  tor  this  purpose  is  advisable  heating  low-pressure 
gas.  By  calculation  for  fiqal  pressure  dif 1 erentials  on  diaphragm 
this  effect  was  independently  previously  fund  by  n.  1.  Khvostov. 

Lower  cn  the  basis  of  analytical  solution  is  obtained  the  universal 
dependence  of  the  attainable  teaperature  of  stagnation  or  gas  flew  in 
shock  tube  on  teaperature  of  lcw-presscre  gas. 

let  us  give  analysis  for  the  simplest  diagram  of  shock  tube  for  ‘ 

the  perfect  gases  (see  fig.  1).  The  cbtaincc  results  can  be  common 
for  the  cases  of  sore  compound  circuits  and  for  gases.  \ 

V 

Eage  102. 

The  temperature  of  stagnation  10  of  gas  flew  in  region  2 (behind 
shock  wave),  it  is  logical,  it  depends  on  the  rate  of  flow  u2  = u3  = 

V and  the  parameters  of  low-pressure  gas.  rite  limiting  value  of  the 
velocity  of  the  gas  flew  V is  reached  at  an  infinite  pressure 
differential  cn  diaphrags,  value  i ,„„K  depends  only  on  the  parameters 
of  the  high-pressure  gas: 

V'm.*  | "■  . • UKKiT,  (2) 

fee  the  fixed  value  of  V,  we  ccee  to  the  task  of  tue  maximum 
attainable  temperature  of  stagnation  lo  in  flew  behind  shock  wave  in 
the  gas,  compressed  dr i ving/moving  at  a rate  of  V by  the  piston  (role 

i 
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c £ pistcn  performs  contact  surface,  see  Pic.  1). 


Eton  the  equations  of  shock  wave  fropa^aticn,  it  xs  possible  to 
obtain  the  following  expression  fcr  the  relative  rate  of  flow  behind 
tfce  shock  wave: 


v 


2 \ : *■  M'J.  | 

Q I 


i * i 
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Fig.  1.  1 - high-pressure  cbaafcers;  2 - low-pressure  chamber;  J - 
diaphraga;  4 - shock  wave;  5 - contact  surfece;  6 - centered 
rarefaction  wave. 


Page  103 
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A jump/drop  in  the  static  temperatures  T21j  on  shocK  wave  is 
determined  by  the  known  relationship/ratio 


r 

/, 


i t 


M,  i I 


M 


(I) 


Bespectively  the  ratio  of  the  unknown  temperature  or  stagnation 
la  in  flow  behind  shock  wave  to  the  static  temperature  of  the 
lcw-piessure  gas  Tt  is  possible  write  in  tie  term 

/< M,)  ~ ‘ 'I*  (M,).  <5 1 


whence  cne  should  expression  for  the  dimensiccless  value  of  the 
teapenature  cf  stagnaticn  it  flew  the  shock  wave 


T„  /<M|)  *,  — I 

t ’VC'i  ATA-.i  '!•  (M,)  2 


(0) 


One  should  consider  that  according  to  tguation  (3)  to  each  value 
cf  pasameter  H , correspcnds  the  specific  value  cf  the  relative 
velocity  of  flow  behind  shock  wave  (piston  speed)  V/a,. 


i 


fig.  2,  gives  a change  in  the  diaecsiciless  value  or  the 
teaperature  of  stagnation  T0  depending  cn  tie  relative  speed  v/a,.  it 
is  evident  that  at  the  assigned/prescriled  value  of  velocity  v cf 
flew  behind  shock  wave  (piston  speed),  an  increase  in  the  temperature 
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cf  the  low-pressure  yas  T , leads  tc  an  increase  in  the  attainable 
temperature  of  stagnation  Tc.  Ihus,  for  instance,  if  parameter  Mx  wa 
located  ip  the  range  10  N<  M j N<  ^ 20,  then  curing  heatinj  cf 
low-pressure  gas  in  shock  tube  it  is  possible  it  is  more  than  tc 
twice  raise  the  temperature  of  stagnation  lc  in  flow  behind  sIicck 
wa  ve . 

FEE  ERE  KCES. 
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EASE  EBESSU  RE  AFTER  RECT ANGILAR  STEES  R IT H IIEIERENT  RATIUS  OF  HEIGHT 
1C  Ifcl  HIDTH  OF  STEP. 

G.  K.  Lavrufchin. 

Ate  giveq  the  results  cf  the  experimen tal  investigation  of  base 
pressire  after  rectangular  steps  with  different  ratios  of 
height/altitude  to  the  widtb  cf  step  at  subronic  and  supersonic  speed 
cf  external  flow. 

Is  given  the  empirical  formula  for  determining  the  case  pressure 
after  rectangular  steps,  vfcich  considers  the  effect  of  tae  relative 
height/altitude  of  steps  and  Mach  number  of  the  incident  flow. 

At  present  is  investigated  in  sufficient  detail  base  pressure 
after  a x is ym metric  and  flat/place  steps  [1]  - [3].  As  showed  these 
investigations,  in  the  case  of  turbulent  bcindary  layer  base  pressure 
after  the  axisymmetric  step  higher  than  base  pressure  after 
flat/plane  step  on  10-15o/o  with  sutsoqic  apd  cn  30-40o/c  at 
superscnic  speeds  of  external  flow,  however,  the  information  about 
base  pressure  after  the  steps  cf  scie  '.'trassiticnal"  forms  - oval. 


rectangular  with  different  ratios  of  he ig h t/a 1 t itude  to  the  width  of 
step  and  the  like  - is  virtually  atsent. 

Theoretical  studies  of  tase  pressure  after  such  step:*  encounters 
great  difficulties  due  to  the  need  fcr  considering  transverse 
overflowing  for  separation  acnes. 

Article  gives  the  results  of  the  experimental  investigation  of 
tase  pressure  after  rectangular  steps  with  cifierent  ratios  of 
height/altitude  to  width.  I rvestigaticns  were  conducted  with  Mach 
numbers  of  external  flew  C-t4-2.7£. 

Reynolds  numbers  changed  in  the  range  te  = 4.5-14*10*  as  a 
result  of  a change  both  the  length  of  the  model  and  the  total 
pressure  in  external  flew. 

Here  investigated  the  models  of  two  combinations  (wedge  - 
parallelepiped  and  cone  - cylinder).  The  basis  cf  model  was  central 
drain/vented  plate  1,  establish/installed  on  holding  pylons  2 (Fig. 
1)..  With  the  aid  of  interchangeable  rectangular  plates  J and 
segmental  extensions  4#  fastened  to  central  plate,  they  were 
ccmpose/collected  the  body  cf  combinations  indicated  above. 

The  width  of  rectangular  steps  remainec  fcr  all  models  of 


EQC  - 78068006 


PAGE 

constant  (t  = ot  70  an).  The  relative  height/altitude  of  step  h (all 
linear  dimensions  are  referred  tc  the  width  cf  step)  changed  from 
C-  14  3 to  1.0.  The  relative  length  cf  models  1 changed  frcia  3.43  to 
5.15.  The  half-angle  of  the  wedge  cf  the  ncse  section  of  the  models 
was  equal  to  17c. 

The  model  of  combicaticn  cone  - cylinder,  that  was  being 
intended  for  the  joining  the  results  tc  Jcncwn  data,  had  diameter  of 
the  bpttom  section/shear  D ~ of  70  n,  semiapex  angle  of  cone  of  17° 
and  l - 3.43. 

Models  were  establish/installed  near  the  section/snea r of 
two-dimensional  nozzle  of  5 wind  tunnels  with  the  open  test  section. 
In  the  process  of  experiment,  was  measured  the  base  pressure,  static 
pressmre  and  Hach  number  in  extercal  flew,  The  location  of 
static-pressure  probes  and  the  eleven  receivers  of  base  pressure  is 
shewn  cn  Fig.  1. 

Eage  105. 

Measurements  showed  that  the  static  pressure  on  both  sides  cf 
ocdels  virtually  coincides  and  it  is  close  to  the  static  pressure 
external  flow.  This  indicates  the  absence  cf  the  angle  ot  attack  of 
models  and  that  at  the  selected  length  cf  the  models  of 
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distur tance/perturbatic c frcm  spcut  cf  bottcm  section/snear  they 
virtually  attenuate. 

A change  of  Re  number  (one  ana  a half  tines)  as  a result  cf 
changing  the  length  of  nodel  i.  from  3.43  tc  5.15  with  this  Mach 
number  of  external  flow  did  not  virtually  icfluence  the  value  of  base 
pressure.  Apparently,  and  fcr  the  transiticcal  forms  of  steps  there 
is  a region  of  self-similarity  acccrding  tc  Re  number,  discovered  for 
the  axisymaetric  and  flat/plane  steps  (for  example,  see  [4]). 

Fig.  2,  depicts  the  diagr ai/curves  of  pressure  on  the  end/face 
cf  models  with  Mach  number  = 0.84.  With  the  large  Mach  numbers,  the 
form  of  diagram/curve  is  retained,  changes  ccly  the  level  of  base 
pressure.  For  the  models,  close  in  ferm  to  the  flat/plane  step  (with 
value  cf  h < 0.5),  is  ncted  certain  pressure  increase  on  the  edges  of 
step,  connected,  apparently,  with  the  emergence  of  intense  tip 
vortejes.  With  an  increase  k during  transition  tc  step  with  square 
scctign/shear,  the  increase  cf  base  pressure  in  the  edges  of  step 
becomes  less  noticeable.  In  the  middle  part  of  the  steps,  the  base 
pressure  is  retained  cccstart. 

Subsequently  during  the  analysis  cf  test  results  for  value  /’< 
is  accepted  the  base  pressure,  averaged  on  the  height/altitude  cf 
step.  Fig.  3,  gives  the  averaged  values  of  base  pressure  after 

i i 


DOC  - 73068006 


FAGE  -tT"  . 

X/1 

rectangular  steps  with  diferent  ratios  c£  be igfct/altitude 
Here  corrected  values  o£  tase  pressure  for  ilat/plane  step 
for  ayisyanetric  step  />".  averaged  on  tfce  dcta  cf  works  [2 


to  width. 

p”  and 
] and  [3  ]. 
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Key:  41).  Experiment.  (2).  cone  - cylinder.  (3).  calculation 
according  to  fornula  [1]. 

Page  106. 

The  experiaental  values  of  Ease  pressure,  cfctained  by  the  author 
cn  cylindrical  Model,  ace  close  to  valve  //n  £cr  an  axis ya metric 
step,  which  can  serve  as  indirect  proof  of  tbc  fact  that  the  pylcns, 
which  support  aodels,  were  arrange/ located  at  this  distance  from  the 
tottca  secticq/shear  when  tleir  ellect  cn  tie  value  of  pressure 
virtually  disappeared.  Ihe  values  cf  tie  lase  pressure  cylindrical 
■cdel  and  aodel  with  sguare  bcttci  section/shear  (h  = 1.0)  virtually 
coincided  with  all  Nach  nuaters  of  external  flew. 

Nith  the  decrease  of  tie  relative  heiglt./altitude  of  rectangular 
step,  the  base  pressure  decreases  frea  the  value  of  base  pressure 
after  axisyaaetric  step  tc  the  value  of  base  pressure  after 

flat/plane  step  p",  Moreover  with  an  increase  in  Nach  nunber  of 
external  flow  the  lav  of  this  change  ever  acre  approaches  linear 
(curves  for  Nach  nuabers  - 2.50  and  2.76;  fig.  3). 


The  values  of  base  pressure  after  rectangular  steps  with  the 
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relative  height/altitude  cf  step  h and  Macb  nuiter  ot  external  flow, 
cbtaracd  on  the  proposed  eapirical  dependence 

i>,  p\  (/A1  ,0(i  - h")"'.  (t) 

where 

..Tor ' mu.  ■ 

satisfactorily  will  agree  *ith  experimental  values. 
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18  E S1BUCJURE  OF  POWERFUL  SHOCK  WAVE. 

A.  L-  Stasenko. 

It  is  proposed  to  find  density  profile  in  plane  shock  wave  in 
the  fgtm  of  two  branches  of  smooth  function,  which  asimptotically 
tend  for  density  values  at  infinity  (before  the  wave  and  aftet  it) 
with  tte  dyne  of  the  relaxation  of  the  crdei  of  local  mean  tree  path. 
With  arbitrary  viscosit y-te  xperature  dependence,  the  tasx  is  reduced 
to  quadratures,  and  for  rigid  and  Haxwelliar  molecules  it  is  solved 
in  elementary  functions.  Is  ccpducted  the  comparison  with  the 
results,  obtained  hy  otter  xetbods,  which  stows  that  the  proposed 
examination  is  reasonable  with  Hach  numbers  of  the  incident  flew  on 
tte  order  of  three  and  it  is  above. 

The  proposed  simple  method  of  calculation  of  the  parameters 
according  to  shock  wave  thickness  can  render/stew  useful  for  the 
cough  estimates,  for  example,  during  tte  study  cf  the  passage  of 
small  solid  particles,  driving/mevipg  in  jet,  through  shock  waves  or 
during  the  determination  cf  the  effect  of  tie  strongly  spraying  jet, 
which  escapes  into  evacuated  space,  cn  the  xacrcscopic  bodies  whose 
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size/d i«en sions  arc  coaparahle  with  the  thickness  of  wave. 

let  X-axis  be  perpendicular  tc  shock  icvc,  point  x - 0 
corresponds  to  the  position  of  hydrodynasic 

discoctinuity/interrupticn.  Let  us  search  fer  density  distribution 
r.(x)  in  the  fora  continuous  functions,  with  two  tranches: 

n n, 

n„  — n, 


n. ..  ■ n 

n — n„ 

where 

\ 

there  is  a aean  free  path  cf  molecule  .in  the  systea,  connected  with 
wave;  “t  - mean  free  path  in  gas,  determined  through  the  coefficient 
cf  ductility/toughness/visccsity  p;  ■ - molecular  mass;  u,(c)- 
■acroscopic  and  average/acan  theraal  cf  gas  velocity.  Factor  (u  ♦ <c>, 
/ (fc ^describes  passage  fron  cne  systea  tc  another:  in  the  case  cf  an 
intense  shock  wave,  it  is  ef  the  erder  cf  cuaber  H t in  the  incident 
tlroa  before  wave  <(>-,  “<>'  and  the  erder  cf  cne  - for  a wave 

i r > H, i Indices  by  0,  1,  2 are  related  tc  ccnditions  in  the 
"center"  of  wave,  in  the  flew  before  the  wave  and  after  it  at 
infinity. 


«•  *|» 


exp 


I ,u 

j 'i') 
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. « 0. 


II) 


<i  1 (r) 
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C2) 
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Euctilit y/tou g hnes  s/ viscosity  can  be  tfe  arbitrary  function  of 


teiperature  T. 


Eage  1C8. 


In  relationship/ratics  (1)  are  taker)  iitc  account  tna 
conditions,  by  which  oust  satisfy  function  t(x)  at  infinity  and  in 
zeic:  nth  x ■$--  we  have  n * nl#  with  x • wc  have  n^n* 
(asymptote)  , with  x-^O  we  have  n-+n0  (cc  ct  i 1 1 it  y ) . 


Iroi  (1)  after  d if  fere  1 tiat  ic  n we  obtain 


,/.  / ,r 

Requiring  another  equality  the  derived  both  branches  of  function 
.<>>  £or  a (x>  ] at  poiat  x = C,  we  cbtais  [taking  into  account 


continuity  -.to  and  i(x)] 


1 •»..  ....  whence 


(i  ...» 


Constructed  thus  function  belongs  to  class  C, , since  its  second 
derivative  at  point  x=0  suffers  the  discoatinuity/interr upticn 
(analogous  situation  occurs,  fcr  exaaple,  it  the  case  of  the 
artificial  ductility/tcugboese/visccsity  of  Meuaann-Rikhtaayler,  in 
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which  second  derivative  of  velocity  profile  is  disruptive  on  the 
f cent/ lead ing  and  rear  touncaries  cf  wave). 

I he  thickness  of  wave  from  Pcaqdtl  is  determined  by  the 
expression 

. i . 1 

'■  ,/.  ■,  i (5) 

(t  y M 

With  N l the  order  of  the  tendency  of  values  v2  ana  v„  toward  ' 

unity  is  identical;  consequently,  tie  thickness  L0  limited  by 

relaticnship/ratio  (5),  approaches  final  lint.  This  means  that  j 

re laticnsh ip/ratios  (3)  ct  equivalent  to  them  r elationsn lp/rat ics  O) 

are  wrong  for  weak  waves.  * 

Iroa  the  law  of  conservation  cf  nas$,  we  have 

un  U,ll,  II  II  (,,(  , 

The  law  af  conservation  cf  energy  let  is  accept  in  the  form, 
analogous  to  Bernoulli's  integral  in  the  thccry  cf  the  steady  flow: 

II  ^ H j II 

<V  ' ; l 'i  t n fpT,  } a . (7) 

where  r„  - the  haat  capacity  at  a constant  pressure. 

This  relationship/ratic  is  fulfilled  at  any  point  of  wave  only 
in  the  case  of  Prandtl's  nunber  Pr  - 3/ U;  tie  latter  is  most  close  to 
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the  value  o f number  Pc  foe  twe-  and  triatooic  iclecules.  Is 
real/actual,  according  to  Eucken's  fcrsula,  nuater  t'r  4 * where 

- relation  heat  capacity  c£  gas,  sc  that  ratio  4Pr/3  is  equal  £ol 
such  aclecules  56/57  and  64/63  respectively,  for  monatomic  gases  Pr  = 
2/3  and  the  relation  indicated  is  equal  to  6/9. 

Utilizing  relationship/ratics  (6),  (7)  apd  taking  into  account 

that 

* i 

2 ", 

I • ' , ' Mf  ' 

let  us  knowingly  satisfy  all  laws  cf  ccpservaticn  at  infinity  before 
and  after  wave. 

Cage  1C9. 


Iroa  relationship/ratics  (2),  (6)  and  (7)  we  have 


'•  1 , I •’ 
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» * 


<") 


•M, 


M: 


C») 


where  1,  - mean  free  path  in  the  gas  hetore  the  wave.  Sunstituting 
(6)  in  (3),  we  will  obtain  ordinary  differential  equations  for  two 
tranches  of  the  functicn  cf  reverse/inverse  te  •to.  solved  in 


quadratures 
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It  least,  for  two  cases  - rig  id  < i T\  and  Maxwellian  (>j  ~ T) 

■ clecales  - quadratures  for  ,(,)  can  te  expressed  in  the  elementary 
functions: 


for  rigid  aolecules,  »/h  - \ '< 


X 


"3“ 

i 

|| 

<3* 

/,f 

, i : o.  •- 

a - j/  | - n- 
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"2 

v2 

(10) 


for  Maxwellian  molecules,  p/f-i  = 6, 
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Big.  1 and  2,  give  the  airfoil/prcfiles  of  shock  waves  tor  the 
case  of  rigid  molecules,  designed  crj  formulas  (1C),  and  the 
airfoil/profiles,  designed  it  work  [1]  acccnding  to  Na  viec-Stokes 
equations,  the  method  of  flctte-Snith  and  Mcrte-Carlo.  In  Fig.  3,  are 
constructed  the  density  profiles  and  teapeictur es  tor  the  case  of 
Kaiwelliaq  molecules  on  fcrsulas  (11)  c£  present  article  and  the 
a ir toil/ prof iles,  designed  in  work  [2]  acocrdrng  to  Na vrer -S tokes 
equations,  flotte-S mitn *s  method  and  with  tfc<  aic  ot  elliptical 
(two-temperature)  the  distribution  functions  of  molecules  according 
tc  speeds. 
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Fig.  1. 


Fig.  2. 


Fig.  1. 


Fey:  {1).  Formulas  (10).  (2).  Nav  ier-Stckes . (2).  Mott  e-  om  it  h.  (4). 

flcnte  Carle. 


Fig.  2. 

Key:  i 1).  Formulas  (10).  (2).  Nav ier-Stckes  . (2).  Hotte-  jmith.  (4). 

fente  Carle. 

tage  110. 


From  the  given  curve/graphs  it  is  evident  that  for  the  powerful 
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shock  waves  (N  ^3)  of  the  difference  between  the  density  profiles 
and  temperatures,  obtained  ty  the  enumerate c methods  (including 
pccposed  in  this  article)  cf  cne  order.  Furthermore,  is  noticeable 
powerful  asymmetry  (relative  to  the  "center"  of  wave)  obtained 
simple,  is  noticeable  powerful  asymmetry  (relative  to  tne  "center"  cf 
wave)  cf  the  obtained  density  profile  which  was  noted  both  in 
thecretical  and  in  experimental  [3],  [4]  investigations,  in 
particular.  Fig.  4 gives  experimental  data  cn  the  measurement  cf 
density  distribution  in  sbcck  wave  by  the  .method  of  electron  beam 
[3].  Shock  wave  wis  obtained  in  air  (*  M;M,  id  am  Hg)  during  the 

flew  around  disk  and  sphere  by  diameter  1.5«1Q'2  ■ (white  and  LlacK 
small  circles  respectively).  In  Fig.  4,  are  plotted  also  the  results 
cf  the  calculation  of  the  a irf cil/pccf ile  cl  wave,  carried  out  by  the 
author  of  work  [3]  according  tc  Mctte-Smith 's  method  (dot-dash  line). 
Unbroken  curve  - calculation  according  to  formulas  (11)  for  the 
i‘->xwellian  molecules:  in  the  range  of  temperatures  T, ::  35°K  <T<2  fcb°K  = 1 ° 
in  flpw  the  dependence  cf  the  ductility/toucbness/viscosit y of  air  cn 
temperature  is  close  tc  straight  line  p ~ 1,  moreover  • 5 • 1 C” 6 

kg  «oi -*  •£  ”* , whence  mean  free  path  in  gas  is  ettcined  equal  to 
iiXO. 37» 10' 3 to  m.  Fig.  4,  shows  that  the  ccnvergence  experimental 
data  and  theoretical  results,  obtained  cn  tie  tasis  of  the  proposed 
method,  sufficiently  goed. 
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BALLISTIC  TUBE  FOR  0RA6  MEASUREMENT 
ON  MODELS  IN  FREE  FLIGHT  AT 
HYPERSONIC  SPEEDS 


L.  P.  Our'yashkin,  A.  P. 
Krasil 'shchikov,  and  V.  P. 
Pod ob in 


A brief  description  is  given  of  a ballistic 
wind  tunnel,  the  operating  principle  of  which  is 
based  on  the  firing  of  models  towards  the  super- 
sonic flow  in  the  wind-tunnel  working  section. 

The  installation  is  Intended  for  the  measurement 
of  the  drag  coefficient  and  for  studying  flow 
fields  of  axially  symmetrical  bodies  in  the 
range  of  supersonic  and  hypersonic  flight 
velocities . 

The  range  values  of  the  M number  in  the  ballistic  wind  tunnel 
is  equal  to  1.5-15.  High  resulting  M numbers  are  obtained,  in 
the  first  place,  as  a result  of  the  firing  or  models  towards  the 
flow,  and,  in  the  second  place,  as  a result  of  the  lowering  in  the 
speed  of  sound  in  the  working  section  of  the  supersonic  wind  tunnel 
by  means  of  cooling  of  the  air  in  transit  through  the  nozzle. 

For  the  resulting  M number  the  following  expression  is  correct: 

y 

M-M.+ 

where  V - the  speed  of  flight  of  the  model  relative  to  earth; 

a - speed  of  sound  in  the  flow; 

M - number  of  M flow. 
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Th«  values  of  the  Reynolds  number  are  obtained  considerably 
higher  than  In  wind  tunnels  at  the  same  M number.  The  expression 
for  the  Reynolds  number  can  be  written  In  the  form 

R, 


where  d - the  diameter  of  the  midsection  of  the  model; 
Vn  - flow  velocity; 
p - air  density  In  the  flow; 

W - viscosity  of  air  In  the  flow. 


The  first  term  corresponds  to  the  Re  number  of  the  model  at 

zero  velocity  of  the  firing,  and  the  second  term  - to  the  Re 

number  during  motion  relative  to  earth.  The  Re  numbers  realized 

In  the  ballistic  wind  tunnel  under  various  conditions  of  the 

6 7 

experimentation  are  of  the  order  of  10  -10  . 


In  ballistic  Installations  the  stagnation  temperature  Increases 
with  an  Increase  In  speed  of  the  flight  of  the  model  [1].  The 
maximum  stagnation  temperature  in  the  described  installation  Is 
equal  to  ~2800°K. 

Figure  1 gives  a diagram  of  a ballistic  wind  tunnel  which 
consists  of  three  basic  elements:  the  wind  channel,  rifle  stand 
and  electron  optical  equipment. 
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providing  ror  the  compensation  of  the  boundary  layer,  the  nozzles 
have  an  octahedral  cross  section.  The  Inlet  section  of  the 
working  section  3,  whose  length  Is  2^  gauges,  has  the  shape  of 
a correct  octahedron  with  the  diameter  of  an  Inscribed  circle  of 
7**.5  nun.  In  proportion  to  the  distance  from  the  Inlet,  the  cross 
section  of  the  working  section  Increases  because  of  a decrease 
In  the  area  of  the  angular  Inserts,  which  accomplish  a Junction 
from  the  octahedral  cross  section  to  the  square.  In  this  case  two 
pairs  of  opposite  walls  of  the  working  section  of  the  wind  tunnel 
remain  In  parallel  to  each  other.  Installed  on  these  walls  are 
optical  windows  for  photographing  the  model.  The  working  section 
Is  closed  by  a subsonic  diffuser  4 which  Is  connected  with  a 
turning  elbow  5. 


The  braking  of  model  after  the  flight  of  the  working  section 
occurs  In  the  high-pressure  air  In  the  phase  of  trajectory  L,  which 
consists  of  the  Internal  cavity  of  the  antechamber  and  tube  6 with 
the  collector  of  the  model  7 deflected  90°.  The  velocity  of  the 
model  In  this  section  decreases  approximately  5-10  times. 


The  models  were  fired  from  smooth  and  threaded  barrels  8 of 
the  caliber  1^.5  mm  with  flight  velocities  of  500-2000  m/s.  The 
firing  was  conducted  with  Dural  steel  or  brass  models. 


For  the  test  work  In  the  range  of  the  low  supersonic  m numoei 
.erodynamlc  circuit  of  wind  tunnel,  i.e.,  the  nozzle,  working 
on  and  diffuser,  was  dismantled  and  In  Its  place  a thermal 
er  was  Installed  where  the  pressure  could  change  from  1 to 
m.  The  Inlet  of  the  model  Into  the  thermal  chamber  with 
iased  pressure  was  accomplished  with  the  help  of  an  explosive 
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Th*  ballistic  wind  tunnel  was  equipped  with  three  Identical 
measuring  stations  for  photographing  the  model  and  measuring 
the  time  between  the  moments  of  photographing.  Each  station 
Is  Shadow  system  with  a parallel  light  beam. 

In  a direction  perpendicular  to  the  optical  axis  of  the 
shadow  system,  In  the  center  of  the  photographing  field  there 
passes  the  plane  of  photoblocking,  which  Is  a slotted  parallel 
light  beam  which  goes  from  source  9 to  the  photoelectric  pickup  10. 
To  prevent  the  Illumination  of  the  photographic  film  by  diffuse 
scattered  light,  the  source  Is  equipped  with  a light  filter  which 
displaces  the  spectrum  of  light  Into  the  region  Insensitive  for 
photographic  film.  The  model,  flying  through  tho  working  section, 
consecutively  Intersects  the  light  beams  of  the  photoblockings  of 
three  stations.  The  signal,  which  appears  In  the  photosensitive 
device  at  tne  Intersection  of  the  light  ray  by  the  model,  enters 
Into  the  control  unit  c the  spark  light  source  11.  There  occurs 
an  Intense  light  flash,  and  the  model  Is  fixed  on  the  film  of  the 
camera  of  the  first  station.  Simultaneously  the  impulse  of  the 
light  of  the  spark  source  falls  on  photoelectric  head  I".  12,  which 
co  isists  of  the  vacuum  phototube  STsV-'l  and  cathode  follower. 

The  electrical  signal  of  the  photoelectric  head  starts  the  first 
electronic  chronometer  13.  The  second  station  works  similar  to 
the  first  with  the  only  difference  being  that  simultaneously  there 
occurs  the  stopping  of  the  first  chronometer  and  the  starting  of 
the  second  one,  which  Is  stopped  from  the  flash  of  the  light  source 
at  the  third  station. 


The  thus  obtained  space-time  dependence  of  the  flight  of  the 
model  was  used  for  the  calculation  of  the  drag  coefficient  [2]. 
The  accuracy  of  the  measurement  of  time  by  the  electronic 
chronometer  Is  equal  to  0.25,10-f>  s,  and  the  accuracy  of  the 
measurement  of  the  position  of  the  model  on  the  trajectory  Is 
equal  to  0.2  mm. 
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j Given  below  for  an  example  are  some  experimental  reaulta 
obtained  on  the  Installation  desoribed. 


Plgure  2 depicts  depending  on  the  M number  the  values  of  the 
drag  coefficient  of  the  axlsymmetrlc  model  with  the  conical 
expanding  panel  half-angle  of  aperture  of  which  is  equal  to  35°  • 
The  head  section  of  the  model  represents  a cone  with  the  half- 
angle of  aperture  equal  to  15° • The  length  of  the  cylindrical 
part  of  the  model  is  equal  to  1.5d,  and  the  ratio  of  the  diameter 
of  the  cylindrical  part  to  the  diameter  of  the  midsection  of 
the  model  d/D  * 0.565.  A decrease  in  the  drag  coefficient  in  the 
range  of  M numbers  of  1.5  to  -8  are  caused  by  the  change  in  the 
pattern  of  flow  of  the  model.  With  a further  increase  in  M number 
the  drag  coefficient  does  not  change. 


Plgure  3 shows  the  changes  in  the  drag  coefficient  of  the 
ellipsoids  of  revolution  at  the  fixed  M numbers  depending  on  the 
ratio  of  the  semiaxis.  (a  - horizontal  semlaxls,  b - vertical 
semiaxis).  An  equidistant  displacement  of  this  dependence  with  a 
change  in  the  M number  is  observed.  — i — .m/  . — i > — j 
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! Th«  development  of  the  width  of  the  turbulent  nucleus  of  the 
wak^  behind  the  sphere  of  up  to  3500  gauges  In  air  at  ataospherlc 
pressure  Is  shown  on  Fig.  it  In  the  form  of  dependences  of  7 on  x 
whehe  7 ■ 6/d  Is  the  width  of  the  turbulent  nucleus  referred  to 
In  the  diameter  of  the  midsection  of  the  model,  and  x • d/x  Is 
the  distance  from  the  model  to  the  place  of  the  measurement  of 
the  width  of  the  wake  referred  to  the  diameter  of  the  midsection 
of  the  model.  The  turbulent  wake  behind  the  sphere  In  Its  near 

C 

part  at  m ■ 2.8  and  Re  ■ 10  Is  shown  on  Fig.  5.  For  a comparison 
Fig.  i<  shows  test  data  from  works  [3]  - [5]  and  the  theoretical 
dependence  [6]  for  M ■ 8.5.  Despite  the  great  difference  In 
velocities  all  the  test  data  agree  well  with  each  other.  In  the 
Interval  of  the  gauges  from  30  to  3500  the  development  of  the 
width  of  the  turbulent  wake  can  be  approximated  by  empirical 
relation  7 » 0 . 208ilx^  ^ . 
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Fig.  4. 

KEY:  (1)  Data  of  authors. 
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lage  115. 

THE  EtEECT  CF  THE  EROCESS  Cl  CVEPCEAEGING  C N IhE  EFFECTIVENESS  OF 
1C  N 1C  SOURCE  WITH  VOLUHE1RIC  ICNIZA11CN. 

Yu.  E.  Kuznetsov,  V.  P - Rudckcv. 

Is  examined  the  flat/plane  model  cf  icric  source  wua  the 
distributed  over  its  length  neutral  coi^onut  cl  plasma. 

It  is  shown,  that  the  process  cl  overcharging  whose  intensity  i 
otter lined  by  the  physical  properties  cf  gar  (atomic  weight,  the 
section  of  overcharging  and  lcnizatico)  and  by  electron  temperature, 
limits  the  value  of  the  overall  efficiency  cf  the  use  of  a mass. 

The  effectiveness  of  tfce  work  of  icnic  source  is  characterized 
ty  the  value  of  the  coefficient  of  the  tse  cf  mass  n (ratio  of  ion 
flew  in  by  output/yield  section  to  the  contraption  of  the  work 
substance,  expressed  in  unity  of  eguivalent  icn  current). 
Comparatively  high  effectiveness  have  mcurces  *ith  the  ionization  of 
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neutral  particles  by  electron  collision  [1]  - [5].  The  maximum 
effectiveness  of  such  sources  is  examined  it  work  [6],  in  this  case 
the  process  of  overcharging  in  the  volume  cf  discharge  is  not 
examined. 

BnliJse  the  conditions  cf  the  accelerator  cf  the  plasma  where  the 
process  of  the  distributed  overcharging  plajs  positive  role  [7], 
under  conditicns  of  ionic  source  this  process  can  become 
barrier/obstacle  for  achievement  of  the  hich  values  cf  tne 
coefficient  cf  the  use  cf  a mass,  since  the  rapid  neutral  particles, 
formed  as  a result  of  overcharging,  have  the  lew  ionization 
prchabilit  y. 

Jn  the  present  woi;k  are  approximately  examined  the  processes, 
which  cccur  in  the  camera/chamter  cf  iebizatien,  taking  into  account 
overcharging. 

The  camera/chamber  cf  icnizaticn  is  cctsidered  plane.  The 
neutral  and  ionic  components  of  plasma  move  alcrg  coordinate  x, 
perpendicular  to  chamber  walls  which  are  found  under  the  potential, 
negat4ve  with  respect  to  the  potential  cf  plasma.  Is  examined  the 
case,  when  the  potential  difference  near  the  electrode  is  greater 

than  the  electron  temperature,  expressed  in  electron  volts.  Through 

<x, 

the  l«ft  wall  (Fig.  1ft)  cccus  the  irleakage  of  working 
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■ediua/propellant.  The  density  of  the  iccoiia;  flow  of  work  sutstance 
/„.  jest  as  the  densities  ct  the  flews  cf  ceitral  particles  in  the 
ca mera/cha aoer , it  will  he  expressed  in  unity  ct  equivalent  ion 
current,  i.e.,  in  aaperes  tc  square  centimeter. 
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Fig.  1. 


Eage  116. 


The  discharge  of  icnic  and  neutral  coaioqevt  occurs  through  the 
right  wall  - to  riding-crop  with  transparency  t s' , where  - area 

cf  cpening/apertur es,  vK  - ccBmon/g«n«ral/tota 1 area  of  chamber  end. 
Fcr  ions  and  geutral  particles,  the  effective  transparency  of  grid  is 
accepted  identical  and  equal  gecoetric  transparency  y. 

Neutral  component  let  cs  consider  consisticg  of  three  parts 
(Fig.  1b): 

1)  the  flow  of  "cold"  neutral  particles  with  a density  of 
;,n(‘).  dri  wing/mov  ing  of  left  wall  tc  the  right  (always  positive). 
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2)  the  flow  of  "cold"  neutral  particles  wit!  a density  of  /«<*)• 
dr iviog/moving  from  right  wall  tc  the  left  (always  positive). 

3J  the  flow  cf  rapid  nertral  particles  with  a density  of 
termed  as  a result  of  overcharging. 

The  first  two  unidirectional  flows  are  characterized  Ly  certain 
average  constant  thermal  velocity  i„.  deteraiged  by  the  temperature 
cf  walls  the  third  - by  velocity  whose  cverage  value  the  order  or 

the  velocity  of  ions  fi  is  determined  ty  electron  temperature. 

High  velocity  and  the  light  pcrticn  of  the  flow  of  rapid  neutral 
particles  in  comparison  with,  velocity  and  flew  cf  cold  neutral 
particles  makes  it  possible  to  disregard  in  eguations  tue  terms, 
which  consider  secondary  overcharging  and  the  icnization  of  rapid 
neutral  particles. 

let  us  consider  that  the  ions  move  from  certain  section  x0 
within  the  camera/chamber  to  the  left  and  right  walls  (as  this 
follows  from  work  [8]),  moreover  on  both  walls  ion  currents  are  egual 
anc  are  determined  by  Echa's  fotiula  [9]. 
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The  neutral  particles,  which  were  being  formed  because  o i 
overcharging,  move  to  the  side  of  the  acticr  of  ionic  flows,  i.e«,  cii 
the  left  side  of  the  ca aera/chaabe r tc  the  left,  and  in  the  right 
side  gf  the  camera/chamter  to  the  right  (see  f i g.  1b).  la  section  x,, 
the  flow  of  these  neutral  particles  is  eguel  tc  zero. 

Are  accepted  the  following  assumptions:  ir  any  section  x,  all 
ions  have  identical  velocity  t.(u>  which  depends  only  on  x,  moreover 

Jl<x>-  the  mean  free  path  of  neutral  particle  before  collision 

•'flr 

with  neutral  particle  and  with  ion  is  such  acre  the  distance  between 
the  walls  of  source;  electrcnic  cc ncentrat icn  and  the  distribution 
function  of  electrodes  along  X-axis  are  constant. 

Within  the  framework  of  the  adoptee  assumptions,  the  ion  flews 
and  rapid  neutral  particles,  through  tie  arbitrary  section  x taking 
place*  are  determined  cily  by  their  generation  and  are  expressed 
respectively  by  the  following  dependences; 

ilj,( x)  i n,,  n„  i,  tv  dx\  (I) 

en* 11  n 3n«|i  I Vj  (a)  | dx,  (2) 

where  »n  - a concentration  of  cold  neutral  particles, 

jn\  (*)  ‘ j m ( * ) 

rV„ 

n„„|  - the  section  of  overcharging*  accepted  by  constant  in 
value  in  each  section  x 1 ; 
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«i*v  - averaged  according  to  distribution  function  product  of 
ionization  cross  section  to  the  velocity  of  electrons. 

FCC1NC1E  *.  According  tc  the  estimation  of  the  authors,  introduced  by 
this  assumption  error  into  the  flow  value  cl  the  recharged  neutLal 
particles,  for  the  «ajoiity  cf  work  substances  and  in  the  real  speed 
range  cf  iens  does  not  exceed  20o/c  as  a result  of  the  icgarithnic 
dependence  of  the  section  cf  symmetrical  resonance  overcharging  cn 
the  mcdule/mcdulus  cf  relative  velocity.  E M. FC C 1 NOTE. 


the  flows  of  cold  neutral  particles,  ori vi tg/movi ng  from  left  to 
right  and  froa  right  tc  left,  ate  expressed  by  the  respectively  given 
belcw  dependences  and  represent  by  themselves  the  loss/depreciation 
cf  neutral  particles  as  a result  cf  ionization  and  overenarging: 

i<-»)-  - en,  n„  a,  r,  etx  en,  n„  9„t.p  (3) 

ilj'ni  <-*  ) en,  «/,  -i  v,  it\  -I  en,  n„  a,,,.,, ; v,  (*)  <t*  M ) 

Eage  1 17. 


Xguations  ( 1)  — ( 4 ) for  determining  the  ccnstants  it  is  necessary 


tc  supplement  by  the  following  boundary  conciticns: 
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in  ~ ini  (<•)  1 in.  (")  ; ii  (")  i j„  (0) 

(-r 

0); 

(•r>) 

Jn  - j„i  (V,)  i )„  (-tr()  (x 

fi): 

<r>) 

in-  (-fj)  ^ ( 1 ;)  1 (•»)  i ii(  y-i)  , j'n 

(->:•)! 

(') 

7 , < «•.. ) 

1 Jl  l->  ) ! /„;(>.,) 

(«) 

The  physical  sense  of  expressions  |5)-|7)  is  clear:  the  total 
flew  cf  all  particles  is  equal  to  ccnsuapticn,  and  in  section  xlf  the 
icq  flews  and  rapid  neutral  particles  are  equal  to  zero  [see  (5)  and 
(6)].  from  grid  in  section  i2  the  ceutral  particles  reflected  and 
neutralized  ions  give  ccunterflow  of  ccld  ceutral  particles  (7). 
Condition  (8)  determines  the  coefficient  cf  the  use  of  a lass. 

He  will  be  restricted  to  the  exaninaticn  cf  the  Maxwellian 
distribution  function  of  electrons.  Let  us  intrcduce  the 
designations: 

a,  l'r°j  - ionization  cross  secticn,  averaged  according  to  the 

tv 

Faxwelliaq  distribution  function  cl  electrcps. 

v,  - the  aean  arithmetic  velocity  of  electrons. 

/,(<))  J,Uj)  - Bcha's  foraula  [ 9 ]r  ir  which  t i 0,2,  | r"r 
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5V? 

(where  m'  - ratio  of  the  wess  of  electron  to  the  mass  o t ion)  . 

Mi 

a l'a"ep  - the  parameter  of  ove  rc  bar  g ing. 

*•3, 

Under  flews  ynl  (*),  jni(x),  J'„  u>.  /,(*•)  let  us  subsequently  understand 

their  relation  to  the  module/modulus  of  the  raximum  strength  of  icr 
current  A-  while  under  independent  the  variable  x let  us  understand 
the  d iwensicnless  quantity,  which  represents  by  itself  tne  ratic  of 
longitudinal  coordinate  tc  the  mean  free  path  cf  neutral  particle  to 
the  ionization: 

. Vn 

>i  * 

n„  v,.  i, 

Latter  makes  sense  when  the  effectiveness  cl  icrization  does  net 
depend  on  coordinate. 

In  the  adopted  designations  we  cbtain  the  following  system  of 


differential  equations: 

dji(x)  | Jm  (.V)  y„(.r)|  </.«;  <■•) 

d/„(x)  t j,  ( <)  | /„,  ( v)  y„  (a)|  tlx,  (Hi) 

dj„ i(  v) jnl  (r)  1 1 A | y,  (X)  | | dx;  (II) 

<ij„i  ( x)  - y„; (x)  1 1 -j  i ; y,  (•*)  1 1 < |2) 

with  the  boundary  conditions 
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x ^2 


J'n  ■ J„l  ((>)  i in  (0)  ; 4,(0) 

- 1 (x  0); 

(i;t) 

in  x ini  * 'il  in?(xi) 

(«■=  v,); 

(M) 

iiri(x»)  =-  ( 1 7)  I in  (X,)  1 

i„(x j)|  (v  » i, 

(IS) 

1 

j'nl  (X  .)  1 - j'n<  i,| 

(V  Xj) 

(16) 

Page  118. 


Solving  together  (S)  and  (10),  let  as  have 


Key:  i 1)  . with. 


j'„ (<)  1/'2l')  »pn  >-<•')" 

j'n  (X)  --  — AI‘^X)  |®H  j,  (X)  < 0 


(17) 


Then  is  deterained  the  value  cf  the  current  of  rapid  neutral 
particles  cn  the  boundaries: 


.1 


H 0) 


(IK) 


Let  us  replace  of  the  differential  of  independent  variable  in 
equations  (11)  and  (12)  fay  the  expression*  obtained  fron  (9), 


llx 


'III  ( « I 


let  us  introduce  the  designations 

i»\  ( «)  + >»•(-<)  i’ 

Jm  ( O ! *■ 


(19) 


(20) 
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he  Mill  obtain  the  system  of  the  differential  equations 

dy  ■ 1 1 I /»l>,(t)  | <!/,  ( »);  | 

»•  (2)) 
ji  - ■ 1 1 i t y,< o'  i Jjt(x)  j 

Its  solution  taking  into  account  boundary  conditions  and  (18) 
takes  the  fora: 


in  field  j,  u)  .0 

V in  y,  (ft 

l K'\  *Jn  7, <0+0  f Ah,)  j‘  ( » ) + Ajf{  t)  t — | 


Aj‘{x) 

2 


t-y,'(.o 


field  J,  (>)<<) 


V in  i,  < * > 


Aj]  ( i) 


where 


A- )Ux) 

| 2i'  | ( 2/„/,l  k)  i (I  + Aj„)j{(y)  Ajt(x)  4 


(22) 


(23) 


Euncticn  z,  proportional  to  the  ocncee  tr  at  ion  of  cold  neutral 
particles,  represents  by  itself  tbeir  relative  distribution  along  the 
length  of  the  caaera/chaaf er  cf  ionization. 


1 he  fourth  integral  of  systems  (S)-(1i)  Mill  be  located  during 
the  solution  of  differential  equation  (S)  , written  in  tne  form 

, <Ui  (XI 


(24) 
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which  it  is  necessary  separately  tc  integrate  in  two  fields  cf  the 
caaera/chaaber  where  v<0>  >«  and 

Page  119. 


Relative  value  of  the  consumption  cf  verkisg  aedium/prcpellant 
/„  in  (22)  and  (23)  can  fce  expressed  by  dependence  /„-  which 

i 

easuee  froa  the  determination  of  the  coefficient  of  the  use  of  a 
working  mediuu/propellant. 


The  results  of  numerical  count  for  two  different  values  and  in 
different  parameters  of  overcharging  A are  given  to  Fig.  I and  3. 

A system  of  differential  equations  (9)- (12)  can  be  utilized  for 
the  case  of  the  overlap  cf  working  ned ium/ p t epe llant  from  the  side  of 
gr ids^ 


In  this  case,  will  chatqe  only  boundary  conditions: 

y'm(O)  -r/n-m  i I +•  U ");  (25) 

r, );  (26, 


/„<*,)  I /»  (fi)  (> 

A 

/»»(■*  ) (>  - ;<)  )n\  (M  ( I t-  ,, 

A ; L i 

> T'"  t i 

/„!  (<•..)  I I f 

Respectively  will  cbaqce  the  integrals  of  the  system; 


t in  (' 
( v > X..) 


(27) 

<Z«> 


*n  field  M*>># 
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M*)- 


I 


/ 

/I 

1/  /?(-«■) 

1 ' 9 Ji(*)  '' 

4o  field  h (■*)<« 


vhexe 


J\U) 


* = 1 ^ 


|/  a*-1*  j 1 * 2 y,(>) 

O,  2/"  / y"  , 

V V •( 


+ 2C»  • 

/«  i 


f 
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cn  th«  boundary  of  [lasia  (in  the  examined  t ne-dimensional  case)  is 
always  equal  to  value  A/2. 


Recall  that  the  constant  A defends  only  cr  tne  physical 
properties  of  gas  (atomic  weight,  the  secticn  of  overcharging  and 
ionization)  and  of  electron  temperature  and,  that  very  substantially, 
directly  it  does  net  deperd  cn  the  value  cf  the  camera/cha mber . 

Bren  (It)  and  (28)  it  is  evident  that  the  upper  limit  of  the 
value  cf  the  coefficient  cf  the  use  cf  a mars  even  for  the 
sufficiently  long  camera/chambers  is  the  expression 
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>$4 


Fig.  4. 
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1HE  VEBSION  CF  DISPERSICB  HETfiOD  A N L HE  TASKS  CF  BOUNDARY  LAYER 
ST  A E IE  IT  Y. 

V.  fl.  Lutovinov. 

Is  examined  the  version  of  dispersion  aethcd  for  the  systems  of 
the  differential  second  crder  eqaaticns.  Is  noted  its 

ccamunication/connecticn  with  the  task  cf  tie  factorization  of  linear 
differential  expression,  is  shown  tie  pcssifility  of  applying  the 
method  for  the  numerical  solution  cf  the  tasks  cf  boundary  layer 
stability  within  the  framework  of  linear  theory. 

let  us  examine  the  solution  of  hcundar y-va  lue  problem  for  the 
system  of  the  differential  equations  of  the  second  order 

>1"  A M*  /'  (U 

let  on  one  of  the  ends  of  the  cut  cf  integration  ab,  for 
example,  with  y - b the  boundary  conditions  be  arbitrary  relatively 
‘I  and  >r.  tmt  on  other  - »ith  y = a they  take  the  form 

i*  i km  ( ; 


here 
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'1 , 

h 

h 

. 

,1  ■ 

, 

• ; 'l  ’ 

'I ; 

7 

1 

1 • 

dy 

ii 

1 « 

| l 

, f" 

• ‘‘in 

* n • 

*i« 

] 

: 1 ; 

A 

• 

Ii"’"'  • 

* (inn 

"nX  ' 

‘ "(in 

an  and  - conplex  cunters  (i. 


1»  . • • # c)  • 


On  the  basis  of  the  ideas  of  the  dispeisicc  method  of  Gelfand 
and  Lckutsiyevskiy  [1],  let  us  search  fer  octrix/die  B = B (y)  and 
vector  'i>  *i»  ( v)  • where 

1 ■ *•-  I*' 

fl  ; • ''' 

I'll  ' * ' ,f 

which  satisfy 

>r  * II  ( 'I  <l>;  (>) 

/Ihi)  K,  ■!•(«)  ( (3) 

Pace  122. 

■ if ferenbiating  (2)  and  utilizing  (1)  and  (2),  we  will  obtain 
'i"  \ Oi  /r>v  . ■!»  — H‘i>  since  •»’  is  soluticn  (1)  and  satisfies  (1*), 
aatrix/die  B vector  functicr.  «t*  they  erst  satisfy 
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> 

It  It  t.  (4) 

/M*  = / , (5) 

lUni  K.  ( « ) (.  (6) 

Integrating  systems  (k),  (5)  kith  initial  conditions  (6),  we 

bill  attain  at  point  y = fc  of  conditions,  wlich  are  missing  tor  the 
investigation  of  the  solvability  of  boundary- va lue  problem.  Alter 
determining,  if  this  is  possible,  '<</').  let  is  find  via 
re  verse/ in  verse  screw  die  solving  the  problem  of  Cauchy  for 

egtaticn  { 2)  . 

It  is  possible  to  show  that 


here  n . Cf  , there  is  A of  the  linearly  independent  solutions 

-*  -»  -?  -> 

which  satisfy  *+K#=C  with  y = a;  0 1 - mataix/die, 

reverse/inverse  U.  ' 

Ihe  described  method,  as  other  dispersion  methods*  can  be 
obtained  from  the  theorem  alout  the  factorization  of  linear 
differential  expression  [2].  From  this  same  theorem  follows  the 
ccrtinuous  dif fere ntiat i lit y of  matrix  elements  8. 

During  the  study  of  bctndary  layer  stability  within  the 
framework  of  linear  theory,  it  is  recessary  for  the  equation  of 

Cr  r-Sf mmerf eld 

f 

s.,v  2i*  <f"  I-  «‘<f  / <iAP |(i»  — <■)(*"  a* <f)  — v"  f I <K> 
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kith  the  boundary  conditions 


'f 

w 

T 


«<f"  I av)  0; 

!*?"-•*(?’  + h)  - « 

<£> 

f -f  0 up  ii 


I t* 

I II  pH  v-  '<! 

J 

v — 0 


Key:  41).  with 


CJ) 

(10) 


to  determine  among  eigenvalues  R,  «,  , c,  : /<•,  cf  boundary-  va lue 
problem  (8)  — (10)  the  field  where  there  exist  «,  <>. 

Equation  (8)  and  boundary  conditions  4 S ) can  be  presented  in  tne 

fen 


he  re 


Faqe  123. 


if"  .1*1  0; 


HI) 


if-  a;  *r  — o * it  h y s; 


(i-') 


A 


— hUv"  i~  | hk(v  f)  | 

I 1 I! 

K - 1 » 

o «i 


The  dispersion  method  in  questicr  is  led  in  this  case  to  the 


rstegsation  of  the  system 
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with  the  ipitial  conditions 

/(('.)  K '"I 

The  investigation  cf  a guesticq  concerting  the 
a £ fil  Jat  ion/accesscry  cf  the  values  c£  (diaieteis  H,  a,  c,  with  the 
number  cf  its  own  is  reduced  taking  intc  acccuct  (10)  to  checking  cf 
rendition  bt2(0)=0.  The  presence  of  com  nun i c at ic n/conn ec tion  (7) 
gives  grounds  to  rely  c(  the  stability  cf  the  process  ot  integration 
(13).  This  is  confined  ty  numerical  experiments. 

let  us  note  that  the  method  in  guesticr  can  be  used  alsc  tc  the 
stability  analysis  cf  jets. 
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